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FOREWORD 


The work described herein was performed at the Aerojet Liquid Rocket 
Company under NASA Contract NAS 8-33553, with Mr. Fred Braam, NASA-Marshall 
Space Flight Center, as Project Manager. The ALRC Program Manager was 
Mr. Jeffery W. Salmon, and the Project Engineer was Mr. Gregory M. Meagher. 

The technical period of performance for this study was from 9 July 1979 
to 31 October 1980. The overall objective of this study was to upgrade per- 
formance and heat transfer analysis techniques for the dual-throat and dual- 
expander combustion chamber concepts. 

The author wishes to acknov/ledge the efforts of the following personnel 
for their significant contributions; 

Dick Ewen, ALRC Gary Nickerson, SEA Inc. 

Jerry Pieper, ALRC Irwin Alber, Consultant 

Jim Duey, ALRC Joe Hoffman, Purdue University 

This report presents the Summary, Section I - Introduction, Section 
II - Conclusions and Recommendations, and Sections III, IV, V and VI - 
Technical Discussion. 
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SUMMARY 


A. STUDY OBJECTIVES AND SCOPE 

The objective of this study was to expand and extend the analysis 
models and parametric studies previously performed for both the dual -throat 
and dual-expander engine concepts. This was accomplished within five basic 
tasks, including (1) improvements to the existing dual-throat aerodynamic and 
performance prediction computer model, (2) preliminary geometric analysis of 
the dual-expander concept, (3) preliminary flow field analyses of the dual- 
expander concept, (4) further preliminary heat transfer analysis of both con- 
cepts, and (5) engineering analysis of data from the NASA/MSFC hot-fire 
testing of a dual-throat thruster model thrust chamber assembly. As a result 
of this effort, a quantitative comparison of the aerodynamic, performance, and 
thermal design characteristics between the two dual -nozzle concepts can be 
derived. 


B. STUDY RESULTS 

During this study, analytical models to predict performance and 
operating characteristics of dual -nozzle concepts have been developed and 
improved. Aerodynamic models are now available to define flow characteristics 
and bleed requirements for both dual-throat and dual-expander concepts. 
Advanced analytical techniques were utlized to provide quantitative estimates 
of the bleed flow, boundary layer, and shock effects within dual-nozzle 
engines. Thermal analyses were performed to define cooling requirements for 
baseline configurations, and special studies of unique dual -nozzle cooling 
problems (i.e. , lip and secondary throat regions) defined feasible means of 
achieving adequate cooling. Similar analysis techniques for the dual- 
expander concept are now available. However, these techniques are, 
in general, less developed than those for the dual-throat concept 
and have not been substantiated with empirical aerodynamic correla- 
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Summary (cont.) 


tions. A summary of these results is contained in the following paragraphs, 
and the technical details are provided in the technical discussion of this 
report. 


The dual -throat aerodynamic model has been expanded to include 
effects not considered prior to the subject contract. Previously, the aero- 
dynamic model was limited in its ability to predict the base pressure that 
developed within the exterior thrust chamber. As presented in Section 
the aerodynamic model computer program has been modified to include a proce- 
dure for predicting base pressure. This model assumes that the base pressure 
is controlled by the loss in total pressure which occurs when the shear layer 
flow passes through an oblique shock during attachment at the exterior 
thruster wall. An empirical constant, called the Nash factor, is used to 
calibrate the model with experimental data. A further development of the 
computer program has been to provide for dissimilar exhaust gases to be used 
for the interior and exterior (bleed) flow. Each stream is allowed to possess 
a given constant value for specific heat at constant pressure, molecular 
weight, and total temperature. 

As a result, the aerodynamic model has been generalized so that it 
can be applied to studies of various dual -throat concepts, including tri pro- 
pellant earth-to-orbit engines and bi propellant engine applications such as 
the Orbit Transfer Vehicle Engine (OTV-E) and liquid oxygen/hydrocarbon 
boosters for uprating the Space Shuttle. This model was also used to define 
test conditions and predictions in support of the MSEC hot-fire testing of a 
dual-throat thruster model. Subsequent analysis of these test results, as 
discussed in Section VI, has shown excellent agreement between model predic-r 
tions and test measurements. 
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Summary (cont.) 


The capability to optimize the dual-throat engine performance was 
enhanced through the development of a computer program that allows the design 
of an optimum nozzle wall contour for a given approach flow. This was 
achieved by modifying an existing Rao nozzle computer program, as described in 
Section As a result, it is now possible to design the contour down- 

stream of the exterior nozzle throat so that it will provide maximum perform- 
ance during Mode II operation. In order to accurately predict the performance 
of nozzles designed using these procedures, revisions were also made to the 
TDK computer program and to the aerodynamic model computer program. This pro- 
cedure was used to define two secondary nozzle contours for the baseline con- 
figuration: one optimized for booster phase operation (Mode I) and one opti- 

mized for sustainer operation (Mode II). In this case, the performance for 
the operating mode in which the nozzle contour was optimized was about one (1) 
second greater than its performance with its non-optimum nozzle. Although a 
minor difference in this case, it shows that the capability now exists to 
optimize the performance for either mode of operation if this proves to be 
advantageous from a total mission standpoint. 

Another analysis task conducted during this study assessed the 
effects of shock waves on nozzle performance. Shock waves are introduced in 
the supersonic expansion by the shape of the jet boundary and by its attach- 
ment to the exterior thruster wall. The resulting shock structure is inher- 
ently two-dimensional and, thus, a two-dimensional numerical procedure must be 
used to predict the flow field. The VNAP computer program was used for this 
purpose. The results of this study task, presented in Section III,D, indicate 
that shocks will have only a minor influence on the performance and thermal 
operating characteristics with properly designed bleed flows. This analysis 
was also used to define boundary conditions for the thermal analysis of the 
adverse condition of Mode II operation without any bleed flow. 

The boundary layer analysis and performance loss of dual-throat 
engines was investigated further by using the JANNAF Reference computer code 
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Summary (cont.) 


BLIMP, in conjunction with the aerodynamic shear model (described in Section 
IIIjC). This procedure provided a rigorous method of assessing boundary layer 
influences, including the effects of bleed flow addition, and guided an 
improvement to the simplified boundary layer loss procedure. 

Finally, all of the aforementioned improvements were incorporated 
into the existing dual-throat simplified performance model, and results were 
compared to previous performance predictions, as noted in Section III,F. In 
general, the revised model predicts lower specific impulse ('*-!%) for the base- 
line configuration in Mode II, but as a result of this investigation, the 
uncertainty of this prediction has been significantly reduced. Parametric 
trends with operating conditions proved similar; however, the variation of 
specific impulse with thrust split is altered and should be considered in 
future mission analysis studies. 

As a second phase of this study, discussed in Section IV, perform- 
ance models for dual -expander thrusters were developed. The dual expander is 
similar in concept to the dual-throat thruster and consists of an exterior, 
annular thrust chamber that surrounds a conventional thrust chamber.. Two modes 
of operation are treated: Mode I, the booster mode of operation, during which 

both engines are on, and Mode II, the sustainer mode of operation, during 
which the annular thruster operates alone. During Mode II operation, a free 
shear layer develops along the jet boundary and is subsequently turned paral- 
lel to the axis of symmetry and transformed into a turbulent wake. Bleed flow 
can be used to control the shape and properties of this viscous region. 

An aerodynamic model similar in concept to the dual-throat model 
was developed to determine the shear layer and base recirculation flow fields 
during Mode II operation. Use of this model showed that the bleed flow 
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requirements for a dual -expander nozzle were significantly influenced by the 
truncation point (Mach number) of the inner nozzle. The overall bleed flow 
required was, in general, 0 to 3%, which is significantly less than that 
required for the dual throat. In fact, it was concluded that the dual- 
expander engine could be designed and operated without bleed flow while still 
obtaining controllable aerodynamic flow characteristics. 

A model for predicting shock losses and flow effects during Mode 
II operation was also applied to predict shock location and performance 
effects with varying inner nozzle truncation and bleed flow. The results 
indicated little performance differences with nozzle truncation and no signi- 
ficant shock influences over the range of bleed flows considered. However, 
the VNAP results were not sufficiently accurate in its present form to be used 
for performance prediction or for optimum nozzle design. To overcome this 
problem, development of a method of characteristics analysis with shock capa- 
bility is recommended. Most of the point procedures required for this 
approach are available and could be used in assembling the new computer pro- 
gram. 


Procedures were also developed to specify preliminary geometric 
design parameters for a dual-expander nozzle. These procedures are based on 
the criterion that the flow angle and static pressure must be matched at the 
primary nozzle lip during Mode I operation. During this study, the VNAP 
analysis was used to calculate the flow characteristics in the annular nozzle 
region. The subsequent results verified that the design criteria had been 
met. 


Thermal design and parametric analyses, described in Section V, 
were performed for both dual-nozzle concepts. A satisfactory dual-throat 
baseline design was obtained by using three parallel cooling circuits. 
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The primary circuit coolant flows from an area ratio of 1.58:1 in the primary 
nozzle to the primary injector. This circuit has the highest pressure drop 
requirement of the three parallel circuits, so its coolant flow fraction is 
selected to provide the minimum pressure drop. A separate circuit cools the 
inner annulus surface, the lip region, and the rest of the primary nozzle. 

The majority of the channel pressure drop in this circuit is required for the 
lip region, for which a three-dimensional SINDA model was developed. A 
counterflow circuit cools the outer surface of the secondary chamber; this 
circuit is a composite design which considers both Mode I and Mode II oper- 
ation. The channel design over a short region upstream of the secondary 
throat is defined by the primary plume impingement associated with a bleed 
flow malfunction in Mode II. An exact solution for the inviscid plume 
impingement problem was obtained with the VNAP program. This solution 
provided edge conditions for an integral boundary layer analysis which defined 
the wall heat transfer edge conditions. The initial momentum thickness for 
the boundary layer analysis was varied in order to bound the problem, with the 
channel design based on the worst case. The maximum heat flux ranged from 
essentially the Mode I value to almost 50% higher. 

Individual cooling circuit characteristics for the dual-throat baseline 
design have been defined. Although the fuel system pressure drop of 1160 psi 
is set by the primary circuit, the coolant Mach number is slightly higher in 
the secondary circuit when the latter is designed to accommodate the maximum 
postulated Mode II heat flux for the zero bleed flow case. Parametric study 
results indicate that it is the coolant Mach number, rather than pressure 
drop, which limits increases in chamber pressure and Mode I thrust split. 
Increasing the chamber pressures to 4500/3150 psia or the thrust split to 
80/20 results in a primary circuit Mach number of 0.37 (an unacceptably high 
design value). 
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The dual -expander designs, similar to those of the dual throat, are 
cooled with three parallel circuits. However, the inner annulus circuit is 
now a major circuit in terms of length and heat load and also includes a 
throat region. Regenerative cooling of the baselne design point with hydrogen 
only is difficult. Coolant flows in the secondary and inner annulus circuits 
must be selected to avoid minimum channel depth and Mach number limits, 
respectively, and the remaining flow available for the primary is less than 
optimum. As a result, the system pressure drop of 1720 psi is 130 psi above 
the minimum primary circuit value. Using oxygen to cool the secondary circuit 
and reducing the interface area ratio between the primary circuit inlet and 
the inner annulus circuit outlet allows the hydrogen pressure drop to be 
reduced to 1550 psi. Transpiration cooling the throat region of the primary 
chamber did not prove to be effective since the hydrogen flow fractions 
required for transpiration cooling, coupled with the flow requirements of the 
other circuits, severely limit the flow available for regeneratively cooling 
the remainder of the primary circuit. This results in relatively high pressure 
drops. 


Parametric study designs of the dual expander could not be obtained 
with increased chamber pressures or increased Mode I thrust split, or with the 
secondary chamber pressure reduced to 1500 psi a. The latter results from the 
high primary chamber pressure drop associated with reduced coolant inlet pres- 
sure. These results, along with the baseline design details, indicate that 
the baseline design point is near the optimum chamber pressure ratio and is 
very close to the chamber pressure versus thrust split characteristic which 
represents a practical cooling limit. 
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I. 


INTRODUCTION 


Propulsion systems for future vehicles, such as the Single-Stage-to- 
Orbit (SSTO) and Heavy Lift Launch Vehicle (HLLV), may embrace such capabili- 
ties as dual-mode operation and in-flight changes in area ratio for altitude 
compensation. These vehicles benefit from dual-mode operation through reduced 
vehicle volume by taking advantage of relatively high bulk density propellants 
in one mode and lower density, but higher performing propellants in the other 
mode. Area ratio change during flight provides an increase in performance as 
ambient pressure decreases with altitude. 

Dual-nozzle engines combine both operating capabilities in a single 
design. Their dual combustors allow use of either one or two propellant com- 
binations and, with their two separate nozzle throats and a fixed nozzle exit 
area, allow for a shift in area ratio without resorting to nozzle translating 
mechanisms. 

The dual -nozzle concept can be applied to both bi propellant and tri pro- 
pellant engines with resultant advantages. The baseline engine application, 
analzyed in this and earlier studies, is a tri propel 1 ant single-stage-to-orbit 
(SSTO) engine. Dual nozzles could also be applied to bi propellant dual -thrust 
Orbit Transfer Vehicle (OTV) engine and Liquid Rocket Booster (LRB) engine 
designs. For OTV-E, a high-thrust priority payload and a low-thrust cargo 
transfer engine could be combined in a single thrust chamber. For a LOX/ 
Hydrocarbon LRB, mission thrust tailoring and an area ratio increase at alti- 
tude could be achieved. This concept has been analyzed, with favorable 
results, during the recent Advanced Oxygen-Hydrocarbon Rocket Engine Study 
(Ref. 1). 

Two types of dual -nozzle designs have been conceived — the dual throat 
and the dual expander. An engine system preliminary analysis using the dual- 
throat concept was performed during Contract NAS 8-32967 to examine potential 
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I, Introduction (cont.) 


power cycles and generate parametric data for a tripropellant SSTO vehicle 
engine. A preliminary performance prediction methodology based on a new aero- 
dynamic bleed flow computer model was developed. This model was formulated by 
using the results of cold-flow tests conducted with a subscale dual-throat 
thruster configuration. The dual-expander engine concept, prior to the sub- 
ject contract, has received less formal analysis. Some preliminary engine 
system parametric data have been generated, and a conceptual baseline engine 
system has been established. 

A dual-throat thruster consists of a large thrust chamber assembly that 
contains within it a smaller thrust chamber. Because of this novel construc- 
tion, dual-throat thrusters provide the means to obtain a large area ratio 
adjustment within a single thrust chamber assembly without the need for 
extendible nozzles. This area ratio adjustment is accomplished by changing 
the operating mode of the thruster. For example, in the high-thrust mode, the 
exterior and interior thrust chambers are both operating. When the exterior 
thrust chamber is extinguished, the low-thrust interior thrust chamber remains 
in operation, but now utilizes the exterior thrust chamber exhaust nozzle as 
its exhaust expansion contour. In this manner, the interior thruster is 
caused to operate at a high expansion ratio. The attachment of the interior 
thruster exhaust to the exterior thruster wall contour can be controlled by 
injecting mass flow into the exterior chamber, thereby creating a free shear 
layer that develops along the jet boundary. Mass injection will raise the 
base pressure in the chamber of the exterior thruster and move the jet 
attachment point downstream. Proper location of the free jet can serve to 
lower the losses in interior chamber performance that are caused by 
aerodynamic effects. 

The sustainer mode of operation, as described above, is called Mode II, 
and the booster mode of operation is called Mode I. The dual-throat thrust 
chamber configuration and flow fields for Mode I and II of the tri propellant 
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I, Introduction (cont.) 


SSTO engine case are illustrated in Figure I-l. During Contract NAS 8-32967, 
the baseline engine design gas generator flow was employed as the Mode II 
bleed flow. During Mode I, the gas generator flow would be dumped in the 
nozzle exhaust. 

The dual expander is similar in concept to the dual-throat thruster. 

The dual expander consists of an exterior annuluar thrust chamber that sur- 
rounds a conventional thrust chamber. Figure 1-2 illustrates the two modes of 
operation for the SSTO engine application that are treated in this study. 

Mode I is the booster mode of operation during which both engines are on. The 
annular thruster uses a LOX/H 2 propellant combination, and the interior 
thruster uses LOX/RP-1. Thus, in the dual expander, the exterior chamber is 
the sustainer thruster whereas, in the dual-throat design, the interior 
chamber is the sustainer thruster. Mode II is the sustainer mode of operation 
during which the annular thruster operates alone. During Mode II operation, a 
free shear layer develops along the jet boundary and is subsequently turned 
parallel to the axis of symmetry and transformed into a turbulent wake. Bleed 
flow can be used to control the shape and properties of this viscous region. 

Technological investigation of the dual-throat nozzle concept was ini- 
tiated in 1975 by the NASA George C. Marshall Space Flight Center. In July of 
1976, the NASA and the Aerojet Liquid Rocket Company (ALRC) collaborated on an 
analytical/empirical program to evaluate the concept. A bailment agreement 
was established wherein NASA hardware was supplied for use on an ALRC- 
sponsored cold-flow program. Aerodynamic analyses were conducted, and the 
cold-flow test program \^as completed. After the promising aerodynamic fea- 
tures had been confirmed, NASA contracted ALRC (NAS 8-32666, Ref. 1) to expand 
the data base and to initiate development of a dual-throat performance predic- 
tion model . 
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a. MODE I, BOOSTER MODE OPERATION 
DUAL- THROAT THRUSTER 



b. MODE II, SUSTAINER MODE OPERATION 
DUAL-THROAT THRUSTER 


Figure I-l. Dual-Throat Thrust Chamber Configuration 
and Flow Fields for Modes I and II 
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b. DUAL-EXPANDER MODE II OPERATION 
OUTER ENGINE ON, INNER ENGINE OFF 


Figure 1-2. Dual-Expander Modes of Operation 





I, Introduction (cont.) 


The Dual-Fuel, Dual-Throat Engine Preliminary Analysis Program (NAS 
8-32967, Ref. 2) was subsequently conducted to determine dual-throat engine 
system parameters and vehicle/engine integration requirements for the tri- 
propellant SSTO engine application. This study concluded that while the 
dual-throat, dual-fuel engine was a viable SSTO candidate, a number of uncer- 
tainties remaining in the heat transfer and performance areas needed to be 
resolved. Also, as the Dual-Expander Engine had been conceived during this 
period (Ref. 3), there was a need to make accurate trades between the two 
dual-nozzle concepts. As a result, the effort described in this report was 
conducted. 

To accomplish this, the Dual Nozzle Aerodynamic and Cooling Analysis 
study was conducted. The results are discussed in this report. The major 
objectives of the program were as follows: 

(1) Improve the existing dual -throat aerodynamic bleed flow model and 
performance prediction methodology; 

(2) Develop basic geometric and aerodynamic data defining performance 
for a dual -expander chamber; 

(3) Conduct preliminary thermal analyses of both the dual throat and 
dual expander; 

(4) Provide engineering analysis of data from government hot-fire 
testing of a dual-throat thruster model chamber. 
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I, Introduction (cont.) 


The objectives were accomplished through a series of five tasks: 

Task I: Dual-throat aerodynamic model and performance prediction 

improvements. This task resulted in improving the bleed 
flow model to analytically predict base pressure for 
dissimilar compositions of the primary and bleed flows. 
Secondly, a computer program was developed to optimally 
design the secondary nozzle contour for Mode II operation. 
Third, simplified procedures used to predict boundary layer 
and divergence losses were verified by a rigorous analysis. 
This analysis also investigated the influence of shocks on 
engine operation. 

Task II: Dual-expander preliminary geometric analysis. A design cri- 

terion was established, and the interrelationships between 
nozzle geometry and engine operating characteristics were 
defined. 

Task III: Dual-expander preliminary flow field analysis. The flow 

field characteristics of the dual-expander operation were 
investigated, verifying the design approach from Task II and 
identifying the insensitivity of nozzle performance to inner 
chamber truncation and the low bleed flow requirements of 
the dual expander. 

Task IV: Preliminary thermal analysis of both dual throat and dual 

expander. This task showed that the lip regions of the 
inner nozzle in both concepts could be cooled. Completely 
regeneratively cooled designs were established for both 
engines. 

Task V: Test data analysis. Hot-fire data analysis verified the 

improved aerodynamic bleed flow model. 
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II. 


CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The major conclusions formed during conduct of the program are 
listed below. 

1. The dual-throat aerodynamic bleed flow model has been veri- 
fied on the basis of cold-flow and hot-fire data correla- 
tions and is ready for use as a general design tool. 

2. The predictive accuracy of the aerodynamic bleed flow model 
could be bettered by improving the current plume shape calcu 
lational procedure. 

3. Dual-throat engine performance is relatively insensitive as 
to whether or not the secondary nozzle contour is optimized 
for the Mode I or Mode II expansion. 

4. The previous simplified technique for predicting the dual- 
throat Mode II boundary loss was inaccurate. The rigorous 
procedure employed during this program has resulted in 
significant improvement of the simplified prediction proce- 
dure. 

5. Dual -throat nozzle performance in Mode I is comparable to 
that of a conventional nozzle. Mode II performance is 
approximately 1-2% lower than that of a conventional nozzle 
optimized for altitude (i.e., sustainer) operation. 
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II, A, Conclusions (cont.) 


6. Preliminary performance analysis results for the dual- 
expander nozzle indicate that high Mode II expansion effi- 
ciencies can be achieved. The results indicate low bleed 
flow requirements and insensitivity to the primary nozzle 
truncation area ratio- 

7. There is a degree of uncertainty associated with the dual- 
expander performance results. The calculational procedures 
employed during this program (i.e., the VNAP computer model) 
should be replaced by a method of characteristics solution 
that treats shock waves. 

8. Both the dual-throat and dual-expander baseline tripropellant 
SSTO engine designs can be cooled with three-circuit regener- 
ative designs. The primary nozzle lip and secondary nozzle 
throat gas-side heat transfer coefficient assumptions should 
be verified through hot-fire testing. 

B. RECOMMENDATIONS 

The following recommendations for future work are made on the 
basis of the significant results obtained on this program. They are listed in 
the order of highest to lowest priority. 

1. A demonstration program should be conducted to verify the 
major thermal design assumptions utilized during the study. 
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II, B, Recommendations (cont.) 


2. The dual-expander flov^ field analysis results should be veri- 
fied through conductance of a model cold-flow test program. 

3. The dual -throat SSTO engine parametrics developed on Contract 
NAS 8-32967 should be updated on the basis of the improved 
performance prediction results obtained on this program. 

4. A method of characteristics solution that treats flow field 
shocks should be developed for Mode II dual expander 
analysis. 

5. The plume shape calculational procedure contained within the 
aerodynamic bleed flow model should be improved and corre- 
lated with existing test data. 

6. A procedure for optimization of dual -expander secondary 
nozzle contour should be developed. 
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III. DUAL THROAT AERODYNAMIC MODEL AND PERFORMANCE PREDICTION 
IMPROVEMENT 


The basic objectives of Task I were to improve the dual throat aero- 
dynamic model's capability and verify the approaches used in the dual throat 
performance prediction methodology developed by ALRC for NASA/MSFC during 
Contract NAS 8-32666. 

The specific objectives of Task I were as follows: 

° Improve the base flow pressure and shear layer portion of the 
aerodynamic bleed flow model. 

° Assess the approach used to calculate boundary layer loss for Mode 
II operation. 

“ Assess the approach used to calculate nozzle divergence efficiency 
for Mode II operation and determine how shocks affect it. 

° Develop a Rao nozzle optimization program for the Mode II opera- 
tion and determine the effect of Mode II contour optimization on 
the delivered performance of both modes. 

° Assess changes in predicted performance based on the use of the 
improved aerodynamic model and update the performance model as 
required. 

° Review the existing cold-flow data bank and define areas that 
would benefit from further cold-flow test data. 
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Ill, Dual Throat Aerodynamic Model and Performance Prediction 
Improvement (cont.) 

A. LITERATURE SURVEY 

In support of the improvement of the aerodynamic model , a compu- 
terized literature survey was conducted for the purpose of obtaining refer- 
ences that might be useful in formulating the analytical model as well as data 
that could be used for calibration of the model. Conducting a computerized 
literature search has the advantage of significantly reducing the amount of 
man-hours otherwise required for this task, thus significantly reducing cost 
and time as well. 

The literature search was carried out via a remote data terminal 
to gain access to a computerized library. The method used involves the formu- 
lation of a search strategy based on key words or phrases. These may be com- 
bined using bool i an operators to obtain all literature from a given data base. 
The search is conducted interactively on the main frame computer. For 
example, reference titles that contain a given key word, but not another key 
word, may be listed on the terminal. If a particular title appears to be of 
interest, a request by terminal is made and the abstract of the article is 
printed. If desired, the article can then be ordered over the terminal. 

In order to prepare a search strategy, the initial computer ses- 
sion was spent in checking various key author names and key words for fre- 
quency of citation in the NTIS index. Titles obtained in this way were used 
to construct the bool i an operators to be used in the general search. A 
strategy was established to be accomplished in two steps. 

The first step was to identify those papers which make reference 
to the principal works in the field. This was accomplished with the Science 
Citation Index capability of SCISEARCH. Papers which referred to the works of 
Korst, Bauer, and Alber were sought in this way. 
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Ill, A, Literature Survey (cont.) 


The second part of the literature search utilized various key 
words and employed six data bases. Information obtained in this manner is 
summarized in Table III-I, which gives the number of references found by 
category by source. Abstracts were obtained for these references. After 
reading the abstracts, approximately 30 papers were ordered. 

B. AERODYNAMIC BLEED FLOW MODEL 

The aerodynamic bleed flow model's capability of predicting the 
aerodynamics unique to Mode II operation has been expanded. The Mode II (sus- 
tainer mode) operation is illustrated in Figure III-l. The interior thrust 
chamber is called the primary thruster, and its exhaust is called the primary 
flow stream. A secondary flow stream (bleed flow) is injected into the 
exterior thrust chamber. The thruster expands the primary flow to a super- 
sonic condition at the nozzle exit. Further expansion then occurs in the form 
of a Prandtl-Meyer fan at the nozzle lip. The result is an exhaust plume with 
a constant pressure boundary. Flow is introduced into the exterior chamber so 
as to control the location of the plume boundary in order to minimize 
reattachment shocks which occur when the flow impinges on the exterior 
thruster wall. The plume boundary is an exhaust streamline path and acts very 
much like a nozzle wall. However, a shear layer develops along this boundary 
due to viscous interaction of the exhaust jet and the gases recirculating in 
the exterior chamber. Analysis of this shear layer mixing region and of the 
phenomena described above are an important feature of the aerodynamic model. 

The principal assumptions employed in the aerodynamic model are 
listed below: 


The flow leaving the primary nozzle is one-dimensional and 
supersonic. 
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TABLE III-I 


SUMMARY OF LITERATURE SEARCH-AEROJET DUAL NOZZLE ANALYSIS 


DATA BASS 

~T 



SHEAR LAYERS 

16 


f 

22 



SHEAR LAYERS * HEAT TRANSFER 

4 






SHEAR LAYERS * (HT. TRANS. + ISOENERGETIC + REATTACHMENT) 


22 

8 


2 


BASE FLOW + BASE PRESSURE 

42 


22 

13 

20 

20 

(BASE FLOW + BASE PRESSURE) * (SHEAR LAYERS + HT. TRANS. + REATT) 
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Figure III-1. Mode II, Sustainer Mode Operation 



Ill, B, Aerodynamic Bleed Flow Model (cont.) 


The primary nozzle exhaust flow is represented by an ideal 
gas with constant specific heat and molecular weight. 

The bleed flow is represented by an ideal gas with constant 
specific heat and molecular weight. These values can differ 
from those of the primary exhaust stream. The total temper- 
ature of the bleed flow can also differ from that of the pri- 
mary flow stream. 

The shear layer is represented by a model of the type 
developed by Chapman, Korst, Bauer, Chow and Addey, and 
others (see Refs. 5, 6, 7, 8, and 9). The model used here, 
however, is non-isoenergetic, as defined by the assumptions 
given above. 

The shear layer is treated as two-dimensional. This assump- 
tion has been shown to be accurate (see Bauer, Ref. 7), pro- 
vided that the projected thickness of the mixing zone on the 
radius of symmetry is less than .3 of the radius of symmetry. 
This condition is always satisfied in dual throat thrusters. 

The boundary layer at the primary nozzle lip is thin. 

The shear layer is a constant pressure surface, both in cross 
section and along its entire length. At the inner edge of 
the shear layer, the gas velocity is negligible. At the 
outer edge of the shear layer, the velocity is constant and 
equal to the velocity of the plume boundary streamline. 

The shear layer is superimposed on the inviscid jet plume 
boundary. 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


° The plume boundary can be represented by the methods devel- 
oped by Herron (Ref. 10). 

1. Shear Layer Modeling 

Some features of the free shear layer that develops along the 
jet boundary are shown in Figure III-2. The shear layer is assumed to begin 
at the exit of the primary nozzle where its thickness is negligible. It then 
develops as a turbulent shear layer along the boundary of the primary nozzle 
exhaust plume. Viscous interaction of the exhaust plume with gases recircu- 
lating in the chamber of the secondary nozzle form the shear layer. In 
steady-state operation, the system will maintain a stable equilibrium base 
pressure. This base pressure determines the plume boundary. The theory 
assumes that a streamline divides the shear layer such that gases on one side 
of the streamline are recirculated into the base region. This is the 
so-called "dividing streamline" concept. At steady-state operation, all of 
the gases exhausting from the primary nozzle will also exhaust from the secon- 
dary nozzle. Gases in the recirculating region (in this case the chamber of 
the secondary nozzle) have been trapped because they are degraded in total 
pressure to the extent that they cannot penetrate the static pressure rise 
produced by turning off the flow at the point where the shear layer impinges 
on the secondary nozzle wall. It is assumed, in this model, that the dividing 
streamline becomes a stagnation point at the secondary wall. This stagnation 
occurs abruptly either by a compression that is essentially isentropic, or by 
a non-isentropic compression through a shock structure. Details of this pro- 
cess govern the value that will be obtained for a base pressure, i.e., the 
pressure along the plume boundary. 

The introduction of mass flow into the recirculation region 
will tend to raise the base pressure and shift the plume boundary downstream. 
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Figure III-2. Shear Layer Streamlines and Velocity Profile 



Figure III-3. The "d" and "j" Streamlines within the Shear Layer 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


Mass flow introduced into the recirculation region at a given rate must exit 
the nozzle at the same rate. If this flowrate is sufficiently small, all of 
the injected flow will become entrained within the shear layer. For this 
reason, the analysis used here is called a "small bleed flow" model. 

Figure III-3 shows the location of two important streamlines 
within the shear layer. The "d" streamline is the so-called "dividing stream- 
line" that separates the recirculation gases from the gases that exit the 
nozzle. The recirculation gases are assumed to be a mixture or primary and 
bleed flow gases. The "j" streamline, also shown in Figure III-3, is a 
limiting streamline for the bleed flow gas. On the primary side of the "j" 
streamline, there is "pure" primary flow, whereas, on the secondary side, 
there is "mixed" primary and secondary flow. All of the bleed flow introduced 
during steady-state operation must exit the nozzle and is contained between 
the "d" and "j" streamlines. The "d" and "j" streamlines coincide if the 
bleed flow is reduced to zero. 

The velocity profile at any position across the shear layer 
will appear as shown in Figure III-2. The shear layer will be turbulent, but 
for either laminar or turbulent flow the velocity profile can be represented 
by 

0 = (1 + erf n) / 2 


where 

0 = u/u , 0< 0 <1 
00 

and erf is the error function, i.e., 

n 

S 2 

exp (-n ) dn 

0 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


n is the nondimensional mixing length 

n = ay/l 

In the above equation, y is the distance measured normal to the centerline of 
the shear layer, and l is the length of the shear layer measured along the 
centerline from its origin. The variable a is known as the "shape factor" or 
"similarity parameter." The difference between laminar and turbulent mixing 
is reflected in the relationship used for a. It is necessary to use an empir- 
ical expression for the shape factor. Numerous expressions fora are avail- 
able in the literature (see Refs. 6, 7, 9, 11), but the empirical expression 
given by Korst (Ref. 6) has been used here, i.e., 

a = 12 + 2.758 Mj^ 

where Mb is the Mach number along the inviscid primary jet boundary. 

Since the bleed flow gas properties can be different from 
those in the primary stream, it is necessary to assume a variation in pro- 
perties across the mixing zone, i.e., across the shear layer from the base 
recirculation region to the "j" streamline. It is assumed that the distribu- 
tion of mass fraction across the mixing zone is linear with respect to velo- 
city ratio, thus: 


C = C (1-d) + 0C 
p pr ' pp 

where Cp^ is the specific heat in the recirculation region, i.e.. 


^pr ^s ^^“'^pr^ Sp 

Yp^ is the amount of mass from the primary flow entrained in the recirculation 
region. 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 
The gas constant in the mixing region is 


R = R^ (1-0) + 0Rp 

where R^ is the gas constant in the recirculation region, i.e.. 


with 


R, (1-Y ) + R 

s ' pr' p 



Rs 


= R/Mws 


and 

Rp = R/Mwp 

where R is the universal gas constant. The specific heat ratio in the mixing 
zone is 


Using the relations 


it follows that 


Y 


Y = 


v<s 


- R), 


R = C - C 
Y = Cp/C 


Y^ (1-iZi) + Yp 



vr 


0 


( 1 - 0 ) + 


Sr 



0 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


In like manner, the specific heat ratio along the dividing streamline is found 
to be 


Vr C-«'d> * ’'p "d 


pr 


Next, the velocity ratio at the dividing streamline, is sought. The 
energy equation is 


so that 


Cp T + u 2/2 = Cp 
= 2Cp (T^ - T) 


and 




= C . (T 
pd ' 


od ~ ^d^^*"pp ^^op 


Tb) 


where T5 is the temperature along the plume boundary. Defining the Crocco 
number as 




gives 


= (t^d^ Cp, T„,)/(Cp2 Cpp T„p) 

Expressing the Crocco number in terms of pressure gives 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


and 



Cb/Pod) 




where Pod is the total pressure on the "d" streamline. However, at the 
shear layer attachment point, the "d" streamline is stagnated after passing 
through an oblique shock. The maximum strength of this shock corresponds to 
the deflection angle experienced by the plume boundary when turned abruptly by 
the nozzle wall , i .e. , 


Pod^Pb = P'Pb 


where P/Pb is obtained from the oblique shock relations. On the other hand, 
if the shock imparts no loss in total pressure to the "d" streamline, then the 
flow is compressed i sentropical ly and 


Pod = Pb 

It is assumed that the true recompression criterion falls between these two 
extremes, so that 


and 


Pod = + (1-N) Pb 


0 < N < 1. 


The quantity, N, is known as the Nash factor (see Ref. 11, 12). In the 
analysis presented here, N is used to calibrate the shear layer model with 
experimental data. 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


Total enthalpy across the mixing layer is also assumed to 
vary linearly with respect to velocity ratio, i.e., at the "d" streamline 



so that 


Sd ^od ^ 

Sp ^op 


^pp ^op 



C T \ 

pr r \ 

C J I 

PP op / 


0 


d 


From the definition of the Crocco number, and since = Tod» it follows 
that 



Solving the above equations gives: 




"PP 




+ 4 Sr 

Sp Sp 


The expressions given above are sufficient to determine Y^j, Cd^, and 

by means of numerical iteration, provided that conditions in the base 
recirculation region are known. The similarity coordinate of the "d" stream' 
line is then given by 

Tid = erf"^ (2 -1 ) 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


The similarity coordinate for the "j" streamline, n. can be found from the 
mass and momentum conservation equations. Using these equations, it can be 
shown that for a two-dimensional shear layer* 

^ ip/pj) = 



(p/p_) 0 (1-0) dn 


The term (p/p^), which is the ratio of the gas density across the shear layer 
to the density at the primary edge, P[j» can be evaluated as follows: 


(p/p„) = (p/P|j) 


From the perfect gas law, and since the pressure is constant across the shear 
layer. 


(p/Pfj) 


R 


R 

_2. 

R 


Tb/Top 

T/T„ 


ZoR 

To 


Since 


T/Tq = 1 - c2 


then 


(p/p[j) 


R 


1-C 


1-C 



*For example, see Bauer (Ref. 7) and neglect the axi symmetric terms. 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


As shown previously, it follows directly from the energy equation that 


C C T 


r r T 
V op 


Therefore, 


(p/Pk) = 


R 

R 


1-C 

? 2 T T 

1-C^ 0 pp op 

P 0 


Zo£ 

To 


(p/Pb) - 


1-C. 


To „ 2 .2 ^pp 

op p 


and 


R_ C_ 1 - C, 


(p/pk) = / ^ 


2“X 


pp A- C. 0 


where A is the enthalpy ratio, 

A - \ ’^0 


T C 
op pp 


and is assumed to vary linearly with velocity across the shear layer as follows 


Sr Tr + n - Tp,. Ty, 
^PP Sp \ ^pp Tqp 


) 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 
where 

C = C (1-0) + 0 C 
p pr ' ' PP 


R = (1-0) + 0 Rp 


The similarity coordinate, ns is then found from the following relationship: 

0 3 


J 

-3 


C 

^PP 


R 

_£ 

R 


A - 


r 2 ,2 

Cjj 0 


dn 


^ Rp 0 (1-0) 

^nn R 
PP 


A - Cjj 3 02 


dn. 


The far edges of the shear layer are taken as +3 and -3. The value for 
is obtained by integrating the left-hand side of the above expression until it 
is equal in value to the right-hand side. For computational convenience, the 
followng functions are defined; 


A 





c 



-3 

PP 



A 

I3 (A) = ^ 

c 



\ 

C 



/ 

PP 

so that 


-3 




3 

d dn = C 

R 

Cnn R 
PP 

-3 Sp 

R 

- Cb ^ J 


R 

R 


R 


0 

A - 0^ 


2 — ? 

A - C " 0^ 


dn 


dn 


2^ 


A - c. ^ 


dn- ( "d ‘p 

C„'- R . . 2 ^2 


dn 


-3 


PP 


A - C.^ 


Il (nj) = (3) - I 3 (3) 

Next, a method is presented for determining the species concentrations in the 
base recirculation region. The shear layer to the left of the "d" streamline 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


is turned into the recirculation region at a height r^^ above the primary 
flow centerline. By applying the mass and energy balance equations to a 
control volume surrounding the recirculation region, it can be shown that 
(Ref. 13) equation [18]) 

"s 'VBleed - Vr) = ^ p U(Vp^ - Y^) dy 


The bleed flow is entirely made up of the secondary flow species, so that 


^PBleed " ^p/^ bleed ° 


Thus, 


0 ■ v„ + ^ ^d pu (V„„ - Y„) dy 


pr 




pr p' 


For unity, Schmidt and Lewis No. 


Y - Y = j5(l-Y) 
p pr ^ ^ pr' 


So that 


2irr ^d 

0 = Y r + 

W. 


$ 


pU0dy (1-Yp^) 


Changing to nondimensional coordinates, and proceeding as before, gives 


0= ( V - b - ] O-Cb^) (n,) (l-Yp^)] 


Ws \ a 


-) 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


It can be shown that the mass flow in the secondary stream is* 




W. = ^ [I] (nj) - Ii (n^j)] 


Combining the above two equations gives 

^3 ^'"d^ 


where 


tt = 


II (Oj) - (Oj) 


a = Y /Y 
pr sr 


0 ) 


therefore 


The total enthalpy in the base recirculation region is 


Ss ^os 

^PP ^op 


Sri 

Sp ^op 


Si 

^PP ^Op 


( 2 ) 


The total mass flow from the primary nozzle is 


"op V 


^p*^p^op 


(■^r) 


3!L 
2^) 


Thus, the ratio of the secondary mass flowrate to the primary mass flowr‘ate 


Yp+1 


W./W^ 
s p 


2 
a 


(r 




2YvT) 


[I, (lj) - I, (rij)] (3) 


*For example, see Bauer (Ref. 7) and neglect the axisymmetric terms. 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


Equations (1), (2), and (3) above form a system of three nonlinear simultan- 
eous equations in three unknowns: Yp^, Pb/Pop> ^nd T^/Top* The 

aerodynamic bleed flow computer program numerically solves for these three 
unknowns, using input initial estimates of their values. 

2. Aerodynamic Bleed Flow Model Computer Program 

The small bleed flow model for the shear layer derived at in 
the previous subsection is used in the aerodynamic computer bleed flow pro- 
gram. A flow chart for this computer program is presented in Figure III-4. 

The procedures shown in Figure 1 1 1-4 are discussed below. 

At the beginning of the program, initial and default values 

are set and the input data is read, using FORTRAN NAMELIST. (See Appendix A 
for a description of the program input and a sample problem.) The primary and 
secondary wall geometries are then constructed. Both the primary and the 
secondary nozzles are assumed to have geometries that can be described ana- 
lytically by the following parameters: 

Chamber contraction ratio 

Circular wall radius connecting the chamber and the nozzle 

i nl et 

Conical nozzle inlet 

Circular wall radius upstream of the nozzle throat 

Circular wall radius downstream of the nozzle throat 

The nozzle wall geometries are defined by simple geometrical inputs such as 
wall slope, throat radius of curvature, etc. Complete shapes are specified so 
that the wall contours can be plotted as well as used in the analysis. 
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Figure IH-4. Flow Chart for the Aerodynamic Bleed Flow Model Computer Program 












Ill, B, Aerodynamic Bleed Flow Model (cont.) 


In order to be used easily by the computer programs, the wall geometries are 
converted to finely spaced tables of points. 


Next, an iterative procedure is initiated which consists of 
the steps enclosed in the dashed line rectangle shown in Figure III-4. The 
steps taken in this procedure are individually discussed below. 


(1) Conditions at the exit of the primary nozzle are deter- 
mined in the MAIN subroutine by using the one-dimensional 
relation between Mach number and area ratio, i.e., 


e = M 
P e 


-1 








9 


2 - 1 ^ 1 ) 


(2) The exhaust flow is turned at the nozzle lip and expands 
to equal the boundary Mach number M^. The angle of this 
turn, Au, is determined in the PLUME subroutine by using 
the Prandtl-Meyer relations for supersonic expansion 
at a point, as follows: 
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IIIjB, Aerodynamic Bleed Flow Model (cont.) 

(3) The plume boundary is determined by using the modified 
method of Herron (Ref. 10). Since the plume boundary is 
a circular arc, all of its geometric properties (such as 
its position, arc length to any point, and slope at any 
point) are known analytically. This step is 
accomplished by using subroutine PLUME. 

(4) The intersection of the plume with the secondary nozzle 
wall is determined geometrically by using subroutine 
INTERS. 

(5) The shear layer model allows the shear layer profile to 
be determined exactly at the point of intersection of 
the plume boundary and the wall. All properties of the 
profile can be evaluated from various integrals of the 
error function, as described in Section III,B,1. 

These integrals are evaluated in subroutines lONEV, 

ITHRV, IlINVR, and IZEROS. 

(6) The plume boundary is a streamline. The streamline will 
be turned abruptly (i.e., deflected) to follow the 
secondary nozzle wall, resulting in a shock structure. 

The cases of most interest correspond to solutions of 
the shock relations of the weak family. Cases where a 
lambda shock structure develops, so that the flow is 
shocked subsonically in part of the region, are of less 
interest since they correspond to aerodynamical ly poor 
designs. The obliquie shock relations are used to deter- 
mine the rise in static pressure across the shock and 
the loss of total pressure across the shock. Since the 
development of shock structure in the nozzle is funda- 
mentally a two-dimensional phenomenon, no attempt is 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


made to treat the flow in the nozzle downstream of the 
origin of the shock. The oblique shock relations are 
calculated by subroutine WAVANG. 

The calculations described above require a first estimate for 
the variables P^, T^, and Ypr* It is necessary to find converged values 
for these three variables that satisfy equations (1), (2) and (3) of Section 
III,B,1. Newton's method is used to obtain the next estimates for P|j, and 
Tp. Equations (2) and (3) are used to obtain the four partial derivatives 
required by Newton's method. The next estimate for Yp is obtained by the 
method of iteration, using equation (1). Solutions are sought for each bleed 
flow given in an input table. Converged values for a solution are automatic- 
ally used as first estimates for the next case. 

The primary nozzle contour, the plume shape, and the secon- 
dary nozzle contour can be plotted. Punched card output is available for the 
outer streamline for the inviscid flow, i.e., the primary wall, the plume 
boundary, and the secondary wall downstream of the plume attachment point. 

The punched cards are in the input format required by the TDK program and/or 
the VNAP program (Refs. 3 and 4, respectively). 

3. Model Calibration 


Three empirical parameters exist in the bleed flow model 
which can be adjusted, as necessary, to improve the correlation between the 
Mode II bleed flow and the recirculation pressure and temperature. Two of 
these parameters were adjusted by using the cold-flow data and hardware 
geometry from NAS 8-32666. 

One parameter, the plume scaling factor (FCTR), alters the 
plume shape. This is accomplished by using the plume scaling factor as a 



Ill, B, Aerodynamic Bleed Flow Model (cont.) 


multiplier for the similarity parameter, in the Herron plume correla- 

tion. (Mb is the Mach number along the plume boundary, and y is the ratio 
of specific heats.) Altering the plume scaling influences the plume impinge- 
ment angle and, hence, the resulting shock strength. 

A plume scaling factor of 1.2 resulted in model predictions 
that agreed well with the experimental data for small nozzle spacing as shown 
in Figure III-5. However, for tests with larger nozzle spacing, the model 
predictions became less accurate when a constant scaling factor, as shown in 
Figure III-6, was used. Decreasing the plume scaling factor to 0.9 improved 
the accuracy for large nozzle spacing (Figure III-7). A correlation between 
plume scaling factor and nozzle spacing, as shown in Figure III-8, can be 
inferred, but the validity of applying it to other dual throat geometries 
needs to be determined. 

The second parameter, the Nash factor (N), is used in the 
recompression criteria relating the pressure behind the plume attachment shock 
to the recirculation pressure. This parameter adjusted the stagnation pres- 
sure loss to match the cold-flow data. While the plume scaling factor shifts 
the base pressure versus bleed flow curve up or down, the Nash factor changes 
the actual shape of the curve (Figure III-9). A Nash factor of 0.4 was found 
to result in model predictions that agreed well with the cold-flow data. 

As shown in Figure III-5, the bleed flowrate for the "blow- 
off" condition matches the cold-flow data to within 1/2%, but the base pres- 
sure prediction becomes asymptotic after 3.5% bleed. At 3.5%, the plume 
impingement results in a 3% loss in stagnation pressure according to the 
oblique shock calculations. The Nash "recompression criteria" indicate the 
actual stagnation pressure loss to be approximately 0.4 of the 3%. So the 
situation that exists is that increasing the bleed flowrate from 3.5% to 5.5% 
results in a small change in stagnation pressure loss and in the corresponding 
base pressure. 
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COLD-FLOW TESTS 107-110 




Figure 1 1 1-5. Model Improvements Allowed More Accurate 
Predictions for Smaller Nozzle Spacing ' 




Figure I I 1-6. Model Predictions became less Accurate 
with Increasing Nozzle Spacing 



Figure III-7. Model Accuracy Improved by Changing 
Plume Scaling Factor 
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PLUME SCALING FACTOR (FCTR) 


1.4 


1.3 


1.2 


1.1 



NOZZLE SPACING (XLE) 


Figure III-8. Nozzle Spacing Influence on Plume Scaling Factor 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


The third empirical factor is the jet spreading parameter 
(sigma), which is used in the shear layer calculations and defines the rate at 
which the shear layer jets spread. Presently, the bleed flow model uses a 
correlation developed by Korst which states that the spreading parameter is a 
function of plume boundary Mach number. From the derivation of the shear 
layer calculations presented in Section III,B,1, it can be seen that the jet 
spreading parameter will influence the mass fraction Yp^ and recirculation 
temperature Tf/To* Calibration to the cold-flow data was not possible 
since the recirculation temperature was not measured during cold or hot-fire 
testi ng. 


4. Parametric Analysis 

A parametric analysis was performed to investigate the 
effects of varying bleed flow thermodynamic properties and primary nozzle 
geometry on bleed flow requirements. Changes were made in molecular weight 
(MW) and specific heat (Cp) for both primary and bleed flows, bleed flow to 
primary flow stagnation temperature ratio, nozzle axial spacing, primary 
nozzle area ratio, primary nozzle size and contour. 

Table III-II lists the effects of bleed flow thermodynamic 
properties on bleed flow requirements. The nozzle configuration simulated was 
that used for the dual throat cold-flow tests (NAS 8-32666). Varying the 
ratio of specific heats (vg) for the bleed flow appears to have little 
effect on the required bleed flowrate. Variations in bleed flow molecular 
weight and stagnation temperature appear to have significant effects on bleed 
flow requirements. Decreasing the molecular weight of the bleed gas results 
in a lower required flowrate, whereas a decrease in stagnation temperature 
results in an increase in the required bleed flow. An increase in the primary 
flow gamma (yp) will also result in requiring less bleed flow to turn the 
plume for a shock-free attachment. 
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TABLE III-II 


EFFECT OF BLEED FLOW PROPERTIES ON 
REQUIRED BLEED FLOWRATE 


■^P 

^B 

MWp 

MWb 

^OB^^O 

P 

y Ws 
* u~ 
”p 

Pb/Po 

p 

1.23 

1.4tf 

20 

29 

1.0 

11.93 

.0483 

1.23 

1.15 

20 

29 

1 .0 

11.54 

.0483 

1.23 

1.40 

20 

29 

0.5 

15.22 

.0483 

1.23 

1.40 

20 

15 

1.0 

8.60 

.0483 

1.40 

1.23 

29 

20 

1.0 

4.95 

.0310 

Yp 

- 

Primary Nozzle 

Flow Ratio of 

Specific 

Heat 


^B 

- 

Bleed Flow Ratio of Specific 

Heat 




MWp - Molecular Weight Primary Nozzle 

MWg - Molecular Weight Bleed Flow 

Tqb/Tq - Ratio of Stagnation Temperatures 

^ ~ Bleed Flowrate 

Pb/Pqp - Bleed Flow Pressure Ratio 

NOTE: Configuration Simulated was Dual Throat Cold-Flow Hardware 


47 



Ill, B, Aerodynamic Bleed Flow Model (cont.) 


The effect of geometric design parameters on bleed flow 
requirements was investigated by single parameter variations about a baseline 
configuration. A baseline configuration was generated from the operating spe- 
cifications presented in the final report of the Dual-Fuel, Dual Throat Engine 
Preliminary Analysis (NAS 8-32967). Table III-III lists the main geometric 
design parameters for this design, and Figure III-IO shows the baseline 
geometry and plume contours for various bleed flowrates. For the baseline 
design, the bleed flow model calculated that a 5.94% bleed flowrate was needed 
for the plume to attach shock-free to the secondary nozzle, i.e., the blowoff 
condition. Reducing the bleed flow from the blowoff condition to approxi- 
mately 3% showed that the resulting plume boundary pressure, plume contour, 
and shock strength were practically the same in the blowoff condition. 

The results of varying the nozzle spacing are presented in 
Table III-IV. Figure III-ll shows the configurations used and the resulting 
plume contours for 0% bleed flow. For this case, the plume attachment point 
moves further upstream as the nozzle spacing increases. The plume boundary 
pressure also increases with increased nozzle spacing, resulting in a 
decreased slope of the plume at the attachment point. Further, as the nozzle 
spacing increases, the bleed flow requirement for flow attachment at the 
secondary throat also increases. The net effect of these trends is that the 
total flow angle decreases with increased nozzle spacing, resulting in reduced 
shock strength. 


Table III-V lists the bleed flow requirements for varying the 
primary nozzle area ratio. As the primary nozzle area ratio increases, the 
bleed flow requirement decreases. For 0% bleed flow, increasing the primary 
nozzle area ratio causes the plume attachment point to move upstream and the 
plume boundary pressure to decrease. The decrease in plume boundary pressure 


48 



TABLE III-III 


DUAL THROAT BASELINE GEOMETRY 


Mode I Sea Level Thrust - 600,000 IbF 
Mode I to Mode II Thrust Ratio - 3.2 


SECONDARY NOZZLE: 


Chamber Pressure 2800 psia 

Mixture Ratio 2.8 

Area Ratio 45.6 

Contraction Ratio 2.5 

Throat Radius 6.41 in. 

Percent Mode I Thrust 70% 

Inlet Radius 6.41 in. 

Upstream Throat Radius 6.41 in. 

Downstream Throat Radius 3.25 in. 

Nozzle Length 123.7 in. 

PRIMARY NOZZLE : 

Chamber Pressure 4000 psia 

Mixture Ratio 7.0 

Area Ratio 2.7 

Contraction Ratio 2.5 

Throat Radius 2.90 in. 

Percent Mode I Thrust 30% 

Inlet Radius 2.90 in. 

Upstream Throat Radius 2.90 in. 

Downstream Throat Radius 1.80 in. 

Nozzle Length 6.3 in. 

Nozzle Spacing 2.75 in. 
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Figure III-IO. Dual Throat Baseline Configuration 


TABLE III-IV 

VARIATION IN NOZZLE SPACING 
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0 
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- 

- 

- 

- 

- 
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.06239 

2.41 
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.08187 

2.27 
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-8.0 
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* 
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2.23 

7.9 

2.4 
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* 
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F Ratio 
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XLE 
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V^o 


“Plume 


‘^Flow 

®Wall 

^S.L. 


tp 

F Ratio 
P«/P 


Spacing Between Primary Nozzle Exit and Secondary Nozzle Throat 
Bleed Flowrate 

Plume Boundary to Chamber Pressure Ratio 
Mach Number Along Plume Boundary 

Plume Angle at Attachment Point 
Total Flow Turning Angle 
Wall Angle At Attachment Point 
Shear Layer Thickness 
Primary Nozzle Area Ratio 
Primary Nozzle Throat Radius 
Mode I to Mode II Thrust Ratio 
Pressure Ratio for Primary Nozzle 
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TABLE III-V 


VARIATION IN PRIMARY NOZZLE AREA RATIO 



W7W 
s' p 

Pe/f’o 


"b 

“Plume 

“^Flow 

^^Wall 

*^S.L. 

1.62 

0 

.06569 


2.38 

23.6 

36.4 
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.24237 


3.0 

.09087 


2.21 

18.4 

14.4 
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- 
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2.707 
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3.0 
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- 
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‘'p 
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.04782 
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= 2.753 



^cp 

F Ratio 
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K 

tp 

= 2.90 



= 3.2 




XLE - Spacing Between Primary Nozzle Exit and Secondary Nozzle Throat 

W /W - Bleed Flowrate 

s P 

Pb/Pq - Plume Boundary to Chamber Pressure Ratio 

- Mach Number Along Plume Boundary 
Upiume - Plume Angle at Attachment Point 

Flow - Total Flow Turning Angle 

®Wall " Angle at Attachment Point 

- Shear Layer Thickness 

- Primary Nozzle Area Ratio 

R^p - Primary Nozzle Throat Radius 

F Ratio - Mode I to Mode II Thrust Ratio 

Pg/Pp - Pressure Ratio for Primary Nozzle 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 

is due to the decreasing nozzle exit pressure from the increase in nozzle area 
ratio. Figure III-12 shows the resulting plume contours for 0% bleed flow. 

Changing the primary nozzle size was accomplished by keeping 
the same primary nozzle area ratio but changing the throat area. A change in 
throat area corresponds to a change in the Mode I to Mode II thrust ratio. 
Table III-VI lists the change in bleed flow requirements for change in throat 
radius. As the throat radius increases, decreasing thrust ratio, the bleed 
flow requirement increases slightly. For the 0% bleed flow case, the plume 
attachment point moves upstream, and the plume boundary pressure increases 
with increasing throat radius. Figure III-13 shows the plume contours for the 
0% bleed flow case. 

The effect of primary expansion nozzle contour was evaluated 
by comparing aerodynamic model predictions for the baseline Bell contour noz- 
zle with an exit angle of 16® to results calculated for a 10° conical nozzle. 
The low exit angle conical nozzle is predicted to require less bleed flow 
(4.5% vs 5.9%) which subsequently results in a lower base pressure ratio 
(0.0717 to 0.05489). These decreased values are due to the change of nozzle 
exit angle which, when combined with the Prandtl-Meyer expansion angle at the 
nozzle lip, results in a smaller plume expansion angle. Having a plume bound- 
ary with a shallow angle requires less bleed flow and a smaller base pressure 
to turn the plume and allow it to attach tangentially to the secondary nozzle 
wall. 

In summary, the dual throat aerodynamic bleed flow model 
developed during NAS 8-32666 has been updated to include a procedure for pre- 
dicting the base pressure generated in Mode II operation with zero bleed flow 
and provide for dissimilar exhaust gases to be used for the primary and bleed 
flows. These changes have enhanced the analysis capability of the Mode II 
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TABLE III-VI 


VARIATION IN PRIMARY NOZZLE SIZE 
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Ill, B, Aerodynamic Bleed Flow Model (cont.) 


dual throat operation and provide a design tool to permit tradeoffs of various 
geometric and operating parameters. 

The improved model has been used to evaluate bleed flow 
requirements as a function of both bleed flow properties and geometric design 
parameters. These results are summarized in Table III-VII. This model, 
together with other dual -throat design and analysis improvements, has also 
been utilized to assess changes in predicted performance of the dual -throat 
concept as discussed in later sections of this report. The bleed flow model, 
with suitable modifications, also provided the basis for prediction of the 
base flow region of the dual expander concept as described in Section V of 
this report. 


C. MODE II BOUNDARY LAYER LOSS VERIFICATION 
1. Previous Simplified Approach 

During Contract NAS 8-32666, the Mode II boundary layer loss 
was obtained by using the turbulent boundary layer chart procedures outlined 
in Appendix B of CPIA No. 178 with some additional considerations to the aero- 
dynamics associated with Mode II operation. The primary assumptions for this 
calculation are that the boundary thrust decrement is additive and propor- 
tional to the momentum thickness. The boundary layer loss was obtained by 
using the TBL chart procedures for a conventional Rao nozzle having an area 
ratio defined by the point of plume attachment and the primary throat. A 
ratio of the shear layer momentum thickness (at the point of attachment given 
by the aerodynamic bleed flow program) to the momentum thickness given by the 
TBL chart procedures was used to proportion the conventional nozzle boundary 
layer loss to approximate the actual thrust layer loss obtained for the nozzle 
from the attachment point to the secondary nozzle exit, yielding the total 
Mode II boundary layer loss. This procedure is outlined in Figure III-14. 



TABLE III-VII 


BLEED FLOW THERMODYNAMIC PROPERTIES EFFECTS 


PROPERTY 


Bleed Flow Yg 
Temperature Ratio 



MW Bleed Flow 


EFFECT ON BLEED FLOW REQUIREMENT 
Not Significant 

Increases with Decreasing 
Decreases with Decreasing MWg 



GEOMETRIC PARAMETER EFFECTS 


PROPERTY 

Nozzle Spacing (XLE) 

Primary Nozzle Area Ratio {c^) 


EFFECT ON BLEED FLOW REQUIREMENT 

Decreases with Decreasing XLE 

Increases with Decreasing e 

a p 


Primary Nozzle Size (R. ) 

"cp 


Decreases with Decreasing R 


tp 


Primary Lip Angle (THEp) 


Decreases with Decreasing THE 


P 
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Figure 111-14. NAS 8-32666 Boundary Layer Loss Methodology 



Ill, C, Mode II Boundary Layer Loss Verification (cont.) 


These calculations were made by assuming an adiabatic wall to simplify the 
calculation of loss due to heat extraction from the boundary as well as the 
gain due to the addition of this heat to the incoming regenerative coolant 
fuel . 


2. BLIMP Analysis 

The BLIMP computer program was used to provide a rigorous 
analysis of the boundary layer loss during Mode II operation for comparison to 
the simplified approach described above. BLIMP was successfully run using the 
Kendall and binary diffusion options for the Mode II dual throat baseline con- 
figuration with full bleed flow. These options are consistent with the stan- 
dardized JANNAF performance methodology. To accomplish this, it was necessary 
to modify the aerodynamic bleed flow model to generate additional points used 
in the spline fit to describe the effective Mode II "nozzle contour" and to 
calculate the momentum thickness at the plume attachment point, including the 
contribution due to the bleed flow. Also, the TDK program was modified to 
accept the increased spline fit size needed to calculate the pressure profile 
to be used in BLIMP. 

The BLIMP run was started at the plume attachment point (Sta- 
tion 1) and continued throughout the secondary nozzle. The BLIMP input param- 
eter S(l), which represents the length of a flat plate required to generate 
the desired momentum thickness, was varied until the momentum thickness at 
Station 1 equaled the momentum thickness calculated by the bleed flow program. 
In this manner, the conditions for the first calculation station within the 
BLIMP solution are representative (at least in terms of momentum deficit) of 
the conditions calculated for the shear layer using the bleed flow model. In 
other words, the BLIMP analysis was started with initial boundary layer pro- 
files which are representative of the boundary layer growth within the primary 
nozzle and shear layer, including the effects of bleed flow addition. 
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Ill, C, Mode II Boundary Layer Loss Verification (cont.) 


Starting with these initial conditions, the BLIMP computer program calculated 
the additional boundary layer growth in the secondary nozzle downstream of the 
plume attachment to the nozzle exit. Standard BLIMP parameters at the final 
station were then used to calculate the net boundary layer loss, as shown in 
Table III-VIII. 


Table III-VIII also shows a comparison between the boundary 
layer loss calculated from the BLIMP analysis and that determined by using the 
simplified procedure from NAS 8-32666. As can be seen, the net loss based on 
the BLIMP results is more than twice as great as the loss calculated by using 
the simplified procedure. Since the BLIMP solution is by far the more rigor- 
ous procedure, it is considered the best available current estimate for Mode 
II boundary layer losses. It should be noted, however, that even the net 
boundary layer loss from the BLIMP solution includes the performance loss 
resulting from the addition of the bleed flow. 

As a result of the poor comparison between the simplified and 
rigorous boundary layer loss procedures, a new simplified procedure was devel- 
oped in order to provide loss estimates close to those calculated from the 
BLIMP procedure without extensive and detailed analysis. While the revised 
procedure still uses the TBL chart as a basis for the loss calculation, it 
scales the results to account for shear layer effects in a different manner 
than the previous method. 

Briefly, the TBL chart method is used to calculate the momen- 
tum thickness for a conventional nozzle expansion from the primary nozzle to 
the plume attachment point. The momentum thickness calculated from the shear 
model is then divided by this conventional momentum thickness to define a 
parameter, RTHETA, which is used in the simplified procedure to account for 
upstream influences on the boundary layer. Thus, 
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TABLE III-VIII 


COMPARISON OF BOUNDARY LAYER LOSS CALCULATIONS 


Mode II Baseline Configuration (See Table III-III) 
% Bleed Flow = 5.3 


Boundary Layer Performance 
Parameters 

Total Flowrate (Ibm/sec) 

Gross Thrust Loss (IbF) 

Gross Isp Loss (sec) 

Total Heat Loss (Btu/sec) 

Isp Gain From Regen Heat (sec) 

Net BLL (sec) 

^^^Based on TW/TS =1.0 RTHETA = 

Based on TW/TS =0.2 RTHETA = 

^^^Based on Difference Between Boundary 


BOUNDARY LAYER ANALYSIS TECHNIQUE 


BLIMP 

NAS 8-32666 
Simpl ified 
Procedure 

Revised 

Procedure 

483.3 

483.3 

483.3 

8095 

3335 

7974 

16.8 


16.5(2) 

37,050 

0 

N/A 

3.3 

0 

3.4^^) 

13.5 

6.9 

13.1 


1.00 

2.22 

Layer Loss for TW/TS = 0.2 & TW/TS = 1.0 


= Wall Temperature 

Stagnation Temperature 
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Ill, C, Mode II Boundary Layer Loss Verification (cont.) 


RTHETA 


®shear‘^^TBL 
model chart 


Plume attach point 


Previously defined trends, such as those shown in Figure III-15, are used to 
estimate the change in the RTHETA value for further expansion in the secondary 
nozzle section. This final RTHETA value (i.e., at the secondary nozzle exit) 
is used to ratio the momentum thickness calculated for a conventional nozzle 
expansion from the primary nozzle throat to the secondary nozzle exit, i.e.. 


^Exit ' »Exit, TBL "™ETA 


The boundary layer calculations are made by using the actual 
wall temperature conditions. In the event of regenerative cooling, the new 
approach also takes into account the increase in the IspoDE 
increased propellant enthalpy transferred to the regenerative coolant. 


This approach was incorporated into a computer program called 
DT*BLL0SS. The inputs required for the program are 1) characteristic velocity 
(Cp*), 2) chamber pressure (Pep), 3) throat radius (Rtp), 4) ratio of 
specific heats (yp), 5) Mode II area ratio (£eng)» ^eng» 

7) area ratio at plume attachment point (E-rlm), and 8) momentum thickness 
(0Bieed) at the plume attachment point. The first six parameters are 
obtained from ALRC's Dual Throat Geometry and Performance Prediction Program 
(FD 0169). The last two parameters come from the Aerodynamic Bleed Flow 
Program (FD 0170). 

The results of this "revised procedure" for the baseline case 
are compared to both the BLIMP rigorous and previous simplified cases in Table 
III-VIII. As noted, the new procedure provides loss predictions which are 
similar to the BLIMP calculations. The impact of this new procedure on the 
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Ill, C, Mode II Boundary Layer Loss Verification (cont.) 


performance predictions for Mode II dual throat operation is identified quan- 
titatively in Section III,F of this report. 

D. SHOCK MODEL 

The objective of this subtask was to develop a method for 
obtaining two-dimensional axially symmetric flow field solutions that include 
the effects of shocks to estimate shock effects on performance and thermal 
operating characteristics. Mode II operation at design conditions will neces- 
sarily contain shock waves, the effect of which may be negligible at low area 
ratios but which are significant at high area Mode II conditions. 

Since shock wave development in a nozzle is fundamentally a two- 
dimensional problem, a two-dimensional computer analysis must be used to pre- 
dict the effect of shocks on performance. The computer code selected for this 
purpose is the VNAP code (Ref. 4). It has been modified for use in calcu- 
lating both the flow field within the dual throat nozzle during Mode II opera- 
tion and the resulting performance. 

The VNAP program was selected for performing the dual throat noz- 
zle calculations because it is the only computer code readily available that 
is capable of treating supersonic nozzle flow with imbedded shocks in a gen- 
eral manner. It is capable of calculating viscous as well as inviscid steady 
and unsteady internal flow. The approach used is to solve the Navier-Stokes 
equations for two-dimensional, time-dependent, compressible flow. The fluid 
is assumed to be a perfect gas. An explicit artificial viscosity is used for 
the shock calculations. The steady-state solution is obtained as the asympto- 
tic solution for large time. The nozzle walls are transformed to a rectangu- 
lar grid so that the resulting partial differential equations can be solved by 
using the second-order MacCormack finite-difference scheme. 
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Ill, D, Shock Model (cont.) 


Modifications to VNAP were required to add a thrust calculation 
for obtaining thrust chamber performance and converting the plotting sub- 
routines for use with the Cal-Comp plotter. Three cases were investigated: 

(1) A shock-free nozzle configuration. Output from this case was 
then compared with output from TDK to verify and calibrate 
the accuracy of the VNAP code. 

(2) A baseline design configuration where bleed flow is used to 
control the plume attachment. 

(3) A baseline design configuration with zero bleed flow. Output 
from this case was used to define the freestream boundary 
conditions after the shock. 

Previously, VNAP had been applied mainly to transonic analyses. 
This study was the first time that the program was used to investigate super- 
sonic flow fields such as those occurring in large area ratio nozzles. Sys- 
tematic errors were expected because of the large solution grid used in the 
finite-difference approach. 

To determine the magnitude of these errors, VNAP was used to ana- 
lyze the flow field for a conventional shock-free nozzle with the same Mode II 
area ratio and envelope as that of the dual-throat baseline. A comparison of 
the performance calculated by VNAP with that calculated by TDK for the conven- 
tional nozzle provided a performance correction factor which presumably could 
be applied to the VNAP solutions obtained with shock flow. The results of 
this calibration are as follows: 

^^^VNAP ” 484.2 sec 

Isp.|.pi^ = 473.5 

^Isp corr ~ _ _gyg 

484.2 
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Ill, D, Shock Model (cont.) 


Therefore, it appears that the VNAP calculation results in a performance which 
is approximately 2% greater than an accurate method of characteristic solution 
(TDK) and that the .978 correction factor should be utilized for the two addi- 
tional VNAP cases which contain flow shocks. 

Plots were generated showing velocity, Mach number, pressure, tem- 
perature, and density profiles to help visualize the flow field situation. 
Examples of these plots for the conventional nozzle analysis are shown in 
Figures III-16 through III-21.’ To help facilitate the computer storage 
requirement, the nozzle's flow fields were arbitrarily divided into two 
segments, referred to as upstream and downstream on the referenced figures. 

The VNAP simulation of the dual throat baseline Mode II configura- 
tion predicted a nozzle performance of 480.2 seconds. This was corrected to 
469.6 seconds by using the correction factor obtained from the conventional 
nozzle simulation. The resulting Mode II nozzle divergence efficiency was 
found to be 0.988, compared to a conventional nozzle efficiency of .996. 

Thus, the divergence efficiency of the dual throat Mode II expansion is 
approximately 0.8% less than that of a conventional nozzle expanding to the 
same area ratio. This is only 0.3% greater than the value used during NAS 
8-32666. 


Plots of the flow field (Figures III-22 to III-27) show a right 
running shock wave forming in the nozzle due to the compression caused by the 
highly contoured secondary nozzle wall. This shock structure, however, never 
reaches the downstream secondary wall and thus has only a slight influence on 
performance and heat transfer. 

The VNAP program provided a means for determining the gas-side 
boundary conditions used in the thermal analysis of the secondary throat 
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Figure III-16. Conventional Nozzle Upstream Velocity Vector Plot 
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VELOCITY VECTORS (IX) 


N= 601 T=*.1198E-04 SEC 



Figure lH-17. Conventional Nozzle Downstream Velocity Vector Plot 
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Figure III-18. Conventional Nozzle Upstream Mach Number Plot 
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Figure II I -19. Conventional Nozzle Downstream Mach Number Plot 
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Figure III-20. Conventional Nozzle Upstream Pressure Plot 
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PRESSURE N= 601 T=*.1198E-04 SEC 

LOW VRLUE=33.5293E-02 HIGH VflLUE=45 .3076E+0Q 

LOWC0NTOUR=. 4.8325E+QQ HIGH CQNTQUR= 4.0810E+D1 DELTA CONTOURz 4.4972E+D0 


Figure III-21. Conventional Nozzle Downstream Pressure Plot 
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Figure I I 1-24. Dual Throat Baseline Upstream Mach Number Plot 




MflCH NUMBER N= 601 T=«.3S29E-0H SEC 

LOW VflLUE=m .2089E-01 HIGH VPILUE=56.2791E-01 

LOHCONTQUR=. 4.2716E+00 HIGH CONTOUR^ 5.^^72E+00 DELTR CONTOUR= 1.5D70E-01 



Figure IH-25. Dual Throat Baseline Downstream Mach Number Plot 




PRESSURE N= 601 T=*.3059E-04 SEC 

LOW VflLUE=73.6069E-02 HIGH VflLUE=30 .0390E+01 

L0WC0NT0UR=. 3.0701E^-01 HIGH CONTOUR^ 2.70H2E+02 DELTR CONTOUR^ 2.99B5E+D1 



Figure lH-26. Dual Throat Baseline Upstream Pressure Plot 




Figure ill -21 . Dual Throat Baseline Downstream Pressure Plot 



Ill, D, Shock Model (cont.) 


region for the case having zero bleed flow (the case expected to be the most 
adverse Mode II operating condition). Pressures and mass velocities v/ere 
obtained at different axial locations and used as boundary conditions for the 
design thermal analysis (discussed in Section IV). 

The flow field plots (Figures III-28 through III-30) show the for- 
mation of the plume attachment shock. The plots are indicative of a weak 
shock formation. This can also be seen in the wall pressure profile shown in 
Figure 1 11-31. The static pressure rise is approximately 15% larger than the 
value calculated by the aerodynamic bleed flow model, but it is approximately 
45% of what would be expected with a discrete shock. This value appears to be 
consistent with the Nash factor value of 0.4 chosen during calibration of the 
bleed flow model, as described in Section III,B,2. 

E. SECONDARY NOZZLE CONTOUR OPTIMIZATION 

During Mode I operation, the dual throat nozzle operates as a con- 
ventional nozzle. The secondary nozzle contour is designed by using the opti- 
mization technique developed by G. V. R. Rao. This technique, which is 
derived from the calculus of variations and employs the method of characteris- 
tics, is used in ALRC's nozzle design program (FD 0002) and will calculate 
contours to provide maximum thrust for a specified length and/or area ratio. 
The Rao method is applicable when, together with an initial wall contour, a 
shock-free supersonic data line is given across the nozzle. In ALRC's pro- 
gram, this initial data line is automatically provided by using the modified 
Hall method to approximate the solution of the transonic flow equations in the 
nozzle throat region. In general, however, the method can be applied to any 
supersonic approach flow that is shock-free. This fact has been utilized in 
the study presented here to provide the design computer program with the 
capability of designing optimum nozzle contours for dual throat thrusters. 

The method used to design nozzles for Mode II operation is described below. 
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Figure Hl-29. Dual Throat Baseline Zero Bleed Flow Pressure Profile 
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Figure III-30. Dual Throat Baseline Zero Bleed Flow Mach Number Profile 
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Ill, E, Secondary Nozzle Contour Optimization (cont.) 

1. Design Method 

The aerodynamic model presented in Section III,B can be used 
to design a dual -throat thruster assembly which, by use of bleed flow, oper- 
ates with a plume shape and attachment position such that the flow field is 
nearly shock-free. Once the approach flow is specified up to the plume 
attachment position, the Rao procedure can be used to calculate an optimal 
contour for the exterior wall downstream of the attachment position. This 
wall contour is typically quite long because the thruster mission requires a 
high expansion ratio. 

The design computer program has been modified to allow an 
optimum nozzle contour to be calculated where the approach flow is determined 
by an input nozzle contour. This feature of the program is called the 
"Extended Nozzle Wall Option." The steps used to implement this option are as 
follows: 

First, the v/all geometry subroutines of the TDK method of 
characteristics program were transferred to the ALRC program and checked out. 
These are subroutines WALL, CUBIC, and XSLP. The maximum number of points 
allowed for inputting the nozzle wall using the spline fit option (IWALL = 4) 
was increased from 20 to 40. This same modification was also made to the 
aerodynamic model computer program. 

Next, subroutines WALLl and WALL2 were revised to allow the 
circular arc portion of the wall contour to be centered at a predetermined 
position. Previously it had been necessary to position the circular arc at 
the nozzle throat and center it at (1, 0). A new wall point MOC subroutine, 
called WALLPT, was also written and checked out. 
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Ill, E, Secondary Nozzle Contour Optimization (cont.) 


A logic subroutine, called RRC, was written to construct a 
nozzle flow field for a given nozzle wall geometry (i.e., as constructed by 
subroutine WALL) using Right Running Characteristics (RRC's). The flow field 
constructed (shown in Figure 1 1 1-32), is bounded as follows: 

(a) The upstream boundary is the initial data line RRC. 

(b) The upper boundary is the nozzle wall contour. 

(c) The lower boundary is the flow axis. 

(d) The downstream boundary is the RRC extending from the 
wall end point. 

Finally, the main program, FD002, and the logic program, 
LOGIC, were modified to accept the above. The program was checked out for the 
baseline case for the dual throat nozzle. That is, an optimum nozzle contour 
of fixed length and area ratio was designed for the secondary nozzle wall, 
assuming the approach flow to be that given by the primary nozzle and plume 
(aerodynamic wall) with 5.3% bleed flow. 

2. Optimization Analysis Results 

A parametric analysis was performed to investigate the effect 
of optimizing the secondary nozzle contour for either Mode I or Mode II. 

Table III-IX and Figure III-33 show the results of this analysis, including a 
comparison of the Mode II nozzle contours according to the following param- 
eters: 



Ill, E, Secondary Nozzle Contour Optimization (cont.) 

Contour (1) is a conventional Bell nozzle with an expansion 
ratio of 222:1. The chamber throat and exit 
point of this nozzle are those of the primary 
nozzle of the dual throat baseline contour. 

Contour (2) results from optimizing the secondary nozzle for 
Mode I operation. The plume attaches to a 
secondary nozzle contour designed for an 46:1 
area ratio in Mode I operation. 

Contour (3) has been designed in terms of "Extended Nozzle 
Wall Option" discussed in this section. The 
approach flow passes through the dual throat 
primary nozzle and the plume expansion, subse- 
quently attaching at the secondary nozzle 
throat. Contour (3) is the optimal contour 
given this approach flow. This nozzle has an 
expansion ratio of 222:1. 

Contours (1) and (3) have the same throat and exit point, 
but, as can be seen from Figure 1 11-33, they are very different in shape. 
Contours (2) and (3) have different throat locations but the same exit point. 
Contour (2) is designed to yield maximum thrust during Mode I operation. 
Contour (3) is designed to yield maximum thrust during Mode II operation. As 
can be seen from Figure 1 1 1-33, Contour (3) has less curvature than Contour 
(2) in order to turn the high Mach number approach flow isentropically. The 
net performance effect of the two optimization options is shown in Table 
III-IX. 
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Figure III-33. Dual Throat Secondary Nozzle Contour Comparison | 



Ill, E, Secondary Nozzle Contour Optimization (cont.) 


TABLE III-IX 

COMPARISON OF NOZZLE DIVERGENCE EFFICIENCIES FOR 
TWO TYPES OF OPTIMUM NOZZLE CONTOURS 


Operating 

Mode: 

Mode I 
e= 45:1 

Mode II 
6= 222:1 

Secondary Nozzle 

MODE I 

MODE II 

MODE I 

MODE II 

Optimizied For: 

(Contour 2) 

(Contour 3) 

(Contour 2) 

(Contour 3) 

'’^DIV 

.991 

.988 

.982 

.985 


Source: TOE Option TDK Program 


The absolute values of the divergence efficiencies shown in 
Table III-IX are somewhat questionable because of the simplified treatment of 
the flow recompression used in the TDE analysis. The relative efficiency 
values appear to be valid, however. Therefore, optimizing the secondary noz- 
zle contour for Mode I or Mode II will only slightly change the performance. 

F. DUAL THROAT PERFORMANCE 

Many of the model improvements described previously have a direct 
impact on the performance prediction methodology for dual throat engines. 
Improvements in the aerodynamic model have resulted in less restrictive calcu- 
lational assumptions, tying it directly to the shear model and overall bound- 
ary layer model via the evaluation of the momentum thickness at the plume 
attachment point. Results of detailed two-dimensional flow analyses (VNAP and 
TDE) have further verified the fluid conditions within the dual throat 
thruster and the resulting performance efficiencies. 
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Ill, F, Dual Throat Performance (cont.) 

Presented within this section are the results of an assessment in 
the changes in the predicted dual throat performance based on the use of the 
improved aerodynamic model and updated simplified performance model. First, 
an analysis of the baseline design condition for Mode II is presented and com- 
pared to a conventional gas-generator cycle engine. Secondly, the performance 
of a limited set of parametric cases calculated during NAS 8-32967 is compared 
to values calculated by using the updated simplified performance model. 

1. Baseline Design Performance 

The baseline operating point, shown in Table III-X, is the 
same as that recommended in Contract NAS 8-32967. For this design condition, 
a bleed flowrate of 5.3% is predicted to be required for shock-free (zero 
pressure loss) plume attachment. However, the bleed flow can be reduced to 
3.5% with only a small loss in total pressure (^1%), and thus this test condi- 
tion was selected as a reasonable design tradeoff for this comparison. 

The calculated specific impulse for Mode II operation is 
shown in Table III-XI and compared to a conventional gas-generator cycle 
engine operating at the same condition. The results of this comparison show 
that the performance of the dual throat engine in Mode II is approximately 6 
seconds (1-1/4%) lower than that of the conventional engine. As noted, most 
of this difference in performance is attributed to the less efficient expan- 
sion process of the dual throat engine in Mode II which results in a greater 
divergence efficiency. This performance penalty must, of course, be traded 
off with other dual throat advantages as a result of its unique capabilities. 

2. Comparison of Simplified Performance Parametrics 

The effect of the previously described model improvements on 
performance parametric trends was briefly investigated by using the updated 
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TABLE III-X 


DUAL THROAT BASELINE OPERATING POINT 


MODE I SEA LEVEL THRUST 2.669 MN (600K IbF) 
MODE I TO MODE II THRUST RATIO 3.2 


PRIMARY NOZZLE 

CHAMBER PRESSURE 
MIXTURE RATIO 
PROPELLANT 

SECONDARY NOZZLE 

CHAMBER PRESSURE 
MIXTURE RATIO 
PROPELLANT 


27.58 MPa (4000 psia) 
7.0 

LOX/LH2 


19.31 MPa (2800 psia) 
2.8 

LOX/RP-1 
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TABLE III-XI 


COMPARISON OF THE DUAL THROAT AND CONVENTIONAL 
GAS-GENERATOR ENGINE CONCEPTS 



PROPELLANTS: 

LO 2 /LH 2 



MIXTURE RATIO: 

6.0 (Overall) 



CHAMBER PRESSURE: 

27.58 MPa (4000 psia) 


AREA RATIO: 

222:1 



THRUST (VAC): 

956 KN (215K IbF) 



TOTAL FLOWRATE: 

215.6 Kg/sec (475 Ibm/sec) 

DUAL THROAT 

- MODE II 

CONVENTIONAL 

GG CYCLE 

BLEED FLOW 

= 3.5% 

GG FLOW = . 

3.5% 

PRIMARY 0/F 

= 7.0 

TCA 0/F = 

7.0 

Isp ODE 

= 480.4 sec. 

^SP ODE 

480.4 sec 

’^DIV 

= .988 

'’DIV 

.996 

m 

m 

= .990 

^ERE 

.990 

\IN 

= .999 

'^KIN 

.999 

^^SP BL 

= 13.4 sec. 

^ ^SP BL 

5.7 sec. 

^SP TCA 

= 456.2 sec. 

^SP TCA 

467.5 sec 

^SP ENG 

= 456.2 sec. 

^SP ENG 

462.2 sec 
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Ill, F, Dual Throat Performance (cont.) 


dual throat simplified model. The results for five parametric cases examined 
during NAS 8-329(i7 were compared with the present model predictions for both 
Mode I and II operation. These cases ucre used to investigate tfie effects of 

chamber pressure and thrust ratio (Fv^c^ Mode l/l^vac Mode 2 ) engine 
performance. 


In Mode I, the vacuum specific impulse predictions and trends 
are essentially unchanged from the previously reported results. However, 
because the cal culational technique was changed to require equal strearntube 
static pressure at the secondary throat during Mode 1 operation, a significant 
difference in the sea level specific impulse trend with thrust ratio v/as 
noted. This is illustrated in Figure III-34. The current model predicts less 
variation in the sea level specific impulse at the extremes of the thrust 
ratio. Note that at the baseline thrust ratio of 3.2, the present model pre- 
dicts a sea level specific impulse which is approximately 6 sec greater tfian 
the previous model predictions. Engine geometry trends with thrust ratio were 
also impacted, as noted in Figurfe HI-35. The secondary throat size, rather 
than increase as previously predicted, is now expected to actually decrease 
slightly with thrust ratio (while holding sea level thrust constant). 

> In all cases, the Mode II S|)Gcific imi)ulse is lower with use 

of the updated performance model because of larger boundary layer and diver- 
gence losses. The magnitude of this difference varies from case to case, but 
is generally 1 to 2% less. While the magnitude of the Mode II performance is 
less than previously predicted, the certainty of the prediction is much 
greater as a result of the work completed during this study. In addition, as 
noted, the Mode II performance predicted for the recommended baseline is only 
1-1/4% less than that of a conventional gas-generator engine, and thus the 
previous simplified model predictions appear optimistic in this light. 
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MODE I DELIVERED Isp, SECONDS 






RADIUS, INCHES 


MODE I SEA LEVEL THRUST = 2.669 MN (600K IbF) 
Pc^ = 1448 N/cm^ (2100 psi) 

Pc = 2068 N/cm2 (3000 psi) 

MRg = 2.8 (RPl) 



j- SECONDARY 
\ PRIMARY 


Figure III-35. Dual Throat Engine Throat Radius as a Function of Thrust Ratio 


Ill, Dual Throat Aerodynamic Model and Performance Prediction 
Improvement (cent.) 


G. RECOMMENDATIONS 


1. Design Approach 


The previous methodology determines primary nozzle area ratio 
by the following correlation: 



A. EPPM 


where EPPM is an empirical factor. A value of EPPM =1.8 was chosen on the 
basis of the cold-flow test results from NAS 8-32666. Nozzle spacing was 
chosen to be 0.8, the length of the conical portion of the secondary conver- 
gence section. 


The primary nozzle area ratio and axial spacing between the 
primary and secondary throats should be chosen to ensure that the secondary 
flow is choked at the secondary throat, rather than at the primary lip, by 
allowing enough flow area on the secondary side of the primary nozzle lip. 

This design approach requires small primary nozzle area ratios and large noz- 
zle spacing. Secondly, the nozzle spacing chosen must result in the ambient 
pressure at the primary nozzle lip being low enough to cause the primary flow 
to be choked. This implies that the nozzle spacing should be small. Analyses 
of bleed flowrate trends for Mode II indicate that bleed flow requirements 
decrease with increasing primary nozzle area ratio and decreasing nozzle 
spacing. These trends are shown in Figure III-36. The primary nozzle area 
ratio and nozzle spacing chosen must be a compromise of all these require- 
ments. 
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Ill, G, Recommendations (cont.) 


2. Test Programs 

The aerodynamic bleed flow model has three empirical param- 
eters (identified in Section III,B). Two of these parameters have been cali- 
brated: the plume scaling factor and the NASH factor, using cold-flow and 

NASA/MSFC test data. The third parameter, the spreading factor, still needs 
to be calibrated. This parameter will influence the value of the mass frac- 
tion of primary flow entrained in the recirculation region and the resulting 
recirculation temperature calculated. This may impact the cooling require- 
ments in the plume attachment region. A test design is required that allows 
thermal data acquisition in these regions. This may be done either through a 
calorimetric convergence section or placement of thermocouples. 

3. Model Improvements 

The aerodynamic bleed flow model computer program currently 
determines the plume boundary shape by using a correlation of the form sug- 
gested by Herron. The intersection of the plume boundary and the secondary 
nozzle wall is determined and, together with the slope of these two surfaces, 
used to calculate recompression criteria for the shear layer. Because a 
correlation is used, the plume boundary is limited in range and is necessarily 
somewhat inaccurate. The effect of specific heat ratio is grossly simplified. 
The point of intersection and the slope of the plume and wall at the inter- 
section are also in error. The error in streamline slope change, i.e., the 
deflection angle of the plume streamline, can be more than 20%. 

Figure III-37a shov/s isobars calculated by the MOC for a 
cold-flow nozzle configuration. The streamline shape was calculated by the 
correlation discussed above. The plume streamline should be constant pres- 
sure, but is actually correct to within only about 15%. (The isobar labeled 
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Ill, G, Recommendations (cont.) 

X approximately follows the plume streamline.) Figure III-37b shows charac- 
teristic lines calculated by the MOC interior to the nozzle and plume boundary 
streamlines. In Figure III-37c, the flow field has been divided into two 
regions that are separated by the RRC that comes down from the primary nozzle 
lip. Region I, upstream of this RRC, is independent of the plume shape. 

Region I can be calculated by using the TDK Program and then input to the 
Bleed Flow program. 

Next, Region II can be calculated and used to replace the 
plume correlation. It is defined as follows: 

Upstream boundary: RCC ending Region I 

Left-hand corner: Prandtl -Meyer fan expanded to the 

base pressure 

Upper boundary: Constant pressure boundary (= base 

pressure) 

Downstream boundary: LRC from the upstream boundary to 

the plume attachment point on the 
secondary nozzle 

It would be necessary to calculate Region II using the MOD 
within the inner search procedure of the Bleed Flow program, i.e., directly 
replace the plume correlation with a MOC calculation. As can be seen from 
Figure III-37b, however, it is not necessary to calculate many points (only 
perhaps 300). The advantage would be that the MOC solution will be highly 
accurate. Thus, it is recommended that the aerodynamic model be modified to 
incorporate a MOC calcualtion for the plume boundary. 
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IV. 


DUAL EXPANDER ANALYSIS 


A. PRELIMINARY GEOMETRIC ANALYSIS 


This study was the first effort in technically analyzing the dual 
expander nozzle concept. The analysis deteniiined the dependence of perform- 
ance, throat area, nozzle and steamtube area ratios, mass flow, and nozzle 
length with the variation in Mode I sea level thrust, thrust ratio (Mode 1/ 
Mode II), primary and secondary chamber pressures, chamber pressure ratio 
(Pcs/Pcp). aiid nozzle exit pressure. The parametric range analyzed was as 
follows: 


Thrust 

Thrust Ratio 

Chamber 

Pressure 

Pc Ratio 

Exit Pressure 


2.22 - 8.90 MN 
1.0 - 5.0 
6.89 - 55.16 MPa 

0.25 - 1.0 
34.5 - 101 KPa 


(5 X 105 - 2 X 1Q6 IbF) 
(1000 - 8000 psia) 

(5.0 - 14.7 psia) 


To ensure a valid comparison with the dual throat analysis, LOX/ 
RP-1 was chosen as the primary nozzle propellant, and LOX/LH 2 was chosen for 
the secondary nozzle. Mixture ratios used were 2.8 and 7.0 for the primary 
and secondary chambers, respectively. 


Table IV-I lists the parametric values used for each case. The 
resulting performance, area ratios, throat areas, flowrates, and efficiencies 
are listed in Table IV-II. 


Mode I performance calculations ivere based on the simplified 
JANNAF methodology and vi/ere similar to the streamtube approach used during the 
Dual Throat Geometry Program (FD 0169). Mode II performance was determined by 
applying the simplified JANNAF methodology to a Bell nozzle having the same 
throat area, area ratio, and length as the secondary nozzle during Mode II 
operation. The procedures were programmed into a Dual Expander Preliminary 
Geometry and Performance Prediction Program. 
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TABLE IV-I 


DUAL EXPANDER CASES ANALYZED 


CASE 

SEA LEVEL 
THRUST 
MN (IbF) 

THRUST 

RATIO 

p 

^CP 

MPa (psia) 

p 

^CS 

MPa (psia) 

p 

*^EXIT 
KPa (psia) 

Pcs/Pi 

1 

2.67 (6 

X 

10^) 

2.37 

41.37 (6000) 

20.68 (3000) 

34.5 (5.0) 

.5 

2 

2.67 (6 

X 

10^) 

2.37 

41.37 (6000) 

10.34 (1500) 

34.5 (5.0) 

.25 

3 

2.67 (6 

X 

10^) 

2.37 

41.37 (6000) 

31.03 (4500) 

34.5 (5.0) 

.75 

4 

2.67 (6 

X 

10^) 

2.37 

41.37 (6000) 

41.37 (6000) 

34.5 (5.0) 

1.0 

5 

2.67 (6 

X 

10^) 

1.25 

41.37 (6000) 

20.68 (3000) 

34.5 (5.0) 

.5 

6 

2.67 (6 

X 

10^) 

3.50 

41.37 (6000) 

20.68 (3000) 

34.5 (5.0) 

.5 

7 

2.67 (6 

X 

10^) 

4.50 

41.37 (6000) 

20.68 (3000) 

34.5 (5,0) 

.5 

8 

2.67 (6 

X 

10^) 

5.0 

41.37 (6000) 

20.68 (3000) 

34.5 (5.0) 

.5 

9 

2.67 (6 

X 

10^) 

2.37 

13.79 (2000) 

6.89 (1000) 

34.5 (5.0) 

.5 

10 

2.67 (6 

X 

10^) 

2.37 

27.58 (4000) 

13.79 (2000) 

34.5 (5.0) 

.5 

11 

2.67 (6 

X 

10^) 

2.37 

55.16 (8000) 

27.58 (4000) 

34.5 (5.0) 

.5 

12 

2.22 (5 

X 

10^) 

2.37 

41.37 (6000) 

20.68 (3000) 

34.5 (5.0) 

.5 

13 

4.44 (1 

X 

10^) 

2.37 

41.37 (6000) 

20.68 (3000) 

34.5 (5.0) 

.5 

14 

6.67 (1, 

,5 

X 10®) 

2.37 

41.37 (6000) 

20.68 (3000) 

34.5 (5.0) 

.5 

15 

8.90 (2.0 

6 

X 10 ) 

2.37 

41 .37 (6000) 

20.68 (3000) 

34.5 (5.0) 

.5 

16 

2.67 (6 

X 

10®) 

2.37 

41.37 (6000) 

20.68 (3000) 

68.9 (10.0) 

.5 

17 

2.67 (6 

X 

10®) 

2.39 

41.37 (6000) 

20.68 (3000) 

101 (14.7) 

.5 

18 

2.67 (6 

X 

10®) 

1.0 

- 

20.68 (3000) 

34.5 (5.0) 

- 
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TABLE IV-II 


PARAMETRIC ANALYSIS RESULTS 


MODE 2 


. CASE 

THRUST 

SPLIT 

FI VAC 

ISPo 

VAC 

ISPo 

S.L. 

'^i 

tp 

^s 

^‘noz ^ 

HPF 

HPOX 

u 

Wjp 

WSF 

WSOX 

«TS 

’’overall 

F2 VAC 

ISPq 

^11 

n 

OVFHALL 

1 

59.3/40.7 

676440 

390.4 

346.3 

63.2 

34.13 

48.2 

82.33 

290.4 

813.2 

1103.6 

78.6 

550.3 

628.9 

.9669 

285533 

454 

107.9 

.9’ 37 

i 2 

59.?/40.7 

683403 

386.2 

339.1 

42.4 

34.53 

99.42 

133.95 

293.8 

822.7 

1116.5 

81.6 . 

.571.3 

652.9 

.9676 

288306 

441.6 

57.1 

.9743 


59.4/40.6 

673341 

391.9 

, 349.2 

75.9 

34.02 

31.72 

65.74 

289.5 

810.5 

1100 

77.3 

541.1 

618.4 

.9649 

284192 

459.6 

157.3 

.9725 

A 

59.4/40.6 

671416 

392.9 

351.2 

84.5 

33.94 

23.57 

57.51 

288.8 

808.7 

1097.5 

76.4 

534.8 

611.2 

.9638 

283209 

463.4 

206.1 

.9718 

1 s ' 

20.2/79.8 

682766 

422.4 


52.8 

11. 7T- 

94.9 

106.61 

99.67 

279.1 

378.77 

154.7 

108.3 

1237.7 


549583 

AAA.O 

59.4 


1 6 

72.9/27.1 

674274 

380.4 

338.5 

68.3 

41.82 

32.21 

74.03 

355.9 

996.4 

1352.3 

52.5 

367.5 

420 

.9653 

192676 

458.5 

156.9 

.9714 

7 

1 

79.2/20.8 

673266 

376.0 

335.1 

71.0 

45.37 

24.82 

70.19 

386.0 

1080.9 

1466.9 

40.5 

283.4 

323.9 

.9647 

149604 

; 461.9 

200. E 

.9709 

; 8 

j 

81.3/18.7 

672933 

374.5 

333.9 

72.0 

46.58 

22.31 

68.89 

396.3 

1109.7 

1506 

36.4 

254.7 

291.1 

.9644 

134593 

462.4 

222.4 

.9692 

I 9. 
1 

59.7/40.3 

702529 

370.5 

316.4 

26.2 

111.37 

154.84 

266.21 

319.9 

895.7 

1215.6 

85.1 

595.6 

680.70 

.9635 

296 336 

435.3 

45.0 

.5722 

, 10 

59.4/40.6 

684471 

383.8 

336.4 

A5.A 

52.53 

73.98 

126.31 

299.4 

838.3 

1137.7 

80.7 

565.2 

645.9 

.9658 

288 803 

447. J 

77.7 

.5726 

; 11 

59,2/40.8 

671761 

394.6 

352.4 

80.1 

25.20 

35.72 

60.92 , 

285.0 

797.9 

1082.9 

77.5 

542.3 

619.8 

.9667 

283 410 

457.3 

136.7 

.’’5725 

12 

59.3/40.7 

563712 

390.4 

346.2 

63.2 

28.45 

40.18 

68.63 

242.1 

677.8 

919.90 

65.5 

458.7 

524.2 

.9667 

237950 

453.9 

107.9 

.4735 

’3 

59.3/40.7 

1127339 

390.6 

346.5 

63.2 

56.87 

80.28 

137.15 

483.9 

1354.9 

1838.8 

130.9 

916.6 

1047.5 

.9673 

475857 

454.3 

107.9 

.''«2 

U 

59.3/40.7 

1690938 

390.7 

346.6 

63.2 

85.28 

120.27 

205.65 

725.7 

2031.8 

2757.5 

196.3 

1374.2 

1570.5 

.9677 

713750 

454.5 

107.9 

.9747 

15 

59.3/40.7 

2254520 

390.8 

346.7 

63.2 

113.69 

160.44 

274.13 

967.4 

2708.6 

3676 

261.7 

1831.6 

2093.3 

.9679 

951634 

454.6 

107.9 

. 5748 

16 

59.5/40.5 

643040 

380.4 

355.0 

36.5 

33.39 

46.81 

80.20 

284.1 

795.5 

1079.6 

76.3 

534.4 

610.7 

.9684 

271292 

444.2 

62.5 

.3739 

17 

59.6/40.4 

631733 

374.4 

355.6 

27.0 

33.38 

46.62 

80.00 

284.0 

795.1 

1079.1 

76.0 

532.2 

608.2 

.9692 

266506 

438.2 

46.3 

.5731 

‘ 18 

0/100 

686078 

441.0 

385.6 

49.1 

- 

119.25 

- 

- 

- 

- 

194.5 

1361.4 

1555.9 

.9736 

686078 

AAl 

49.1 

.9768 
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IV, A, Preliminary Geometric Analysis (cont.) 

Figures IV-1 through IV-3 show the dependence of engine perform- 
ance and area ratio on thrust split. Decreasing Mode I performance as thrust 
ratio increases is due to the increased contribution of the lower-performing 
primary flow to the overall performance. The Mode II delivered Isp increases 
with thrust ratio because the Mode II area ratio is increased. This increase 
in area ratio occurs because the secondary throat decreases with respect to 
the nozzle exit area as the thrust ratio increases. 


The Mode I area ratio increases with thrust split because of 
increases in primary throat area with thrust split. From geometry, the Mode I 
area ratio can be defined as: 





■^ts 


+ e 



The primary (Ep) and secondary (€$) streamtube area ratios are constant 
for a given ratio of specific heats (r) and pressure ratio (Pg/Pc)* 

Therefore, increasing Atp results in a larger ^i* 

Figures IV-4 and IV-5 show the effect of variation in thrust on 
specific impulse and throat area. Throat area will increase with increasing 
thrust, but the specific impulse will essentially be unchanged. 

Figures IV-6 through IV-8 show the effects of chamber pressure on 
throat area, area ratio, and specific impulse. Increasing chamber pressure or 
chamber pressure ratio results in decreasing throat areas and increasing area 
ratio, in turn resulting in an increase in specific impulse. 


Figures IV-9 and IV-10 show how decreasing nozzle exit pressure 
will result in increasing nozzle area ratio and specific impulse. 
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ENGINE ISP EFFICIENCY 
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Figure IV-8. Engine I „ versus Chamber Pressure 
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IV, A, Preliminary Geometric Analysis (cont.) 


The previously listed geometric parameters were used in deter- 
mining the dual-expander nozzle design. The primary nozzle contour was found 
by truncating a bell nozzle optimized for the primary flow field area ratio to 
some specified length or area ratio. An annular nozzle with a conical inter- 
nal expansion section was positioned such that the flow angle and static pres- 
sure at the exit of the internal expansion section matched the flow angle and 
static pressure at the primary nozzle exit. 

These steps are part of the procedure developed for designing the 
dual-expander nozzle. The nomenclature for this procedure is shown in Figure 
IV-11. Derivation of the equations is given in Appendix B. 

The design procedure is as follov/s; 

1. Calculate the overall engine area ratio, primary and secon- 
dary flov/s streamtube area ratios, throat areas, and equivalent nozzle lengths 
for a given set of operating conditions. 

2. Generate a Rao nozzle contour for the primary nozzle by using 
the calculated streamtube area ratio and nozzle length. Truncate this contour 
to the desired area ratio value. The primary nozzle exit radius (R£p), 

flow angle (&l)» 2-D pressure are now known. 

3. Once the secondary nozzle chamber pressure and the desired 
exit pressure are known, calculate the required conical section area ratio. A 
conical section is used because the 2-D to 1-D Mach number ratio will be close 
to one. The conical section exit area (Ag^) can now be calculated for the 
secondary flow. 
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Figure IV-11. Dual Expander Nomenclature 


IV, A, Preliminary Geometric Analysis (cont.) 


4. Pick a value for either the gap (<^Ec) radius 

(Rg;^) at the conical section exit and solve for the other value by using 
Equation 1. 

Ag = 2 tt 6r Rp (1 ) 

^c ‘‘c ‘"c 


5. Choose an inner lip thickness (Dl) and calculate the secon- 
dary flow inclination angle (Os) from Equation 2. 


Cos 0 


2 (Rp - Rp - D, Cos 0, ) 
c P 


( 2 ) 


*SE, 


6. The exit geometry for the conical section is now known. The 
conical half-angle can also be found by 


®c = ®s - ®L 


The next step is to find the proper secondary nozzle throat 
radius (Rts) 9^P (^ts) occur for the inclination angle, 0$* 
and secondary nozzle throat area (A^^) calculated in Step 1. 

7. The secondary nozzle throat radius (Rtg) is calculated 
from Equation 4. 



2 sec 8^ 


sin 9L 

sin 0^ cos 0^ 


Rt, 


2 sec 0 

s 



+ cos 0|^ 


RE^ sin 0. 
c L 

sin 0^ cos 0 
s c 



(4) 


Equation (4) is of the form 
2 

ax + bx + c = 0 
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IV, A, Preliminary Geometric Analysis (cont.) 


! = -b±[ 


b - 4ac 


2a 


sin 0. 

a = 2 sec M - sip 9 cos 0 


Rp sin 0, 

b = -2 sec 0„ / Rr + D, Cos 0, - ^c ^ 

^ ' '"p *- sin 0g cos 0^ 


8. Once R^^ is known, the throat gap (6^ ) can be found 
by using Equation 5: ^ 


2uS . R. 


(5) 


The length of the conical section (Xc) is found by the following equation: 


( 6 ) 


Sin 0 


9. The remainder of the secondary nozzle contour is approximated 
with a third-order parabola: 


R = A + BX + CX^ + DX^ 

The constants A, B, C, D are determined by specifying the coordinates and 
angle of the outer conical section wall and specifying the coordinates and 
angle at the overall nozzle exit. The exit radius is found from Step 1, and 
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IV, A, Prel i mi nary Geometric Analysis (cont.) 

the exit angle is assumed to be the same as calculated for the complete pri- 
mary nozzle contour. 

Equations (1) through (6) show that once the exit radius 
(REq) and gap (5^ ) g^e picked, there is only one set of values, 

®s> Rts» '^ts> ^c> satisfy these equations and result in 

a symmetrical conical nozzle section. Table IV-III lists different values of 
Re , , and D|_ and their influence on throat gap and location, 

V w 

throat orientation, and length of the conical section. For a constant lip 
thickness, increasing the exit radius results in decreasing the exit gap, the 
throat orientation angle, and the distance of the secondary throat from the 
central axis. The secondary throat gap and the length of the conical section 
both increase with increasing Increasing the lip thickness for a 

given Rg^ and 5^^ causes the secondary throat inclination angle to 
increase and moves the throat further from the central axis. Both the throat 
gap and conical section length get smaller. 

B. DUAL-EXPANDER AERODYNAMIC BLEED FLOW MODEL 
1. Model Development 

An aerodynamic model has been developed for Mode II operation 
and is presented in this section. This model differs from the dual-throat 
thruster Mode II aerodynamic model in several respects. First, the input 
geometry for the dual expander aerodynamic model is very simple. It is only 
necessary to specify the throat area, A^^, and lip angle, ^l* The plume 
model is also very simple. The flow at the exit of the annular nozzle is 
assumed to be one-dimensional in a cross section normal to the nozzle lip. 

The flow is assumed to follow a Prandtl -Meyer expansion around the nozzle Tip 
until an angle is reached such that the pressure of the base region is 
achieved. The plume boundary is taken to be a straight line inclined at this 
angle. 
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TABLE IV-III 


DUAL EXPANDER GEOMETRIC INFLUENCES 



•"c 

^E 

‘'c 




s 

X 

o 

J 


0.5 

10. 

3.2986 

81.2 

50.2 

8.9978 

0.8526 

1.014 

O 

00 

II 

u 

0.5 

10.5 

3.1416 

61.4 

30.4 

8.8027 

0.8715 

1.933 

PtRUNC 

0.5 

10.75 

3.0685 

49.3 

18.3 

8.2898 

0.9254 

3.245 

.C = 4.3 

0.5 

10.85 

3.0402 

43.6 

12.6 

7.695 

0.9969 

4.575 

Ap = 207.26 in.^ 

0.5 

11.0 

2.9988 

33.4 

2.4 

3.9776 

1.9286 

12.76 

C 

0.5 

10.75 

3.0685 

49.3 

18.3 

8.2898 

0.9254 

3.2451 

A* = 48.2 in. ^ 

‘'s 

0.75 

10.75 

3.0685 

59.1 

28.1 

8.9741 

0.8548 

2.0697 

R^p = 9.32 in. 

1.0 

10.75 

3.0685 

68.1 

37.1 

9.4266 

0.8138 

1.4263 


0.5 

11.0 

2.988 

33.4 

2.4 

3.9776 

1.9286 

12.76 

=31 = 

0.75 

11.0 

2.988 

46.2 

15.2 

8.2536 

0.9294 

3.8051 


1.0 

11.0 

2.988 

56.7 

25.7 

9.1272 

0.8405 

2.2407 




IV, B, Dual -Expander Aerodynamic Bleed Flow Model (cont.) 


The shear layer is superimposed on the plume boundary. With 
one exception, the assumptions made in formulating the shear layer model are 
the same as those used for the dual throat (described in Section III). The 
exception is that the shear layer is assumed to have axial symmetry. Assuming 
axial rather than planer symmetry significantly complicates the formulation; 
however, as discussed in Section III, this assumption is necessary for the 
dual expander because of the proximity of the shear layer to the flow axis. 

In Reference 18, Mueller has presented an analysis for an axisymmetric shear 
layer for a gas with constant properties. The analysis presented here allows 
for the addition of bleed flow into the base recirculation region. As was the 
case with the dual throat aerodynamic bleed flow model, the bleed flow gas can 
have properties different from those of the primary stream. The dual expander 
shear model is illustrated in Figures IV-12 and IV-13. Figure IV-12 shows the 
shear layer developing along the jet boundary until a position is reached 
where the shear layer is turned and transformed into a wake. The position at 
which the shear layer is turned is called the "neck." Turning the flow 
abruptly at the neck results in a recompression shock (shown in Figure IV-12). 
The distance from the exit of the interior nozzle to the neck, Z3» is shown 
in Figure IV-13. The streamline deflection angle at the neck, 83, and the 
trailing wake radius, r^, also shown in Figure IV-13. The calculational 
procedure is illustrated in Figure IV-14. As was the case for the dual throat 
aerodynamic bleed flow model, the problem requires solving three nonlinear 
simulateous equations for three unknowns. The three unknown quantities are 
Pb/Pop» 8r/Hop and Yp^* The three equations are (A), (B), and (C) 
in Figure IV-14. These equations are: 


7T 2„ 
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. 2 ^ 1 ) 


W^/Wp 


A„ cos 0Q 
P 3 
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- A (I-j (rij) - I-| (t1j|) + (J-| (rij) - J-| = 0 (A) 
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HEIGHT OF BASE REGION 



Figure IV-13. Dual Expander Aerodynamic Model Geometry 
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Pq for outer stream 

Mach NO. 

Recompression Factor^ N 
% Bleed. 
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Estimates 
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METHOD: 

Newton's Method, given above estimates find 
converged values that satisfy condition of 

(A) CONTINUITY FOR TOTAL FLOW 

(B) ENERGY CONSERVATION 
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OUTPUT: 

Pb/Po . H/H, . Y 

CONSTANT PRESSURE BOUNDARY SHAPE 
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Figure IV-14. Calculational Procedure for the Dual Expander Aerodynamic Model 
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IV, B, Dual -Expander Aerodynamic Bleed Flow Model (cont.) 
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IV, B, Dual -Expander Aerodynamic Bleed Flow Model (cont.) 


An input manual for the dual expander aerodynamic bleed flow 
model is presented in Appendix C. 

2. Model Calibration 


The dual expander base flow model contains the following 
three parameters that must be determined from experimental data: 

a = The jet spreading parameter 

N = The Nash recompression factor 

k = A proportionality factor relating the radius of the 
trailing wake to the shear layer velocity thickness 

The approach taken in selecting each of these parameters is 
discussed below. Note that all of the data used is for air (y = 1.4) without 
bleed flow. 


a, the Jet Spreading Parameter 

The classical jet spreading expression measured by Korst 
(Ref. 6) for air was used, i.e.. 


a = 12 + 2.758Mjj 

For incompressible flow (M^ = 0), this expression gives the well-known 
result that 


Many alternate expressions for a are available in the literature (see Ref. 7, 
11, 18). 
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IV, B, Dual -Expander Aerodynamic Bleed Flow Model (cont.) 


N, the Nash Recompression Factor 

The data presented by Mueller (Ref. 18) was used to determine 
N. This data is for the flow of air past a cylinder. Measured values of base 
pressure ratio vs the Mach number of the approach flow are presented in Figure 
IV-15.* The theory predicts that base pressure depends only weakly on the 
trailing wake radius ratio. Measured values of the trailing wake radius ratio 
are shown in Figure IV-16. 

The computer model was calibrated by inputting the flow 
conditions that correspond to the experimental conditions and by inputting the 
trailing wake radius from Figure IV-16. It was found that 

N = .4 

gave the best fit to the data. This is consistent with the treatment used in 
the dual throat aerodynamic model computer program. It should be noted that 
Mueller (Ref. 18) used a value of N = 1 (i.e., the full value from the shock) 
and a different expression for a a (i.e., an expression that gives values that 
are greater than those given by the Korst expression). 

K, Shear Layer/Wake Proportionality Factor 

For the case of air without bleed flow, the trailing wake 
radius ratio can be obtained from Figure IV-16. In order to obtain a similar 
result for situations with bleed flow and with gases other than air, the 
following assumption was made. 


*Figures IV-15 and IV-16 are from Reference 18. 
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Figure IV-16. Trailing Wake Radius versus 
Mach Number for 6 = 1.4 




IV, B, Dual -Expander Aerodynamic Bleed Flow Model (cont.) 


It was assumed that the radius of the trailing wake is 
proportional to the velocity thickness of the shear layer at the neck, 5. 
Since this value is proportional to 

{x/r^)/a 

where x/r[j is the nondimensional shear layer length, it follows that 

Vb ' 

where 

K is the proportionality constant. 

It was found that 


K = 4 

i.e., this approximate value v/as consistently calculated when 

a = 12 + 2.758Mb 
N = .4 

and the cases shown in Figure IV-16 were calculated. 

C. DUAL EXPANDER FLOW FIELD ANALYSIS 

Performance calculations were carried out for both Mode I and Mode 
II operation. In order to carry out these calculations, it was convenient to 
divide the dual expander flow field into regions that could be analyzed 
separately. These regions, numbered I through V, are shown in Figures IV-17a 
and IV-17b. 
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a. DUAL EXPANDER MODE I OPERATION 
BOTH ENGINES ON 



b. DUAL EXPANDER MODE II OPERATION 
OUTER ENGINE ON, INNER ENGINE OFF 


Figure IV-17, Dual Expander Modes I and II Operation 
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IV, C, Dual Expander Flow Field Analysis (cont.) 


Region I is the flow field inside the interior thruster during 
Mode I operation. The downstream boundary of Region I is the exit plane of 
the primary thruster. Region I can be conveniently, accurately, and inexpen- 
sively calculated by using the TDK computer program (Ref. 3). 

Region II is the flow field inside the annular thruster. This 
flow field is the same during both Mode I and Mode II operation. The down- 
stream boundary of Region II is the exit plane of the annular thruster. 

Region II can be calculated by using the VNAP computer program (Ref. 4) that 
was discussed in Section III. VNAP is particularly well suited to calculating 
the transonic flow field in an annular nozzle. The VNAP calculations are 
somewhat lengthy and require dividing Region II into subregions in order to 
obtain sufficient accuracy. 

Region III is the dual expander flow field during Mode I operation 
(both engines on) that is downstream from Regions I and II. In order to mini- 
mize shock losses, the engine is designed to operate so that the static pres- 
sures in Regions I and II are equal at the point where these two flow streams 
come together. This point is the origin of a contact discontinuity, or slip 
line, between the two streams. The TDK program has been modified to calculate 
the flow field in Region III. In order to use the option, it is necessary to 
input flow conditions along the exit planes of Regions I and II. The two flow 
streams that comprise Region III are allowed to have different values of stag- 
nation pressure, specific heat ratio, and molecular weight. 

Region IV is the inviscid portion of the flow field that occurs 
downstream of Region II during Mode II operation (inner thruster off, outer 
thruster on). The interior boundary of Region IV is determined by the shape 
of the viscous-inviscid boundary of the base recirculation flow and wake. 
Region IV can contain shock waves. The VNAP computer program (Ref. 4) has 
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IV, C, Dual Expander Flow Field Analysis (cont.) 


been used to calculate the flow field in this region so that the effect of 
shock waves can be included. In order to maintain accuracy, it is necessary 
to divide Region IV into subregions. 

Region V is the base flow recirculation and wake region that is 
calculated by the dual expander aerodynamic model described in the previous 
section. 


A prototype design for the dual-expander thruster was generated by 
using the procedure defined in Section IV, A, and performance calculations were 
carried out by analyzing the five flow regions described above. These calcu- 
lations are described in the following subsections. 

1. Calculations for Mode II Operation 

Calculating the dual expander performance for Mode II opera- 
tion requires calculating Regions V, II, and IV that were shown in Figure 
IV-17b. Results obtained for each region are discussed below. The purpose of 
these calculations was to determine the influence of primary nozzle truncation 
and/or the introduction of bleed flow. In order to isolate these influences, 
it was necessary to keep the same secondary nozzle contour and vary the pri- 
mary nozzle exit area only. 

Region V 


Region V was calculated for the prototype engine for three 
bleed flow values and for three truncation positions for the dividing wall. 
Results from these nine cases are presented in Table IV-IV, labeled as Cases 2 
through 10.* The bleed flow values used were zero bleed, "moderate" bleed, 
and "maximum" bleed. 

*Case 1 was used as a sample case for VNAP for program checkout prior to com- 
pleting Cases 2 through 10. 
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TABLE IV-IV 


BASE FLOW STUDY, SUMMARY OF RESULTS 


Case 

No. 

cm 

«B 

in. 


Z3/Rb 

1 

24.97 

(9.83) 

21.° 

1.3 

2 

25.18 

(9.913) 

14.64° 

2.296 

3 

25.18 

(9.913) 

8.62° 

3.300 

4 

N/A 

N/A 

N/A 

N/A 

5 

22.47 

(8.8446) 

12.80° 

2.203 

6 

22.47 

(8.8446) 

7.02° 

3.249 

7 

22.47 

(8.8446) 

4.92° 

3.487 

8 

19.93 

(7.8474) 

10.38° 

2.730 

9 

19.93 

(7.8474) 

6.17° 

2.775 

10 

19.93 

(7.8474) 

3.53° 

3.240 


Vb 

% 

Bleed 

"p 


Truncation 

.5 

N/A 

N/A 

N/A 

Minimum 

.4 

0% 

2.79 

23.4° 

Minimum 

.5 

.05% 

2.79 

19.2° 

Minimum 

N/A 

N/A 

2.79 

N/A 

Minimum 

.5 

0% 

2.19 

25.2° 

Moderate 

.6 

.25% 

2.19 

21.7° 

Moderate 

.7 

.5% 

2.19 

20.7° 

Moderate 

.5 

0% 

1.016 

34.2° 

Maximum 

.7 

1% 

1.016 

32.4° 

Maximum 

.8 

3% 

1.016 

35.3° 

Maximum 



IV, C, Dual Expander Flow Field Analysis (cont.) 


The aerodynamic model is a small bleed flow model. If the 
bleed flowrate becomes too large, it will lead to the following: the base 

pressure region will not close; there will be no recirculation, the model will 
no longer be applicable. "Maximum" bleed was defined as an amount just short 
of causing this condition. "Moderate" bleed flow was taken as half of the 
"maximum" bleed flow value. 

"Minimum" truncation corresponds to arbitrarily choosing to 
truncate the primary nozzle to a geometric area ratio of 8:1. "Maximum" wall 
truncation was obtained by truncating the dividing wall just downstream of the 
Mach one position on the annular nozzle side of the wall. The "moderate" wall 
truncation position was taken as halfway between "maximum" and "minimum" trun- 
cation. 


The parameters tabulated in Table IV-IV were illustrated in 
Figure IV-13. Inputs to the dual expander aerodynamic model that differ in 
each case are 1) the exit radius of the interior nozzle, r^,, 2) the Mach 
number of the approach flow at the end and upper side of the dividing wall, 
Mp, and 3) % bleed for inner engine. The approach flow angle is 31°. The 
operating conditions for these cases are listed in Table IV-V. 


TABLE IV-V 


Flow 

Mixture Ratio 
Y 
Cp 

Molecular Weight 
Total Temperature 
Total Pressure 

Throat Area 


MODE II OPERATING CONDITIONS 


Outer Nozzle 
LOX/LH 2 
7.0 
1.18 

.81371 cal/gm °K 
16. 

3778°K (680p°R) 

20.68 MPa 
(3000 psia) 

281 cm2 (43.5 in. 2) 


Bleed to 
Inner Nozzle 

GG 

.72 

1.389 

2.04454 cal/gm °K 
3.467 

712°K (1281°R) 

N/A 

N/A 
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IV, C, Dual Expander Flow Field Analysis (cont.) 


The geometry of the shear layer and wake as calculated by the 
dual expander aerodynamic model is defined by the parameters 0^ 

Z 3 /rb, and r^^/rb. Calculated values for these parameters are listed 
in Table IV-IV. Also listed is the value of the shock angle, corres- 
ponding to a deflection angle of Q^. 

The approach Mach number at the three truncation positions 
was input from the VNAP results for Region II. The shear layer and wake shape 
calculated for Region V was used to define the lower boundary for Region IV. 
The Region II and IV calculations are discussed below. 

Region II 

Figure IV-18 shows an enlargement of Regions I and II. The 
upper wall of Region II is steep, having a maximum wall angle of 67°. To 
obtain accuracy. Region II was divided into an upstream subregion and a down- 
stream subregion, and these were calculated separately by using VNAP. 

Mach number contours calculated for the upstream subregion 
are shown in Figure IV-19. Inlet conditions for this subregion were pre- 
scribed to be one-dimensional parallel flow along a vertical cross section at 
the extreme left of Figure IV-19. The calculations show that, along the lower 
wall, the flow is turned inward abruptly while subsonic (M = .4). Along the 
upper wall, the nozzle throat, a circular arc with a small radius of curva- 
ture, causes the flow to be choked. The mass flowrate through the nozzle is 
found to be only about 91% of the one-dimensional flow value. At the exit of 
this subregion, the flow is found to be supersonic, and, therefore, indepen- 
dent of the conditions downstream. 
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IV, C, Dual Expander Flow Field Analysis (cont.) 


Mach number contours calculated for the downstream subregion 
of Region II are shown in Figure IV-20. Inlet conditions for this subregion 
are the exit conditions of the upstream subregion.* An interesting feature of 
this plot is that the flow near the lower wall is found to expand like a one- 
dimensional flow. This tends to validate the design approach of using a coni- 
cal section in the annular portion of the secondary nozzle so that the pres- 
sure may be determined by using one-dimensional calculations. As was dis- 
cussed previously, it is important to be able to estimate this pressure value 
when designing the nozzle. It can also be seen from Figure IV-20 that the 
flow is being compressed by the upper wall. Downstream, this compression 
develops into a shock. 

The Region II calculations are used to supply inlet condi- 
tions to Region IV and also to the Region III, Mode I, calculations. 

Region IV 

Inlet line conditions for Region IV were supplied by selecing 
the appropriate flow cross section from the VNAP results for Region II. The 
upper wall boundary was determined by using the design approach described in 
Section IV, A. The lower boundary was defined by the parameters 03» 

Z3/>"b. and rw/rb* 

Region IV was divided into an upstream subregion and a down- 
stream subregion. The following 5 plots were prepared for each case: 

*Each region is actually overlapped to increase the accuracy. The last line 
or two calculated by VNAP is not accurate because of the outflow boundary 
condition that must be imposed by the method to close the problem. (See 
Ref. 4.) Thus, the VNAP calculations have been overlapped and the last 
two lines were discarded. 
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IV, C, Dual Expander Flow Field Analysis (cont.) 


(1) velocity field, (2) pressure, (3) density, (4) temperature, and (5) Mach 
number contours. Only the Mach number contours are discussed in this report. 
Mach number contours for these subregions are presented in Figures IV-21 and 
IV-22 for case 2, i.e., zero bleed and no truncation. Thus, Figures IV-19, 
-20, -21, and -22 show the entire Mode II flow field predicted for the base- 
line nozzle contour and zero bleed flow. 

The lower boundary in Figure IV-21 is the shear layer. The 
lower left-hand corner in the figure shows the Prandtl-Meyer expansion fan at 
the lower lip of the annular nozzle where the flow turns through 45.64°. The 
finite-difference method used by VNAP cannot accurately calculate this expan- 
sion fan because the Mach number gradients are too great for the relatively 
coarse mesh spacing required by VNAP. Thus, this expansion fan is a source of 
error in the VNAP calculations. 

The lower boundary of Figure IV-22 is the shear layer and 
wake. A compression shock is induced at the corner where the flow is 
deflected through an angle of 14.64°. This deflection angle corresponds to a 
shock angle of 23.4°. Figure IV-22 clearly shov/s this "neck" shock; however, 
the shock angle appears to be less than 23.4°. VNAP predicts that the neck 
shock will exit the nozzle far from the upper wall. 

A shock can also be seen in the upper half of Figure IV-22. 
This shock is induced by compression of the flow by the upper wall. 

Results 


Results of the Mode II performance predictions are summarized 

in Table IV-VI. 
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TABLE IV-VI 


SUMMARY OF RESULTS, MODE II PERFORMANCE PREDICTION 


l“> 

l/> 
f n 

BASE 

TRUNCATION PRESSURE 

MR?. (psia) 

% BLEED 

Wg/Wj X 10^ 

n 

INITIAL LINE 
Kq/sec (Ibm/sec) 

F 

INITIAL LINE 
KN (K IbF) 

EXIT LINE 
Kq/sec (Ib/sec) 

>^1 

EXIT LINE 
MN (K IbF) 

^2 

EXIT 

MN 

LINE 
(K IbF) 

* 

^SP 

1 

sec 


Min. 

.0183 

(2.65) 

0 

262.7 

(578.7) 

704.7 

(158.4) 

260.0 

(573.2) 

1.089 (244.8) 

1.107 

(248.9) 

423.1 

3 

,Xin. 

.0322 

(4.67) 

.05 

262.7 

(578.7) 

704.7 

(158.4) 

270.8 

(597.1) 

1.154 (259.6) 

T.151 

(258.8) 

443.5 

• 

Min. 

.0534 

(3.92) 

^ .07 

262.7 

(578.7) 

704.7 

(158.4) 

260.7 

(574.7) 

1.132 (254.5) 

1.111 

(249.4) 

439.8 

5 

Mod. 

.123 

(17.89) 

0 

263.1 

(579.7) 

665.4 

(149.6) 

282.5 

(622.9) 

1.212 (272.4) 

1.198 

(269.3) 

469.9 

c 

*‘od. 

.131 

(27.7) 

.25 

263.1 

(579.7) 

665.4 

(149.6) 

279.6 

(616.5) 

1.215 (273.2) 

1.191 

(267.7) 

471.2 

7 

"cd. 

.222 

(32.2) 

.5 

263.1 

(579.7) 

665.4 

(149.6) 

282.0 

(621.7) 

1.222 (274.6) 

1.198 

(269.3) 

473.7 

3 

Max. 

1.727 

(250.5) 

0 

259.5 

(571.7) 

672.9 

(151.3) 

382.2 

(842.7) 

1.675 (376.6) 

1.631 

(366.6) 

558.7 

9 

Max. 

2.14 

(31C.4) 

1 

259.5 

(571.7) 

672.9 

(151;3) 

374.3 

(825.2) 

1.678 (377.2) 

1.598 

(359.3) 

559.8 

17 

Max. 

2.437 

(353.5) 

3 

259.5 

(571.7) 

672.9 

(151.3) 

377.3 

(831.8) 

1.678 (377.3) 

1.604 

(360.5) 

659.5 


* (Isp) , = Pi/?nit1a1 

line 

“ (Up), = F,/m 

^ exit 
line 



2 

see 


433. e 

434. ’ 

422.3 

434.1 

433.2 

435. :■ 

435.4 

432.5 



IV, C, Dual Expander Flow Field Analysis (cont.) 


For cases 2, 3, and 4, where the wall dividing the upper and 
lower nozzles were as shown in Figure IV-18, the approach Mach number was 
found to be 2.79. With zero bleed flow, the predicted Prandtl -Meyer turning 
angle at the nozzle lip was 31 + 14.64 = 45.64®. This turn gives a Mach num- 
ber of 4.62 for the shear layer boundary. For a Mach number this high, very 
little bleed flow is required to open the wake. The aerodynamic model pre- 
dicted a trailing wake radius ratio approaching one for a bleed flow of only 
.07%. The values of base pressure predicted for cases 2, 3, and 4 are shown 
in column one of Table IV-VI, and are found to be quite low. For cases 5, 6, 
and 7, which have moderate v/all truncation, and cases 8, 9, and 10, which have 
maximum wall truncation, the approach Mach numbers are 2.19 and 1.016, respec- 
tively. Consequently, the predicted values of base pressure and bleed flow 
required to open the wake are larger than for the zero truncation cases. 

The remainder of Table IV-VI summarizes the value of mass 
flow, thrust, and specific impulse that were predicted by VNAP. Mass flowrate 
and thrust were calculated by using two different methods: 

1) Initial Line Method 

'^initial ~ Integral across the VNAP start line of 
line Figure IV-21. 

Fi = Integral of axial momentum + pressure 

across the VNAP start line of Figure 
IV-21, + integral of axial pressure along 
both the upper and lower boundaries 
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IV, C, Dual Expander Flovi/ Field Analysis (cent.) 


2) Exit Line Method 


= Integral across the VNAP line at the 
end of the nozzle 

= Integral of axial momentum + pressure 
across the VNAP line at the end of the 
nozzle 


Two different values of specific impulse were calculated by using the above 
quantities: 


'^'exit 

line 

f^2 


dsp)l = 

line 

dSP)2 = F2/mexit 
line 

Referring to Table IV-VI, it can be seen that there is substantial error 
growth in both the mass flow and thrust calculations, especially for cases 8, 
9, and 10. This is primarily because of the inability of VNAP to properly 
calculate the large gradient region in the vicinity of the Prandtl-Meyer 
expansion fan. These gradients are especially large in the last three cases 
because of the low Mach number (1.016) of the flow upstream of the fan. 

Nearly all of the observed error growth occurs in this region, i.e., in cal- 
culating the first several VNAP solution lines. The Initial Line Method 
cannot be used to predict the specific impulse, because a thrust calculation 
with error growth should not be ratioed to a mass flow value with no error 
growth. The Exit Line Method, however, seems to give a consistent prediction 
of specific impulse for all 9 cases. This is because the mass flow and thrust 
integrals are very similar and experience nearly the same error growth. All 
nine cases are found to yield a specific impulse of 434 + 1.5 seconds. Cor- 
recting these values by using the VNAP correction factor (.978) obtained 
during the dual throat analysis, the Mode II divergence efficiency was found 
to be 100% + .2%. 
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IV, C, Dual Expander Flow Field Analysis (cont.) 

Although the results for Mode II operation obtained with 
VNAP appear to be qualitatively correct, it is clear that a more accurate 
method must be used to predict the nozzle performance and that a cold-flow 
test program is needed to verify these results. 


Performance Calculations For Mode I Operation 


Calculating the dual expander performance for Mode I opera- 
tion required calculating Regions I, II, and III shown in Figure IV-17a. 
Results obtained for each region are discussed below. 



Results of the Region II calculations were available from the 
Mode II analysis. Region II results were discussed in the preceding subsec- 
tion and are applicable to both modes of operation. Next, Region I was cal- 
culated by using the TDK computer program, with operating conditions for the 
inner nozzle as shown below in Table IV-VII. 


Flow 

Mixture Ratio 
Y 

Cp 

Molecular Weight 
Total Temperature 
Total Pressure 
Throat Area 


TABLE IV-VII 

MODE I OPERATING CONDITIONS 
Outer Nozzle Inner Nozzle 

LOX/LH 2 LOX/RP-1 

7.0 2.8 

1.18 1.199 

.8137 cal/gm°K .44985 cal/gm°K 

16. 26.6 

3778°K (6800°R) 3972°K (7150°R) 

20.68 MPa (3000 psia) 41.37 MPa (6000 psia) 

280.6 cm2 (43.5 in. 2) 220.6 cm2 (34.2 in. 2) 
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IV, C, Dual Expander Flow Field Analysis (cont.) 


It was intended to adjust the chamber pressure of the inner nozzle so that the 
calculated exit pressure would match the VNAP exit pressure. However, a 
satisfactory match was obtained on the first attempt. Pressures calculated 
along the dividing wall for Regions I and II are shown in Figure IV-23. Pres- 
sures at the downstream end of the wall are nearly equal, about 105 psia on 
both sides. It is important to match the static pressures at the end of the 
wall to avoid forming a strong shock wave at the position where the two flow 
streams come together. 

Figure IV-23 demonstrates that a large, negative pressure 
gradient occurs along the upper wall and that a much lower gradient occurs 
along the lower wall. Thus, if the wall is truncated, a large difference in 
static pressure will result at the end of the wall. For example, at the 
"moderate" truncation position, the upper and lower wall pressures are found 
to be 1970 (286) and 979 (142) KPa (psia), respectively. For "maximum" trun- 
cation, the values are 11606 (1683) and 1248 (181) KPa (psia). In order to 
match the pressures at the lip for these truncations and still operate at the 
desired chamber pressure requires a complete redesign of the nozzle. 

Region III 

Region III, shown in Figure IV-17a, is the dual expander flow 
field during Mode I operation (both engines on) that is downstream from 
Regions I and II. The TDK program has been modified to calculate Region III. 
Previously, TDK would calculate the flow field only if the difference between 
the streamtube and stagnation pressure was within 10% of each other. For the 
dual expander this difference is 20%. 

Region III is divided into an inner region and an outer 
region. The inner region is composed of exhaust gases from the LOX/RP-1 
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Figure IV-23. - Inner and Outer Wall Pressures (Predicted) 


IV, C, Dual Expander Flow Field Analysis (cont.) 


propellants at a mixture ratio of 2.8. The outer region is composed of 
exhaust gases from the LOX/LH 2 propellants at a mixture ratio of 7.0. 

Operating conditions for the two engines are as shown in Table IV-VII. 

Inputs along the start line for the inner region were input 
from the exit plant of the TDK calculation for Region I. Inputs along the 
start line for the outer region were input from the last vertical line of the 
VNAP calculation of Region II. Values input are pressure, density, flow 
angle, and velocity at each point. The two regions are divided by a slip 
line, i.e., a double valued streamline with equal values of pressure and flow 
angle, but different values of velocity, density, Mach number, and temperature 
on each side. 


The TDK calculations were made by removing the thickness of 
the dividing wall and then correcting the TDK calculations for mass flow and 
thrust. This calculation showed the divergence efficiency to be .995, 
matching the simplified prediction. 

D. RECOMMENDATIONS 

1. Mode II Performance Prediction Program 

After careful analysis of the existing Mode II performance 
prediction methods for the flow field of the dual expander propulsion concept, 
it is now felt that a method of characteristics analysis and computer program 
capable of more precise prediction techniques is needed. A sketch illus- 
trating the geometry and the major features of the flow field is presented in 
Figure IV-24. In view of the complexity of this flow field, a finite- 
difference method that smears Prandtl-Meyer expansions and oblique shock 
waves, such as the VNAP program, is considered inadequate to fully determine 
the potential performance of the dual expander concept. 
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Figure IV-24. Dual Expander Geometry and Flow Field Model 
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The general features of the proposed flow model are pre- 
sented, along with a discussion of the overall numerical algorithm. The pro- 
gram would be less complex, faster and cheaper to run than VNAP. Such advan- 
tages allow the program to be run parametrically to determine an optimum Mode 
II design. The possibility also exists to expand the program to analyze dual 
throat Mode II operation. 

2. Flow Model 


The flow model should be developed for steady axi symmetric 
flow. The flow should be allowed to be rotational to account for the entropy 
gradient created by the oblique shock wave emanating from the inner wall. The 
gas would be assumed to have constant specific heats and molecular weight. 

The supersonic initial value line (IVL) would be input to the analysis in 
tabular form. Both the outer wall and inner wall would be input in tabular 
form. A Prandtl -Meyer expansion would emanate from point L on the inner wall. 
A discrete left-turning obliquie shock wave would emanate from point N on the 
inner wall and be reflected from the outer wall if it intersects that wall. 

All unit processes for calculating the flow field would be 
based on the method of characteristics. The interior point unit process would 
be constructed by following the two Mach lines. Outer and inner wall unit 
processes would follow the appropriate Mach line and the streamlinb (i.e., the 
wall itself). 


3. Overall Numerical Algorithm 

The overall logic would follow left-running Mach lines. The 
first left-running line would emanate from the second point on the initial- 
valve line (IVL) below the outer wall, point 12. That Mach line would be 
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extended through the flow field until it intersects the outer wall. Subse- 
quent left-running Mach lines would emanate from the remaining points on the 
IVL, in decreasing order, until the Mach line from point L has been extended 
to the outer wall. Left-running Mach lines would then be generated through 
the Prandtl-Meyer expansion fan at point L until the required amount of flow 
turning at point L has been achieved. Each of these left-running Mach lines 
would be extended across the flow field to the outer wall. A right-running 
Mach line would then be extended from the second point on the last left- 
running Mach line of the expansion fan, point A, to the inner wall. A left- 
running Mach line would be extended from that wall point across the flow field 
until it intersects the outer wall. This procedure would be repeated until 
the left-running Mach line emanating from point N on the inner wall has been 
extended across the flow field. A left-running oblique shock wave would then 
be generated at point N and extended across the flow field. This shock wave 
would either intersect the outer wall or pass beyond the nozzle exit. If it 
intersects the outer wall, a shock reflection unit process would be employed 
to determine the flow properties across the reflected shock wave. Left- 
running Mach lines would be extended from points on the inner wall downstream 
of point N until a left-running Mach line passes beyond point K on the outer 
wall (the exit lip point). If the inner wall is a constant radius section, 
the solution is complete. If not, the solution would be continued until the 
flow field is determined on the inner wall at point W. The nozzle thrust 
would then be determined by calculating the momentum flux and pressure thrust 
across the tabular initial-value line, line IL, and adding to it the integral 
of the pressure forces along the outer wall. 
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V. 


PRELIMINARY HEAT TRANSFER ANALYSIS 


The objective of Task IV was to conduct a preliminary thermal analysis 
for both the dual throat and dual expander configurations. The dual throat 
studies represent a continuation and refinement of the work initiated in Ref. 
2. Vehicle studies following that effort resulted in a change in the base- 
line design point: primary and secondary chamber pressures, respectively, 

were increased from 20.68/14.48 MPa (3000/2100 psia), to 27.58/19.31 MPa 
(4000/2800 psia), and the Mode I secondary (LOX/RP-1) thrust was increased 
from 60 to 70 percent of the total. Both changes make chamber cooling more 
difficult and result in a baseline design point not explicitly considered in 
the parametric studies of Ref. 2. A new baseline design was developed which 
included special emphasis on the lip region formed by the termination of the 
primary chamber and on the effect of primary chamber plume impingement in the 
secondary throat region during Mode II operation. New parametric design 
studies investigated the effects of independent changes in chamber pressure 
and Mode I thrust split. 

Previous thermal analysis of the dual expander concept was limited to a 
very preliminary IR&D study (Ref. 21) which indicated problems in apportioning 
the limited hydrogen coolant flow among the three heated surfaces. Three can- 
didate baseline designs were investigated in the present study: (1) regenera- 

tive cooling with hydrogen only; (2) regenerative cooling with hydrogen and 
oxygen; (3) a transpiration-regenerative cooling combination with hydrogen. A 
special study of the supersonic lip region formed by the termination of the 
primary chamber was also included in the baseline design effort, conducted for 
a primary chamber pressure of 41.37 MPa (6000 psia) (LOX/RP-1), secondary 
chamber pressure of 20.68 MPa (3000 psia) (L0X/H2)> and a Mode I primary 
thrust equal to 59 percent of the 698 MN (600K IbF) total. Parametric studies 
investigated the effects of changes in thrust, thrust split, and both chamber 
pressures for the first cooling concept noted above. 
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A. SPECIAL STUDIES 
1. Lip Regions 

One of the heat transfer problem areas unique to the dual 
throat and dual expander concepts is the lip formed by the termination of the 
primary chamber. During Mode I operation, the primary and secondary chamber 
flows create a wake region behind this lip. Therefore, a special study of the 
lip regions was conducted to identify the boundary conditions established by 
the wake flow and the design problems associated with regenerative cooling of 
a corner heated on both sides. 

Wake region heat transfer data in the literature are primar- 
ily for supersonic flow over axisymmetric bodies, with the sidewall boundary 
layer thickness small in relationship to the base diameter. The concept of 
interest to us incorporates what is, essentially, a plane surface with a rela- 
tively thick boundary layer and subsonic flows for the dual throat and super- 
sonic flows for the dual expander. Because of these differences, a conserva- 
tive design approach has been adopted, using the literature data as a guide. 
The latter indicate that heat fluxes at the center of the wake region are 
typically 80-90 percent of the sidewall heat flux immediately upstream of the 
wake region (see Refs. 22 and 23), but that a significant reduction in heat 
flux is observed as it moves away from the center of the wake toward the 
sidewall (Refs. 22-24). Figure V-1 shows the radial heat flux variation 
observed in Ref. 22; Ref. 24 indicates a greater reduction, with the local 
base heat flux being proportional to (1 - r/RB)l*2 for r/Rg ^ 0.8. 

Preliminary mechanical design studies indicate that a conven- 
tional rectangular channel construction, as illustrated in Figures V-2 and V-3 
for the dual throat and dual expander, respectively, can be utilized in the 
lip regions. Therefore, the emphasis in the present study was on the lip 
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region corners, vvhere a severe cooling problem occurs since, at a land center- 
line, heat input from two surfaces must be conducted in three directions to 
reach cooled surfaces. A three-dimensional SINDA model of such a corner was 
developed and is shown in Figure V-4. The model does not include any part of 
the electroformed nickel closure and thus is limited to that part of the 
corner adjacent to the coolant channels. Since the electroformed nickel can 
be relatively thin in the lip region, the corner model covers a significant 
fraction of the lip half-thickness. Therefore, the procedure used in the pre- 
sent study to provide a conservative design relative to the boundary condi- 
tions of Figure V-1 was as follows; 

(1) Apply a uniform heat transfer coefficient over the end 
wall equal to 80 percent of the sidewall coefficient. 

(2) Set the lip channel depth equal to the sidewall channel 
depth (symmetrical corner model), thereby accommodating 
an even higher heat flux at the center of the lip. 

It should be noted that the sidewall boundary layer was assumed to be 
unaffected by the throat region shock structure. In addition, the coolant 
heat transfer coefficients were considered to be unaffected by the turn at the 
lip corner; the actual enhancement due to curvature was omitted to compensate 
for the lack of a channel corner radius and increased wall thickness in the 
wall conduction model. 

The corner model was used to define the lip region channel 
depth required for each baseline design to limit the corner temperature to 
approximately 756°K (900°F). Figure V-5 shows a typical gas-side wall temper- 
ature distribution over a land centerline in the vicinity of a lip corner for 
the dual throat concept; temperatures over the channel centerline are only 
slightly lower. 
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Figure V-4. Lip Region 3-D Thermal Model 
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In order to provide lip region channel depths for the para- 
metric design studies, we used a simple one-dimensional conduction model with 
an effective wall conductance and coolant surface area defined from results of 
the three-dimensional SINDA model. The validity of this model was demon- 
strated by successfully predicting a dual throat baseline design, using a dual 
expander design to define the one-dimensional model parameters. In all cases, 
the total pressure drop for the two lip region corners was taken as three 
velocity heads. Results from the 3D baseline and ID parametric models are 
presented in subsequent sections as part of the overall results for each 
concept. 


2. Secondary Throat Analysis for Dual Throat Mode II 

During normal Mode II operation of the dual throat concept, 
the secondary chamber bleed flow causes the primary chamber plume to attach to 
the secondary wall downstream of the throat such that the plume flow is nearly 
parallel to the wall. In this case, throat heat fluxes are expected to be 
well below those during Mode I operation. Of interest herein is the possi- 
bility of bleed flow failure and impingement of the primary plume upstream of 
the secondary throat, with resultant heat fluxes that might be more severe 
than Mode I values. 

A number of studies of rocket exhaust plume impingement on 
adjacent vehicle structures were found in the literature (Refs. 25-27). 
Analyses and correlation of these data split the problem into strong and weak 
shock regions, with the former treated as a stagnation point region as far as 
heat transfer is concerned. However, because of the effect on the plume 
boundary of the back pressure in the secondary chamber, formation of a strong 
shock region was not considered likely in the present problem. Therefore, the 
VNAP program was used to define the inviscid flow field for zero bleed flow. 
This program provides an exact inviscid solution, and its use for the nominal 
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bleed flow case was already planned as part of Task I. A boundary layer 
analysis based on the edge conditions determined by VNAP was used to define 
the heat transfer to the secondary wall. 

Figure V-6 shows the Mode II pressure and mass velocity pro- 
files at the secondary wall from VNAP for the baseline design with zero bleed 
flow. The maximum pressure and mass velocity due to the plume impingement 
occurs 4.3 cm (1.7 in.) upstream of the throat. The freestream velocity, 
shown in Figure V-7, is essentially constant in the region of interest, i.e., 
from the attachment point through the zone of maximum heating. This feature of 
the inviscid flow was used to advantage herein, since it greatly simplifies 
the integral momentum equation for the boundary layer. For constant velocity, 
this equation becomes 

^ + 1 0 = C ./2 

dx Pg 'dx 

This form of the momentum equation was solved by using the Van Driest skin 
friction law for compressible flow based on the Von Karman turbulent mixing 
length (Ref. 28): 

Re, - 0.19, V exp (o.A . ^ 

in which ^ is a function of Mach number and wall temperature ratio. The Mach 
number and freestream static temperature distributions from VNAP are nearly 
uniform. A design wall temperature of 811°K (1460°R), consistent with the 
Mode I design studies described later, was used. 

The momentum thickness at the attachment point must be spe- 
cified to start the above boundary layer analysis. Two cases were considered. 
In one case, the momentum thickness of the plume shear layer defined by the 
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Figure V-7. Freestream Velocity Distribution Mode II Secondary Throat 
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Task I bleed flow aerodynamic model was used. This approach assumes that the 
shear layer is not disturbed by its interaction with the wall. In the second 
case, the shear layer was assumed to be nearly destroyed by its interaction 
with the wall, such that the initial momentum thickness for the wall boundary 
layer is very small. Results reported for this case are for an initial thick- 
ness of .0254 mm (0.001 in.), compared with a shear layer momentum thickness 
of .861 mm (0.0339 in.). Further reductions in initial momentum thickness 
have little effect on the Stanton number at the location of the maximum free- 
stream mass velocity. Therefore, these two cases are considered to bound the 
actual interaction of the shear layer with the wall. 


Stanton numbers were defined from the skin friction coeffi- 
cient, using a Reynolds analogy factor of 1.25. This value is a compromise 
between the value of 1.19, calculated from the theory of Ref. 28 for a turbu- 
lent Prandtl number of 0.9, and the value of 1.32 defined by the Colburn 
factor Pr"2/3. Results from the boundary layer analysis are summarized 
below for the location of maximum freestream mass velocity: 


6 , 

Re^ X 10-4 
St 

hg, KW/m2-°K 
(Btu/i n.2-sec-°F) 

Heat Flux, MW/ m2 
(Btu/in. 2-sec) 

♦Defined by plume 


Initial Momentum Thickness, mm (in.) 


.861 (0.0339) 

.0247 

3.35 

.00179 

30.32 (.0103) 
86.1 (52.7) 


.0254 (.001) 

.00244 

.331 

.00263 

44.45 (.0151) 
126.3 (77.3) 


shear layer 


These results compare with a maximum Mode I heat transfer 
coefficient of 28.38 KW/m2-°K (0.00964 Btu/in.2-sec-°F) and heat flux 
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of 85.46 MW/m2 (52.3 Btu/in. 2-sec) in the secondary throat region. Thus, 
when the initial momentum thickness for the zero bleed flow boundary layer is 
equal to that from the plume shear layer, the maximum Mode I and Mode II heat 
fluxes are essentially equal. A similar result for Mode II is obtained by 
scaling the maximum heat flux from the Task I BLIMP run for nominal bleed flow 
according to the local VNAP wall pressures (a procedure suggested by various 
flow reattachment data): 

0 =0 
no bleed ’"nom. bleed 

BLIMP 

= 57 Btu/in.^-sec (93.14 MW/m^) 

The nominal bleed flow BLIMP was started with an initial momentum thickness 
from the corresponding plume shear layer (D streamline) since the plume flow 
is very nearly parallel to the wall. When the initial momentum thickness for 
the zero bleed flow boundary layer is very small, the maximum Mode II heat 
flux is almost 50% higher than that for Mode I. The ability of the secondary 
chamber cooling system to accommodate this higher heat flux is discussed in 
Section V,C,1. 

B. CHANNEL DESIGN PROCEDURES 

Two computer programs for regenerative cooling were used for the 
present studies: one with the dual throat geometric details built in, the 

other for the dual expander. These programs were derived from the SCALER 
program for conventional geometry which was developed specifically for para- 
metric design studies. With these programs, it is economically feasible to 
generate a relatively large number of parametric design points and still 
obtain a detailed, multi-station analysis of a rectangular channel at each 
design point. The chamber geometry and the local gas-side heat transfer 
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coefficients and coolant heat loads are scaled from reference input to other 
thrusts and chamber pressures. At each station, the programs iterate to 
determine the channel depth required to satisfy a gas-side wall temperature 
limit, which can be specified as a function of closeout wall temperature 
consistent with cycle life and creep criteria (Figure V-8). The only simpli- 
fying assumption is that gas-side wall temperature differences between the 
reference input and scaled cases have a negligible effect on gas-side heat 
transfer coefficients and heat loads. Normally, gas-side wall temperature 
limits are well known in advance so that local reference gas-side heat 
transfer analyses can be run at appropriate wall temperatures. 

Two-dimensional conduction effects in rectangular channels, as 
well as the spatial variations of the coolant heat transfer coefficient, are 
approximated as shown in Figure V-9. The hot wall adjacent to the channel, 
the entire land, and that part of the external wall adjacent to the channel 
are represented as fins; that part of the external wall adjacent to the land 
is assumed to be isothermal. A simple resistance network represents the hot 
wall adjacent to the land and connects the hot-wall fin to the land fin. Dif- 
ferent coolant heat transfer coefficients are associated with each fin. That 
for the hot wall is based on the average of the centerline (TWL2) and corner 
(TCORN) wall temperatures, while the coefficient for the external wall is 
based on its centerline wall temperature (TBS). The land coefficient is 
merely a weighted average of the other two. In this way, it has been possible 
to simulate the wall temperatures predicted by the SINDA network analyses very 
closely. 


Channel geometry constraints imposed in the present study due to 
fabrication and flow distribution considerations were as follows: 
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Figure V-8. Cycle Life/Creep Wall Temperature Criteria 



Figure V-9. Schematic of Wall Conduction Model 
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Minimum wall thickness, mm (in.) 
Minimum land width, mm (in.) 
Minimum channel width, mm (in.) 
Minimum channel depth, mm (in.) 
Maximum channel depth/width 


.635 (.025) 
1.016 (.040) 
1.016 (.040) 
.889 (.035) 
5:1 


Wall thicknesses were also required to satisfy the strength criteria of Figure 
V-10. The number of coolant channels was determined by utilizing the above 
minimum dimensions for both the land and channel widths at either the throat 
or the location in a circuit of the minimum radius. 


Channel layouts were limited to combinations of regions with con- 
stant channel width or constant land width (straddle milling). Increases in 
flow area relative to the throat beyond those possible with channel depth 
increases were accomplished by constant land width regions or by step changes 
in channel width. At the pressure levels associated with the present baseline 
designs, only very small increases in channel width .102-. 127 mmg (.004 - .005 
in.) above the 1.016 mm (0.040 in.) minimum were possible without increasing 
the wall thickness above .635 mm (0.025 in.). Details of the channel layouts 
for each concept are presented in subsequent sections. 


C. DUAL THROAT 


1. Baseline Design 


As shown schematically in Figure V-11, three parallel cooling 
circuits are used in all dual throat designs presented herein. A separate 
circuit is included for the inner annulus, lip region, and part of the primary 
nozzle; in previous dual throat analyses (Ref. 20), these regions were 
included in the primary chamber circuit. Since the primary chamber dictates 
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the overall system pressure drop, use of a separate circuit to accommodate the 
pressure drop required for lip cooling significantly reduces the overall pres- 
sure drop. The primary nozzle area ratio, at the interface between 
these circuits was a parameter in the baseline design study; a maximum value 
of 1.92 is possible from a mechanical design standpoint. Note that the 
secondary circuit extends to an area ratio of 7.8, after which an oxygen- 
cooled tube bundle would be used for the remainder of the secondary nozzle. A 
parametric design study of such a tube bundle v/as presented in Ref. 2. 

Channel layout details for the baseline design are given in 
Table V-I. The primary circuit utilizes straddle milling in part of the con- 
vergent section in order to increase the channel width from the throat values. 
In order to avoid the channel aspect ratio limit in the barrel, it was neces- 
sary to increase the channel width above that consistent with the strength 
criteria for the minimum wall thickness. Therefore, the barrel wall thickness 
was increased to .762 mm (0.030 in.). In the secondary circuit, the minimum 
wall thickness of .635 mm (0.025 in.) was used in the barrel, which limited 
the channel width to 1.143 mm (0.045 in.), i.e., only .127 nm (0.005 in.) 
greater than the throat channel width. In order to increase the flow area 
awcy from the throat region more rapidly than provided for by straddle 
milling, a step change in width occurs in the convergent section. This barrel 
design was not optimized for Mode I, since throat cooling in Mode II dictates 
most of the circuit pressure drop and since the total pressure drop does not 
appear to be a problem. The throat land width is maintained in the first part 
of the secondary nozzle, with the wall thickness increasing to 1.270 mm (0.050 
in.) to accommodate the increased channel widths. A channel width consistent 
with the strength criteria for the 1.270 mm (0.050 in.) wall thickness is 
maintained in the second part of the secondary nozzle. A uniform channel 
width of 1.016 mm (0.040 in.) was maintained in the lip circuit. This is 
necessary in the primary nozzle segment of the circuit in order to provide the 
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TABLE V-I 


'vl 

U1 


BASELINE CHANNEL LAYOUT SUMMARY (SI UNITS) 


DUAL THROAT DUAL EXPANDER 


Primary 

Secondary 

Lm 

Primary 

Secondary Inner 

Primary 





Annulus 

Nozzle 


Barrel 










Channel Width (W) 



1.346 

1.143 


1.118 

1.016 

1 .118 


Land Width (L) 



1.854 

2.108 


2.083 

1 .295 

1.346 


Wall Thickness (t) 



.762 

.635 


.635 

.635 

.635 


Convergent Section 

(1) 

L = 

1 .854 

Step W 

W=1.016 

Step W 

W=1 .016 

W=1.118 



(2) 

W = 

1.016 







Throat or Min. Radius 










Channel Width 



1.016 

1.016 

1.016 

1.016 

1.016 

1.118 

1.016 

Land Width 



1.016 

1.016 

1.016 

1.016 

1.016 

1.575 

1.016 

Wall Thickness 



.635 

.635 

.635 

.635 

.635 

.635 

.635 

Nozzle 


W = 

1.016 

(1) L=1.016 

(2) W=2.007 

H=1 . 01 6 

L=1.016 

L=1.016 

W=1.118 

W=1 .016 





t=1.270 
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TABLE V-: 
BASELINE CHANNEL LAYOUT 

DUAL THROAT 

Primary Secondary 


Barrel 


Channel Width (W), in. 


.053 

.045 

Land Width (L) , in. 


.073 

.083 

Wall Thickness (t), in. 


.030 

.025 

Convergent Section 

(1) 

L = .073 

Step W 


(2) 

W = .040 


Throat or Min. Radius 




Channel Width, in. 


.040 

.040 

Land Width, in. 


.040 

.040 

Wall Thickness, in. 


.025 

.025 

Nozzle 


W = .040 (1) 

L = .040 



(21 

W = .079 




t = .050 


(cont. ) 

SUMMARY (ENGLISH UNITS) 


DUAL EXPANDER 

Lip Primary 

Secondary 

Inner 

Annulus 

Primary 

Nozzle 

.044 

.040 

.044 


.082 

.051 

.053 


.025 

.025 

.025 


W = .040 Step W 

W = .040 

W = .044 



.040 

.040 

.040 

.044 

.040 

.040 

.040 

.040 

.062 

.040 

.025 

.025 

.025 

.025 

.025 

W=.040 

L=.040 

L = .040 

W = .044 

W = .040 
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most effective channel geometry in the lip region. However, the inner annulus 
segment could be straddle-milled, which would provide a small reduction in 
circuit pressure drop. 

Figure V-12 shows the pressure drops required for the primary 
circuit for two inlet area ratios as a function of the fraction of the hydro- 
gen flow used in the circuit. Each point on these curves represents a dif- 
ferent channel depth profile, as defined by the cycle life criterion of Figure 
V-8 and the maximum depth/width ratio of 5. As the hydrogen flow fraction 
initially increases, the local channel depth increases for a fixed coolant 
velocity; however, the increased depth results in a larger effective coolant 
surface area, thereby reducing the required coolant velocity and allowing an 
even greater channel depth. The reduced velocity and increased hydraulic 
diameter result in a significant decrease in required pressure drop. However, 
as the coolant flow continues to increase, substantial portions of the chamber 
reach the channel aspect ratio limit resulting in a pressure drop increase. 
Figure V-12 indicates that reducing the inlet area ratio of the primary cir- 
cuit reduces its pressure drop; the coolant Mach number, which exceeds 0.30 
for higher area ratios, is also reduced. In addition, the hydrogen flow frac- 
tion at the minimum pressure drop decreases slightly as the inlet area ratio 
is reduced, thereby providing more flow for the other circuits. 

Figure V-12 also includes the required pressure drop charac- 
teristics for the secondary circuit for Mode I operation. A hydrogen flow 
fraction of 0.30 was selected for the secondary circuit in order to provide a 
large pressure drop margin for redesign of the throat region for Mode II oper- 
ation with a bleed flow malfunction. The Mode II secondary throat analysis 
discussed previously (Section V,A,2) indicated that local heat fluxes in the 
primary plume impingement region may be as much as 50% higher than for Mode I. 
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A composite design with revised channel depths in the Mode II impingement 
region was developed for the worst-case boundary condition of Section V,A,2; 
Figure V-13 compares the Mode I and II heat flux profiles. The pressure drop 
for the composite design during Mode I operation is 4.45 MPa (645 psi), which 
is well below the primary circuit pressure drop of Figure V-12. A detailed 
design of the lip region circuit was completed for an interface area ratio of 
1.58, assuming a primary circuit flow fraction equal to the optimum value of 
0.40. In all cases, the total hydrogen flow fraction available for cooling 
was assumed to be 0.94, based on the gas generator/staged combustion mixed 
cycle studies of Ref. 2. Therefore, the flow fraction available for the lip 
region circuit was 0.24. Table V-II summarizes the baseline designs for all 
three cooling circuits. Since the lip circuit pressure drop is well below 
that of the primary circuit, it is apparent that further reduction in the 
interface area ratio is desirable in order to minimize the cooling system 
pressure drop. 


2. Parametric Study 


Chamber pressure and thrust split were varied independently 
in the dual throat parametric study. Primary chamber pressure* was varied 
from 20.68 to 34.47 MPa (3000 to 5000 psia) at the baseline thrust split. 
However, at 34.47 MPa (5000 psia), the coolant Mach number was impractical; 
therefore, the highest pressure for v/hich results are reported is 31.06 MPa 
(4500 psia). Mode I thrust from the secondary propellants was varied from 60 
to 80% at the baseline chamber pressures. The design procedure in each case 
was similar to that for the baseline, except that composite secondary designs 
including the effects of Mode II v/ere not generated. Primary and secondary 
circuit designs for Mode I v^ere defined as a function of hydrogen flow 

*Secondary chamber pressure was always 70% of the primary, as defined by 
aerodynamic considerations. 
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Figure V-13. Comparison of Secondary Throat Region Heat Fluxes 
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TABLE V-II 


DUAL THROAT BASELINE COOLING CIRCUITS 
INTERFACE AREA RATIO =1.58 



Hydrogen 

Pressure 

A 

T. 

Max. 

Minimum 


Flow 

Drop 


b 

Mach 

Channel 


Circuit 

Fraction 

MPa (Psi) 

°-i 

K (°F) 

No. 

Depth, 

mm (in.) 

Primary 

0.40 

8.07 (1170) 

115 

(207) 

0.29 

2.184 

(.086) 

Inner Annulus 

0.24 

4.41 (640) 

112 

(202) 




Annulus 


.483 (70) 

36 

(64) 

0.04 

5.080 

(.200) 

Lip 


2.69 (390) 

5 

(9) 

0.11 

2.032 

(.080) 

Primary Nozzle 


1.24 (180) 

72 

(129) 

0.12 

3.200 

(.126) 

Secondary (Composite) 

0.30 

4.45 (645) 

363 

(654) 

0.31 

1.118 

(0.044) 
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fraction, as shown previously in Figure V-12 for the baseline. The secondary 
circuit flow fraction was selected to provide a large pressure drop margin 
relative to the minimum primary circuit pressure drop; this margin is avail- 
able to provide smaller channel depths in the Mode II impingement region. The 
remaining coolant flow was divided between the primary and lip region cir- 
cuits, with the primary flow selected to give the minimum pressure drop 
whenever possible. For the higher thrust split and chamber pressure cases, it 
was necessary to select a primary flow less than optimum, resulting in a pres- 
sure drop up to 11% higher than minimum. Due to the time required to design 
the lip region circuit, the flow split between the primary and lip circuits 
was not optimized by balancing the required pressure drops. 

The parametric designs are summarized in Table V-III. 

Figures V-14 and V-15 show the coolant pressure drop and Mach number as a 
function of primary chamber pressure and thrust split, respectively. At the 
highest pressure and thrust split, the coolant Mach number has reached 0.37; 
therefore, further increases in these parameters are not practical with the 
present channel layout. 

D. DUAL EXPANDER 

1. Baseline Designs 

As shown in Figure V-16, the initial design considered was a 
hydrogen regeneratively cooled system with three parallel cooling circuits. 
Note that these circuits are conceptually the same as for the dual throat 
designs. Although the lip region cooling problem is less severe for the dual 
expander, the inner annulus circuit is much longer and includes a throat 
region. 
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TABLE V-III 


DUAL THROAT PARAMETRIC DESIGN SUMMARY 


Primary Mode I 

Pc Thrust 

MPa (psia) Split 

Circuit 

Hydrogen 

Flow 

Fraction 

Pressure 

Drop 

MPa (Psi) 

°K (°F) 

Max. 

Mach 

No. 

Minimum 

Channel 

Depth, mm (in.) 

27.58 (4000) 60/40 

P 

.40 

7.060 (1024) 

94 (169) 

.24 

2.591 (.102) 


L 

.24 

4.116 (597) 

77 (139) 

.14 

3.048 (.120) 


S 

.30 

1.372 (199) 

26 (466) 

.09 

4.089 (.161) 

75/25 

P 

.35 

9.384 (1361) 

149 (268) 

.33 

1.753 (.069) 


L 

.24 

4.585 (665) 

127 (229) 

.12 

1.600 (.063) 


S 

.35 

2.289 (332) 

382 (688) 

.12 

2.896 (.114) 

80/20 

P 

.35 

10.85 (1574) 

166 (299) 

.37 

1.575 (.062) 


L 

.20 

5.578 (809) 

168 (303) 

.15 

1.016 (.040) 


S 

.40 

2.682 (389) 

402 (724) 

.13 

2.540 (.100) 

20.68 (3000) 70/30 

P 

.30 

3.220 (467) 

146 (262) 

.20 

2.692 (.106) 


L 

.30 

2.020 (293) 

97 (175) 


3.429 (.135) 


S 

.30 

1.145 (166) 

373 (672) 

.09 

4.216 (.166) 

31.06 (4500) 

P 

.40 

13.12 (1903) 

119 (215) 

.37 

1.778 (.070) 


L 

.24 

6.709 (973) 

109 (197) 

.14 

1.600 (.063) 


S 

.30 

3.330 (483) 

364 (656) 

.13 

2.388 (.094) 








in 

o. 


CL 

o 

Q 

LiJ 

or 

=D 

CO 

CO 

UJ 

cc 

a. 



Figure V-15. Dual Throat Parameter Study Effect of Thrust Split on Coolant Pressure Drop 
and Mach Number 
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V, D, Dual Expander (cent.) 


Channel layout details for the baseline design are included 
in Table V-I. The primary circuit layout is similar to the dual throat secon- 
dary circuit. Channel aspect ratios in the barrel are within the design 
limit, so there was no motivation for increasing the wall thickness to provide 
a wider channel. A uniform channel width of 1.016 mm (0.040 in.) was main- 
tained in the barrel and convergent sections of the secondary circuit because 
of the small channel depths required. A uniform channel width was also used 
over the entire inner annulus surface. It was found that a 1.118 mm (0.044 
in.) width provided a lower pressure drop than the minimum allowable width of 
1.016 mm (0.040 in.). The latter value was maintained in the primary nozzle 
section of the inner annulus circuit. This section defined the number of 
channels in the circuit and is shown separately in Table V-I. 

Two area ratios, were considered for the interface between, 
the primary circuit inlet and the inner annulus circuit outlet: 2.34, which 

represents the maximum from a mechanical design standpoint, and 1.66. Figure 
V-17 shows the primary circuit pressure drop characteristics for these two 
interface area ratios, along with the secondary circuit characteristics. The 
coolant flow fraction for the latter is set at 0.24 by the minimum channel 
depth constraint. Inner annulus flow requirements for Mach number control, 
the remaining hydrogen flow fractions available for the primary circuit, and 
the corresponding primary circuit pressure drops from Figure V-17 are as 
f ol 1 ows : 


Interface Inner Annulus 

Area Ratio Flow Fraction 


Primary 
Flow Fraction 


Primary 
^P. MPa.(psi) 


2.34 0.35 

1.66 0.40 


0.35 11.86 (1720) 

0.30 12.86 (1860) 


The reduced flow available for the primary circuit with the lower interface 
area ratio offsets the reduced length and increases the primary circuit 
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Figure v-17. Dual Expander Baseline Cooling Circuit Pressure Drop Characteristics 
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V, D, Dual Expander (cont.) 


pressure drop. In both cases, the inner annulus flow required for coolant 
Mach number control forces the primary flow to be less than optimum. Table 
V-IV summarizes the circuit designs for an interface area ratio of 2.34. In 
view of the Mach numbers in the primary and inner annulus circuits and the 
channel depth in the secondary circuit, it is clear that cooling limitations 
have been reached in both the primary and secondary chambers. The secondary/ 
primary chamber pressure ratio was selected in anticipation of such a result. 
However, a slightly lower ratio would have allowed the primary circuit pres- 
sure drop to be optimized. 

Figure V-17 indicates that the hydrogen pressure drop can be 
reduced to 10.69 MPa (1550 psi) by increasing the primary circuit flow frac- 
tion to 0.40 with an interface area ratio of 1.66. This can be accomplished 
by using oxygen to cool the secondary circuit. Table V-V summarizes such a 
design; it should be noted that the secondary circuit channel design has not 
been optimized for oxygen cooling. 

Figure V-18 shows the effect of transpiration cooling the 
throat region of the primary chamber for an interface area ratio of 1.66. The 
upstream end of the transpiration region was varied, with the downstream end 
fixed at area ratio 1.44, which is 3.74 cm (1.47 in.) aft of the throat. 

Since the secondary chamber was assumed to be hydrogen cooled, only 30% of the 
total hydrogen flow was available for the two primary cooling circuits. Con- 
sequently, the significant coolant flow fractions required for the longer 
transpiration sections severely limit the available regenerative cooling flow 
and result in increased pressure drops. The minimum regenerative pressure 
drop in this case is about 11.03 MPa (1600 psi), which is higher than the 
optimum value of Figure V-17 that can be achieved by oxygen cooling of the 
secondary. 
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TABLE V-IV 


DUAL EXPANDER BASELINE COOLING CIRCUITS 
INTERFACE AREA RATIO =2.34 


Circuit 

Primary 
Inner Annulus 
Annulus 
Lip 

Primary Nozzle 


Hydrogen 

Pressure 


Max. 

Minimum 

Flow 

Drop 

D 

Mach 

Channel 

Fraction 

MPa (Psi) 

U. 

o 

o 

No. 

Depth, mm (in.) 

0.35 

11.86 (1720) 

157 (282) 

0.30 

2.159 (.085) 

0.35 

6.45 (935) 

182 (328) 




4.48 (650) 

161 (289) 

0.29 

1.321 (.052) 


.93 (135) 

1.7 (3) 

0.09 

5.08 (.200) 


1.03 (150) 

20 (36) 

0.10 

5.08 (.200) 

0.24 

4.76 (690) 

237 (427) 

0.27 

9.14 (.036) 


Secondary 



TABLE V-V 


Ci rcuit 

Primary 
Inner Annulus 
Annulus 
Lip 

Primary Nozzle 


DUAL EXPANDER BASELINE WITH OXYGEN-COOLED SECONDARY 
INTERFACE AREA RATIO =1.66 


Coolant 

Flow 

Fraction 

Pressure 

Drop 

MPa (Psi) 

°K 

0.40 H 2 

10.69 (1550) 

132 

0.40 H 2 

8.0 (1160) 

166 


4.62 (670) 

141 


1.48 (215) 

1.7 


1.90 (275) 

23 

0.60 O 2 

10.0 (1450) 

136 


(°F) 

Max. 

Mach 

No. 

Minimum 

Channel 

Depth, mm (in.) 

(237) 

0.28 

2.39 (.094) 

(299) 



(254) 

0.29 

1.65 (.065) 

(3) 

0.11 

5.08 (.200) 

(42) 

0.12 

5.08 (.200) 

(244) 

0.23 

2.79 (.110) 


Secondary 
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V, D, Dual Expander (cont.) 


2. Parametric Study 

The dual expander parametric study was based on the hydrogen 
regeneratively cooled concept of Table V-IV. Independent variations of thrust 
and Mode I thrust split at the baseline chamber pressures were investigated, 
along with various chamber pressure combinations at the baseline thrust and 
thrust split. Table V-VI indicates the design points attempted and identifies 
those for which chamber designs could not be obtained; characteristics of the 
successful designs are summarized in Table V-VII. 

In all cases with increased chamber pressures or Mode I 
thrust split relative to the baseline, designs could not be obtained due to 
insufficient hydrogen flow for the cooling of three circuits. A reduction in 
secondary chamber pressure to 10.34 MPa (1500 psia) with the primary chamber 
pressure unchanged resulted in convergence failure due to large pressure drops 
in the primary circuit. The latter results from the reduction in coolant 
inlet pressure associated with the reduction in secondary chamber pressure. 
From the baseline and parametric results, it appears that the present baseline 
design point is near the optimum chamber pressure ratio and is very close to 
the chamber pressure vs thrust split characteristic which represents a prac- 
tical cooling limit. Future parametric studies should define additional 
points on such a characteristic and consider smaller chamber pressure incre- 
ments to optimize the pressure ratio. 

E. SUMMARY OF RESULTS 

A satisfactory dual throat baseline design was obtained by using 
three parallel cooling circuits. The primary circuit coolant flows from area 
ratio 1.58 in the primary nozzle to the primary injector; this circuit defines 
the system pressure drop, so its coolant flow fraction is selected to provide 
the minimum pressure drop. A separate circuit cools the inner annulus 
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TABLE V-VI 


DUAL EXPANDER PARAMETRIC STUDY DESIGN POINTS 




Mode I 

Primary 


Thrust 

Thrust 

Pc 

MN 

(10® IbF) 

Split 

MPa (psia) 

2.67 

(0.6) 

59/41 

41.37 (6000) 




27.58 (4000) 




48.26 (7000) 



44/56 

41.37 (6000) 



68/32 


4.44 

(1.0) 

59.41 


6.67 

(1.5) 




Secondary 

Pc 


MPa (psia) 

Result 

10.34 

(1500) 

Convergence Failure - Primary 

24.82 

(3600) 

Insufficient Coolant 

13.79 

(2000) 

Successful Design 

24.13 

(3500) 

Convergence Failure - Inner Annulus 

20.68 

(3000) 

Successful Design 


Insufficient Coolant 
Successful Design 
Successful Design 


TABLE V-VII 


DUAL EXPANDER PARAMETRIC DESIGN SUMMARY 


Chamber 

Total 

Mode I 


Hydrogen 

Pressure 


Max. 

Minimum 

Pressure 

Thrust 

Thrust 


Flow 

Drop 

D 

Mach 

Channel 

MPa (psia) 

MN 006 IbF) 

Spl it 

Circuit 

Fraction 

MPa (Psi) 

Li. 

o 

o 

No. 


27.58/13.79 (4000/2000) 

2.67 (0.6) 

59/41 

P 

.25 

3.96 (574) 

219 (394) 

.22 

2.82 (.111) 




lA 

.35 

2.70 (391) 

174 (314) 

.18 

2.14 (0.84) 




S 

.30 

.641 (93) 

172 (309) 

.09 

3.10 (.122) 

41.37/20.68 (6000/3000) 

2.67 (0.6) 

44/56 

P 

.25 

11.30 (1639) 

118 (213) 

.28 

2.03 (.080) 




lA 

.25 

7.67 (1112) 

153 (276) 

.21 

2.13 (.084) 




S 

.25 

1.41 (204) 

146 (263) 

.11 

2.95 (.116) 


4.44 (1.0) 

59/41 

P 

.40 

10.42 (1511) 

122 (220) 

.30 

2.87 (.113) 




lA 

.30 

5.96 (865) 

182 (327) 

.25 

1.52 (.060) 




S 

.24 

2.0 (290) 

194 (350) 

.15 

1.80 (.071) 


6.67 (1.5) 


P 

.30 

10.45 (1515) 

147 (264) 

.29 

2.87 (.113) 




lA 

.30 

7.36 (1068) 

167 (300) 

.22 

1.98 (.078) 




S 

.25 

1.613 (234) 

165 (297) 

.11 

2.82 (.111) 


V, E, Summary of Results (cont.) 


surface, the lip region, and the rest of the primary nozzle. The majority of 
the channel pressure drop in this circuit is required for the lip region, for 
which a three-dimensional SINDA model was developed. A counterflow circuit 
cools the outer surface of the secondary chamber. This circuit is a composite 
design considering both Mode I and Mode II operation; the channel design over 
a short region upstream of the secondary throat is defined by the primary 
plume impingement associated with a bleed flow malfunction in Mode II. An 
exact solution for the inviscid plume impingement problem was obtained with 
the VNAP program. This solution provided edge conditions for an integral 
boundary layer analysis which defined the wall heat transfer. The initial 
momentum thickness for the boundary layer analysis was varied in order to 
bound the problem, with the channel design based on the worst case; the 
maximum heat flux ranged from essentially the Mode I value to almost 50% 
higher. 


Table V-II summarizes individual cooling circuit characteristics 
for the dual throat baseline design. Although the fuel system pressure drop 
of 8.00 MPa (1160 psi) is set by the primary circuit, the coolant Mach number 
is slightly higher in the secondary circuit when the latter is designed to 
accommodate the maximum postulated Mode II heat flux with no bleed flow. 
Parametric study results indicate that it is the coolant Mach number rather 
than pressure drop which limits increases in chamber pressure and Mode I 
thrust split. Increasing the chamber pressures to 31.02/21.72 MPa. (4500/3150 
psia) or the thrust split to 80/20 results in a primary circuit Mach number of 
0.37. 


The dual expander designs are cooled with three parallel circuits 
similar to those of the dual throat. However, the inner annulus circuit is 
now a major circuit in terms of length and heat load, and it includes a throat 
region. Regenerative cooling of the baseline design point with hydrogen only 
is difficult. Coolant flows in the secondary and inner annulus circuits must 
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V, E, Summary of Results (cont.) 


be selected to avoid minimum channel depth and Mach number limits, respec- 
tively, and the remaining flow available for the primary is less than optimum. 
As a result, the system pressure drop of 11.86 MPa (1720 psi) is .896 MPa (130 
psi) above the minimum primary circuit value. Using oxygen to cool the secon- 
dary circuit and reducing the interface area ratio between the primary circuit 
inlet and the inner annulus circuit outlet allows the hydrogen pressure drop 
to be reduced to 10.69 MPa (1550 psi). These designs are sumamrized in Tables 
V-IV and V-V. Transpiration cooling the throat region of the primary chamber 
was not effective since the hydrogen flow fractions required for transpiration 
cooling, coupled with the flow requirements of the other circuits, severely 
limit the flow available for regeneratively cooling the remainder of the 
primary circuit and result in relatively high pressure drops. 

Parametric study designs could not be obtained with increased 
chamber pressures or increased Mode I thrust split, or with the secondary 
chamber pressure reduced to 10.34 MPa (1500 psia). The latter results from 
the high primary chamber pressure drop associated with reduced coolant inlet 
pressure. These results, along with the baseline design details, indicate the 
baseline design point is near the optimum chamber pressure ratio and is very 
close to the chamber pressure vs thrust split characteristic which represents 
a practical cooling limit. 
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VI. 


HOT-FIRE DATA ANALYSIS 


A. PRE-TEST ANALYSIS 

The overall objective of this task was to provide support for the 
hot-fire testing of the dual throat hardware by NASA/MSFC. ALRC's first 
objective was to review the MSFC's test plan and provide whatever consultation 
might be necessary during testing. Secondly, ALRC would review the hot-fire 
data and interpret test results by using the aerodynamic model and previous 
cold-flow data. Test data, would also be used to correlate dual throat nozzle 
Mode II performance. Thirdly, corrections to the aerodynamic model believed 
advisable were to be incorporated into the model, and further model improve- 
ments were to be recommended. 

The hardware tested was the same as that used during cold-flow 

testing (NAS 8-32666), (Ref. 1). For these tests, the small primary nozzle 

(e = 1.72) and the large secondary nozzle (e = 3.0) were used. A layout 
P ^ 

of the test hardware is shown in Figure VI -1. The proposed test plan 

called for testing to be conducted at three different nozzle spacings, using 
GH 2 as the bleed gas. Recommendations were made by ALRC for testing to be 
done with various bleed flow gas compositions and temperatures to verify the 
improvements made to the aerodynamic bleed flow model to simulate such condi- 
tions. Funding and time limitations did not permit testing a survey of bleed 
flow gases, although one series of tests, at the larger nozzle spacing, was per- 
formed with GN 2 as the bleed gas. The tests conducted are summarized in Table VI-I. 

The improved aerodynamic bleed flow model was used to define the 
"blowoff" condition to provide a safe first test operation. The predicted 
trend in bleed flowrate versus nozzle spacing is shown in Figure VI-2, for 
both a GH 2 and a GN 2 bleed. A significant influence of bleed flow 
properties on operating characteristics was predicted. These predictions 
indicated that one and one-half to two times the bleed flow would be 
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TABLE VI-I 


DUAL THROAT THRUSTER 
MS FC IN-HOUSE TEST SUMMARY 


TEST BLEED BLEEDFLOW Pc. DURATION. 

NO. RTP QA8 WsAVp,% "/F P8IA SEC 
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• DEFINED BLOW-OFF CONDITION TO PROVIDE 
SAFE 1ST TEST OPERATION 

• PREDICTED SIGNIFICANT INFLUENCE OF BLEED FLOW 
PROPERTIES ON OPERATING CHARACTERISTICS 


I 

Figure vi-2. Predicted Trend in Bleed Flowrates versus Nozzle Spacing 



VI, A, Pre-Test Analysis (cont.) 

required when using GN 2 as opposed to using GH 2 . This difference is due 
to the influence that GN 2 has on the shear layer mixing equations with its 
larger molecular weight and smaller specific heat (Cp). One item of concern 
was the large bleed flowrates predicted using GN 2 and for the large nozzle 
spacing using GH 2 . The assumptions used in developing the bleed flow model 
were expected to limit the model 's val idity and applicability to small bleed 
flows of approximately 10% or less. As noted in Section III, significant 
reductions in bleed flow from the blowoff condition can be made without 
measurably inducing a loss in total pressure. 

B. TEST RESULTS 

1. Data Evaluation 

Wall pressure profiles for the various test configurations 
are shown in Figures VI-3 through VI-5 for various nozzle spacings and in 
Figure VI-6 for N 2 bleed flow. In general, the pressure data are reasonable 
and similar to the previous cold-flow test results (Ref. 1 ). A comparison of 
the hot-fire and cold-flow pressure profiles is shown in Figure VI-7 for the 
large nozzle spacing (Le/Rtp = 5.06) with GN 2 bleed. The plume 
impingement point is shown to move toward the secondary throat with increasing 
bleed flow, as was seen in cold-flow testing. Also, the plume impingement 
shock strength is decreased with bleed flow, which is consistent with the 
cold-flow data and the aerodynamic bleed flow model predictions. Test results 
also show that increased nozzle spacing requires increased bleed flow to 
attach the plume to the secondary throat which is also consistent with both 
the cold-flow data and the model predictions. 

Although the hot-fire data, overall, were consistent with the 
cold-flow data, some anomalies did exist. All the test configurations had 
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Figure VI-3. Wall Pressure Profile, Nozzle Spacing = 5.062, GH, Bleed Flow 
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Figure VI-.4. Wall Pressure Profile, Nozzle Spacing = 3.979, GH, Bleed Flow 
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Figure VI-6. Wall Pressure Profile, Nozzle Spacing = 5.062, GNo Bleed Flow 
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VI, B, Test Results (cent.) 


some questionable pressure measurements at the secondary throat. The wall 
pressure would increase due to the plume impingement, decrease, and then 
increase again at the secondary throat before decreasing as a result of being 
expanded through the secondary nozzle. Cold-flow data show the secondary wall 
pressure to increase due to the plume impingement and then continually 
decrease as the gas exits through the secondary nozle. 

Another anomaly is seen in Tests 15, 18, 23, and 34, shown in 
Figures VI-6, VI-5, and VI-4, respectively. The pressure in the secondary 
chamber increases and then decreases prior to plume impingement. This trend 
was not seen during cold-flow testing and is probably the result of small 
pressure measurement errors. 

Test 8, shown in Figure VI-3, is also inconsistent with all 
other test results. Its measured pressure profile is consistent with a mini- 
mum bleed flow, yet the measured bleed flow was the largest with that con- 
figuation. This may be due to a leak in the bleed circuit after the flowrate 
measurement or an erroneous flowrate measurement. 

2. Model Comparison 

Although some uncertainties do exist with the test data, the 
experimental results have verified the bleed flow model's predictions. An 
example of this is shown in Figure VI-8. The bleed flow model predicts nearly 
exactly the recirculation pressure for a given bleed flow. Also, the model's 
prediction that more GN 2 bleed would be needed to achieve the same recircu- 
lation pressure as with GH 2 was substantiated with the hot-fire test 
results. 

The hot-fire results also verified the correlation of the 
plume scaling factor (FCTR) with nozzle spacing. Figure VI-9 shows that. 
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Figure VI-8. Experimental Results Have Verified Aerodynamic Model Predictions 
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VI, B, Test Results (cont.) 


for zero bleed flow, the bleed flow model predicts exactly the recirculation 
pressure. Some inaccuracy does exist between the data and predictions at the 
larger bleed flows for small nozzle spacing. This may have resulted due to 
inaccuracies in the plume correlation. Another important point verified is 
that a noticeable decrease in bleed flow from that required for blowoff may 
result in a very small stagnation pressure loss from plume impingement. 

3. Engine Performance 

Two means of measurement were used to determine engine per- 
formance, characteristic velocity (C*), and specific impulse (Isp). Table 
VI-II presents the method used to determine these parameters. The format for 
determining specific impulse is the same as that used during previous cold- 
flow testi ng (Ref. 1 ). 

Tests 1 through 7 used only the primary nozzle to check out 
the facility and determine combustion performance. The resulting C* effi- 
ciency proved to be very low (70-75%) probably due to the addition of film 
cooling holes in the outer ring of the primary injector. Testing with the 
dual throat configuration had similar low primary combustion performance, 
resulting in a correspondingly low Isp efficiency. As a result, this 
unusually large combustion inefficiency tended to mask the relatively smaller 
inefficiencies attributable to the dual throat configuration (i.e., boundary 
layer and divergence losses). 

Nevertheless, the performance determined from the hot-fire 
test was correlated using the simplified prediction methodology. Estimates of the 
divergence, boundary layer, and kinetic efficiencies were used to define a perfect 
injector performance for the test configuration. This value was compared to the 
measured Isp for a typical test to "back out" an energy release efficiency as shown on 
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TABLE VI-II 


HOT-FIRE TESTS PERFORMANCE MEASUREMENT 


CHARACTERISTIC VELOCITY 

MEASURED FLOWRATES 
MEASURED PRIMARY P^. 
MEASURED PRIMARY D^ 


C* 

meas 


meas \ 
W 

T meas 


C* 

meas 


^*ODE MR Test 
Pc Test 
Tp Test 


SPECIFIC IMPULSE 

f MEASURED FLOWRATES 

t MEASURED STAT1[C PRESSURE 

PROFILE 

.• MEASURED PRIMARY CHAMBER PRESSURE 


^vac ” meas ^ ^F 


X A. 

Primary ^Primary 


s 


pdA 


Secondary 

Wall 


sp 


vac 


vac ,meas 


W. 


T meas 



^sp vac, meas 

^sp ODE MR Test 
Pc Test 
T Test 



VI, B, Test Results (cent.) 


Table VI-III. This method resulted in an indicated energy release efficiency 
of 76%, which is of the same magnitude as the characteristic velocity effi- 
ciency (73%). Therefore, while the hot-fire results cannot be used to provide 
Vernier calibration of the dual throat performance model, it appears that this 
measured result is consistent with the model and that the unusually low spe- 
cific impulse is primarily the result of a low energy release (combustion) 
efficiency. 


I 
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TABLE VI-III 


ro 


PERFORMANCE COMPARISON OF HOT-FIRE TESTS/AERODYNAMIC MODEL 


TEST 9 DATA PERFORMANCE ANALYSIS 


Pc 

— 

456 psia 



MR 

- 

0.78 

I^PODE. 

401.8 

% BLEED 

- 

7.5% (H 2 ) 

"^DIV 

.985 

"t 

= 

3.28 Ibm/sec 



^vac 

= 

1261 IbF 

'^KIN 

1.00 

I^Pvac 

= 

295 sec 

aIsPbll, sec 

8.2 

C* . 
Primary 

= 

5570 ft/sec 

^^Pperfect Injector 

387.6 

% Isp 

= 

73 

^^Pmeasured 

295 

% C* 

= 

73 

’^ERE 

.76 
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3/10/80 


PROGRAM INPUT. AERODYNAMIC MODEL 


CARD 1; 


The 1st card to be input is used to provide a printed heading for 
the computer output. Columns 2 through 72 are available for text. This 
card must always precede each case. The remainder of the input data is 
read using NAMELIST with the name $DATA, as described below. 

$DATA 


Geometry; Option 1, Cylindrical Secondary Wall, see Figure 1. 
INPUT ITEM DEFINITION 


I WALL = 1 . 

RS 

RSP 

THETAP 

EPSP 


The secondary nozzle wall is a cylinder 
r*. radius of the secondary nozzle wall 
r*p, radius of the primary nozzle 
0p, exit cone half angle, primary nozzle 
cp, expansion ratio of the primary nozzle 


The items r* and r*p, above, can be input with units of inches, feet, cm, 
etc. The input item ep, above, is used only to determine the Mach number 
at the exit of the primary nozzle (not required if MEPI is input). 



THETAP 


4 


RS 


Figure i . Nozzle Geometry for Option 1 
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Geometry; Option 2, Full Geometry Option, see Figure ?• 


IWALL = 2, 


ECP 

RIP 

THIP 

RWTUP 

RSP 

RWTDP 

THETAP 

EPSP 


EC 

RI 

THETAI 

RWTU 

RS 

RWTD 

THETA 

EPS 


Complete wall geometry specified for both 
_pn ma ry_an[d jseconda ry_noz^l es 


Primary Nozzle 


E(. , primary nozzle chamber contraction ratio 

Rip, wall radius connecting chamber and inlet 
of primary nozzle, divided by 

01 p, Inlet angle for primary nozzle 


Ry , wall radius on upstream side of the 
throat, primary nozzlfe, divided by r* 

r*p, throat radius of primary nozzle , inches or cm 


R<j , wall radius on downstream side of the 
throat, primary nozzle, divided by r* 

0 °, exit cone half angle, primary nozzle 

Cp, expansion ratio of the primary nozzle 


Secondary Nozzle 


ec. secondary nozzle contraction ratio 

Rj, wall radius connecting chamber and 
inlet of primary nozzle , divided by r* 

0®I, inlet angle for secondary nozzle 

Rn, wall radius on the upstream side of the 
throat, secondary nozzle , divided by r* 

r*, throat radius of secondary nozzle, inches or cm 

Rrf. wall radius on downstream side of the 
throat, secondary nozzle , divided by r* 

0°, exit cone half angle, secondary nozzle 

e, expansion ratio of the secondary nozzle 




Geotnetry; Options 3, 4 or 5. Full Geometry Option with Spline Fit 
Exhaust Contours 

For the spline fit wall contour options, contour points are input to 
replace the conical exit cone, of either the primary nozzle, secondary nozzle, 
or both. Other inputs are the same as in geometry Option 2, described above, 
except that the inputs for ep, and/or e are not required. The exit cone 
half angle is interpreted as the spline contour attachment angle. The spline 
options are as follows: 

Primary Wall Secondary Wall 

I WALL =3, Spline Cone 

IWALL * 4, Cone Spline 

IWALL = 5, Spline Spline 

The spline coordinates are input as follows: 


X4P(2)* = 

Y4P(2) = 

N4P 

TH4P 

X4(2) - 

Y4(2) « 

N4 

TH4 


Xpj, axial coordinates for primary wall 
spline contour 

ypj, radial coordinates for secondary wall 
spline contour 

np, 1 *1, 2, ... np above, np ^ 20 

0®ep» spline contour exit angle for pri- 
mary nozzle, degrees 

xi. axial coordinates for secondary 
wall spline contour 

y-j, radial coordinates for primary wall 
spline contour 

n, 1 , ■ 1 , 2 . . . n above. N ^ 20. 

e®e, spline contour exit angle for secondary 
nozzle, degrees 


*The first point on a spline contour is automatically calculated by the 
program from the exit contour half angle, 0p®, or e°, thus each of these 
input vectors begins with item 2. 

All x,y coord inuLe vjlucs are normalized (i.e. divided) by the 
corres|)ondiof!| nozzle threat rsflius. 
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Rela t i ve Positioning of P riii ary and Secondary 
Nozzle^^^ (IWALL > 2 Cases Only) 


XLE Lp/r*p, axial distance from the exit plane 

of the primary nozzle to the throat plane 
of the secondary nozzle, divided by 
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INPUT 

FCTR 

IPNCH 

MEPI 

NW 

PBPO 

WLPRN 

WSWP 


ITEM 


FLOW PROPERTIES AND OP TIONS 
DEFINITION 


ASSUMED 

VALUE 


Factor for scaling plume shape. Used as a 1. 

multiplier for the parameter M./y in the Herron 
method. 


If IPNCH = 1, geometry inputs for TDK will be 0 

punched in the NAMELIST format . 

If IPNCH = 0, no punched output, only printed 
output. 

0 


13 


Initial Ruess for the pressure, ratio f®r the pluirie 
boundary. , If not input, the program will calcu- 
late a value corresponding to the area of the 
secondary nozzle throat divided by the area of the 
primary nozzle throat using the relationship for 
one-dimensional flow. 


M . (REAL type) exit Mach No., primary nozzle. 

If Mr <1 then values will be calculated from 
Ep P 

assuming one-dimensional flow 


n^, number of entries in WSWP 


< 50 


( 1 ) 


If a value other than zero is input, print of the 
secondary and primary nozzle wall tables will be 
suppressed (to be used on successive cases when 
the wall geometry does not change.) 

(Ws/Wp)i * 100, pernent bleed flow given for 

these values, i = 1, — , u for which base 

w 

pressure and recirculation temperature and 
primary mass fraction in the recirculation 
region , are solved. 


0 . 


0 , 1 , 2 , 3 , 4 , 5 , 

6,7,8,9,10, 

12,15 
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INPUT ITEM 


DEFINITION 


SUGGESTED 

VALUE 


TRTO 


Recirculation temperature ratio initial guess 
NOTE: TRTO is between TOBASE & 1.0 

1.0 

TOBASE 

Tob/^o 

Stagnation temperature ratio of bleed flow (secondary 
flow) to primary flow 


DELV 


Perturbation fraction for Newton solution procedure 

.0001 

DELVMX 


Maximum fractional change of Pu/P and T /T during 
solution 

.3 

EPSGGT 


Convergence criteria for bleed flow 

.0004 

EPSTOB 


Convergence criteria for T^j^/T^ 

.0004 

ITERPR 


Iteration print option, set equal to 1 for print 

0 

ITERMX 


Max # of iterations 

20 

CKN 


Nash Factor N 

.4 

LPBPO 

Logical 

Variable 

Find Flow-off Base Pressure (PBPOBF) if true 

True 

LPL0T 

Logical 

Variable 

Draw Plots if true 

False 

XMINPL 


Minimum axial coordinate for plot axis 


XMAXPL 


Maximum axial corrdinate for plot axis 


XDISPL 


Length of axis (for axial coordinate) 


YMINPL 


Minimum radial coordinate for plot axis 


YMAXPL 


Maximum radial coordinate for plot axis 


YDISPL 


Length of axis (for radial coordinate) 


PBPOBF 


Blow-off Base pressure to provide a limit during 
iteration or initial guess for finding blow-off 
pressure if LPBPO = .TRUE. 

.01 


(Note: This is then stored in PBPO). 
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INPUT ITEM 


DEFINITION 


SUGGESTED 

VALUE 


GGTBF 

CPF 

CP0 

MWF 

MW0 

ITYFR 

YFRX 


Blow-off Secondary Mass flow (percent) to provide a 
limit during iteration. (This is automatically set 
if LPBPO = .TRUE.) 

Heat capacity (constant pressure) of primary mass flow 
Heat capacity (constant pressure) of secondary mass flow 
Molecular weight primary mass flow 
Molecular weight secondary mass flow 

# of iterations to be taken on mass flow fraction 2 

during each base pass. 

Primary mass flow fraction in recirculation region .5 

first estimate for YFR. 
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3/10/80 


EXTENDED NOZ Z LE WALL OPTION 

Given an input nozzle wall contoui*, this option will calculate an 
extended nozzle wall contour that gives maximum thrust. This is 
accomplished by first calculating the flow field that is dependent .on the 
given input nozzle wall contour, i.e. the downstream boundary of the flow 
field is the right running characteristic extending from the nozzle lip 
to the flow axis. Shocks are ignored. Next, a circular arc of radius 
is attached tangent to the nozzle lip. The optimization options existing 
in the program are then used to calculate the maximum thrust wall contour. 
The following additional input items are required. 

IWL0PT = 1 for the extended wall option. The program assumes 

IWL0FT = 0 . 

RD = wall radius ratio (R/R^) of the circular arc, Rd, attached 

tangent to the nozzle lip. 


NOZZLE GEOMETRY 

The items RWTU, RWTD, and RSTAR are as previously defined, but apply 

to the input nozzle wall contour. See Figures la,b,c, and d for vjall contour 
options. 

IWALL = wall option flag . The value 

assigned this item determines 
the type of wall to be specified. 

0 => arbitrary contour 

1 => cone; 9, c 

2 => parabolic contour 

3 circular arc contour 

4 => spline fit contour 

5 => cone; r_, 

n n 


If IWALL = 0 (arbitrary contour) 

PW(1) = The wall coordinate points 

r j , Xj , are to be input 

sequentially; 

PW(l) “ ^1* ^2* ^2* • • • 

• • • f 0 • » 0 • , 
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a pair of zero coordinates 
marks the end of the table. 

A value of n S300 Is required. 
This table is normalized to 
the throat radius and the 
origin is at the throat so that 
necessarily r^ = 1, = 0. 


If IWAIJ 

r. 1 

(cone; 6, €) 

- 

RWTD 

= 

as above . 


THJI 


exit cone half angle (degrees). 


EPS 

= 

expansion ratio 


If IWALL 

= 2 or 3 

(parabolic contour or circular 




arc contour) 


RWTD 

.= 

as above 


THJI 

= 

wall Inflection angle (degrees) 


RWMAX 

. = 

nozzle end point radial 

coordinate, r 
n 


ZWMAX 


nozzle end point axial 
coordinate, x^ 


If IWALL 

= 4 

(spline fit contour) 


As in IWALL options 1,2, and 3 the program assumes the throat to begin 


with a circular arc downstream of the minimum point. The final point on this cir- 
cular arc Is located at an inclination of THJI degrees. The contour exit angle, THE 
degrees, must also be input. The contour to be spline fit between THJI and THE 
must be input in tabular form. The first point on the spline fit is automatically 
located by the program on the circular arc at the point of inclination THJI. 


RWTD 

=' 

As above 

* • 

THJI 


RWTD and the angle THJI define 
the end point of the circular arc 
and the beginning of the spline fit 
(degrees) 


NRZS 

= 

n, the number of points to be spline 
fit. n s t|o is required. 


PWRS(2) 

= 

radial wall coordinated to be spline 
fit. 




PWRS(2) - rj, r^ 
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PWZS(2) 


THE 

If IWALL = 5 

RWTD 

RWMAX 

ZWMAX 


axial wall coordinates to be spline 
fit, 

PWZS(2) = X , X . . . , X 

nozicle exit angle 

(cone; x^) 

As above 
As above 
As above 
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RWTD 


RSTAR 

Figure la: IWALL = .1 geometry; cone 
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APPENDIX B 

DERIVATION OF DUAL EXPANDER EQUATIONS 


230 



DERIVATION OF AREA EQUATION 


A = 2it 6 R 



Rj^ = Outer Radius 

R|^ = Inner Radius 

R|j^ = Radius @ Mid Point 

6 = Gap 

0 = Angle Gap Inclined 


The projected area of the annulus can be found by 

\ = . (R„2 - 

The annulus area is then 


A = Ap 

cos 0 


A = 

cos 0 


IT (Rf^ - R|_) (Rj^ + R^) 
cos 0 


2 



6 cos 0 


Substitute; 

A = n(6cos 0) 2 R^ 
cos 0 

A 2*fr 6 R|^ 
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DERIVATION OF INCLINATION ANGLE EQUATION 


Cos 0^ = 2 [R^ - R^^ - cos 0^] 


From Figure B-1 




cos 0, 


"’Em ^ '"El 


Rc ~ ^ ^E " *”e 

‘'c ^ "h '"c ‘'L 


E, = Rp + D, cos 0, 


Tr = 2 Rr - (Rc + D, cos 0. ) 
Eh E ' E L L 


"C cos 0c 


[(2 Rp - Rf - \ cos 0, )' 

c p 


(Rg + D|_ cos 0|_)‘^] 


Multiplying out and cancelling terms results 


Ac = [4 R^c - ^ Rr Re - ^ Rr D, cos 0, ] 

“"c cos 0g C C ‘^P “^c L L 


4 7T Rr. 


cos 0g c 


[Rp ■ R^ " \ ®L^ 


cos 0 , 


4 17 R, 


Ej. “ il|^ cos 0|^1 

c p 
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Substitute: 




234 






From Figure B-1 








AREA RATIO RELATIONSHIPS FOR DUAL EXPANDER 


MODE I: 


Primary flow streamtube area ratio (£„ ) 

^I 



Secondary flow streamtube area ratio (cj. ) 

^I 


'"S 


I 



Overall engine area ratio (e) 
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MODE II: 


Secondary nozzle area ratio ( e,. ) 



^'e 

Ae + 

Pi 

^I 



"ts 



A 

Pi 

^ + 

Si 






'Pi 

\s 
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Figure B-3 


APPENDIX C 

INPUT DESCRIPTION, DUAL EXPANDER SHEAR LAYER MODEL 
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INPUT DES CRI PTION 
DUAL EXPANDER SHEAR LAYER MODEL 


The purpose of this computer program is to calculate the 
base pressure that will occur when a uniform supersonic 
flow expands around the corner of a bluff body. An 
example is shown in Figure 1, which corresponds to the 
Dual Expander engine during Mode II operation. 



The method of Korst is used to predict the base pressure. 

The method applies only to flows that are not fully separated. 
The shear layer must approach the flow axis and be turned at 
the so called "nec-k region". The program allows a bleed 
flow to be specified, provided this secondary flow is not 
so great that the flow becomes fully separated. The pro- 
gram allows the molecular weight, specific heat, and total 
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enthalpy of the secondary stream (bleed) to be different 
than that of the primary stream. The method of using this 
program is described below. 

Solutions can be obtained for the following situations: 

1 ) zero bleed flow, 

2) bleed flow with properties identical to the 
properties of the primary stream, 

3) bleed flow with properties different than the 
properties of the primary stream. 

Case 3), above, is the general case solved by the program 

and requires the determination of three unknown variables 

(P./P H /H „,Y ) of three unknown variables such that 

three functional constraints are satisfied. Cases 1) and 

2), above, are special cases of case 3) and require the 

determination of only one variable (P./P ) such that one 

bop 

function is satisfied. In all three cases, initial esti- 
mates of the unknown variables must be supplied. These 
estimated values must lie within the radius of conver- 
gence Newton's method. Estimates can be obtained from 
solutions previously found. Each run on the computer will 
use its final values as the first estimate for the suc- 
ceeding case. This allows parametric studies to be run 
efficiently by the systematically changing of input. 

The program also uses bounds for the unknown variables 

so as to enhance convergence of the numerical method. 

Variables H /H „ and Y are automatically bounded by the 
r op pr 

program and have caused no difficulties thus far. The 

variable P./P , however, requires carefull consideration, 
b op ^ 

A satisfactory lower bound, 10 , is assumed; but the 

upper bound, max, must be input by the uses as 

described in the following steps. 
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step 1) Using y_ and the approach flow Mach No., M , 
look up^the Prandtl- Meyer angle, v , in tRe 
Gas Tabl es . ^ 


Step 2) Calculate the Prandtl -Meyer angle, v, that 
will turn the flow parallel to the axis by 
adding the approach flow angle to v , i. e.. 


Step 3) Look up the flow pressure ratio for parallel 
flow, P/P , corresponding to v. (^k/P „) ’”3’' 
must be l?ss than this value, i.e., 

(Pb/Pop) max = .95 (P/Pq) 

A good estimate for the base pressure ratio is 



One card is read and printed to label the program output. 
This card is followed by the NAMELIST input data set des- 
cribed below, followed by a listing of cards for a sample 
case. 


NAMELIST 


EXPLANATION AND COMMENTS 

UNITS 

DEFAULT 

NAME 




VALUE 

$DUALEX 





CPS = 

C Q 1 

Specific heat at constant pressure 

ca 1 / gm °K 

- 


Ps 

for the secondary (bleed) flow. 



CPP = 

Cn » 

Specific heat at constant pressure 

ca 1 / gm °K 

- 


Pp 

for the primary flow. 



AMWS = 


Molecular weight for the secondary 

- 

- 


*"s 

(bleed) flow. 

i 


AMWP = 


Molecular weight for the primary 

] 




flow . 

1 

1 


H0S = 


Total enthalpy for the secondary 

cal/gm I 

- 


°S 

(bleed) flow. = C„ T„ 

0^ p 0^ 

s s s 


1 

H0P = 


Total enthalpy for the primary 

- - P P — P 

cal /gm 


P0P = 


Total pressure for the primary 

ps i a 

- 


% 

flow. 



WSWP = 

w /w 

, Ratio of the mass flow rate of 

none 



S 

P the secondary (bleed) flow to 





the mass flow rate of the pri- 
mary flow 

1 



SRB = 

'"h ’ 

Radiusofbase 

inches 

- 

! 

1 

D 


(or cm) 


RWRB = 

j 


. , Wake radius ratio ’ 

D 

1 

none 

1 

.5 

AREA = 

The 

1 

throat area of the primary stream 

. 2 

, 2, 

1/ IT 



1 

(or cm ) 


THETAB = 

®b' 

Approach flow angle, i.e. lip angle 

degrees 




of the lower wall of the annular 
nozzle. 



AMP = 

Mp, 

Mach No. of the primary flow 

none 

- 
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NAMELIST 

N AME 

PBP0P1 = 
PMAX = 

HRH0P1 

YPRl = 

AN = 

GC0N = 
EJCNV = 
WF(1) = 

EPSI = 
NUN = 


PERTV = 


EXPLANATION AND COMMENTS UNITS DEFAULT 

VALUE 

P./P , 1st estimate for the ratio of none 
p base pressure to the total pres- 
sure of the primary flow. 


(Pb/Po )'"3x, upper bound for Pf^/P^ • none 1. 

P Use procedure previously P 
described to calculate this 
value. 


Hr/H , 1st estimate for the ratio of none 1. 

the enthalpy of the gas in the 
base recirculation region to 
the total enthalpy of the pri- 
mary flow. 

Y , Mass fraction of bleed flow gas none 1. 

^r contained in the base recircula- 
tion region. 

N, Nash factor none .4 

Convergence criteria for none .001 

Convergence criteria for n.j none .0002 

Wi , weighting factors used in testing none 3*1. 
the convergence of the varables sought 
by Newton's method. 

- 1 ? 

Convergence criteria for Newton's - 10 

method. Convergence occurs when 

n 2 

E, (WiFi ) < EPSI 
1=1 


n, number of varables to be determined. 

If n=3, then P./P , H /H , and Y 
b o_ r o„ p 

p p f'r 

will be determined. 

If n = l , then Pu/P„ will be determined. 

» % 

A1 ways set H /H = 1 , Y =1 
I” 0 _ p 

Perterbation factor for calculating par- - .02 

Itial derivitives for Newton's method 
^see write-up of subroutine NEWT) 
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NAMELIST EXPLANATION AND COMMENTS UNITS DEFAULT 

.... NA-ME y./yz!iL_ 


SK = k, proportionality factor relating shear none 4 

layer velocity thickness, 6, to the 
wake radius ratio, The pro- 

gram has been calibratea versus data 
for air without bleed flow and it was 
found that k = 4, i.e. it is assumed 
that 

r„/r^ = k (x/r,j)/<, 

where 

x/r is the non-dimensional 
^ shear layer length 

a is the jet spreading para- 
meter. 

0 = 1 2 + 2.758Mjj 

IDBUG Flag for debuging print-out - 0 

If IDBUG = 0, no debug print 
If IDBUG = 1, print debug output 


MAXIT = Maximum number of iterations allowed for - 10 

Newton's method. 

MAXITP = Maximum number of iterations allowed in - 25 

calculating Y,j 

MAXITE = Maximum number of iterations allowed in - 25 

calculating n- 
<3 

LAST = Flag to avoid error message that results - 0 

from reading end of file for input. For 
the last case to be run, set 

LAST = 1 

SEND 


247 





TABLE 1, CARD LISTING FOR A SAMPLE CASE 


** CALIBRATIUN CASE ** GAMMAxl.u Al« 

SDHALEX 

CPSr,2396«R7#CPP=. 2496897 r AMWSc29, , AMWPs29, , 
HUS=l,,HfJP=l,,POP=l,, 

SR8=1., 

THETAHzO, , 

AMP=2.3, 

ANSI., 

WSWPsO., 

PMAXs.079,PBPUPl=.0«, 

NUNS t , 

ANs.a, 

LASTsl, 

SEND 
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FINAL Rl.f'Oia l'l^r^llmTION 1 1ST 

DUAL NOZ7LC ACROUYNAMIC AND COOLINr. ANALYSIS 
CONTRACT NAS 8-33S5.T 


NAME 

cppi.es 

NAME 

National Aeronautics and Space Ad.iiinistration 


Rocketdyne 

Division of Rockwell International 

George C. Marshall Space Flipht Center 


6633 Canoga Avenue 

Marshall Space Flight Center, AL 35812 
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